On asymmetry of magnetic activity and plasma flow temperature in Jupiter’s magnetosphere

Discs of plasma around giant planets are natural laboratories that contain within mechanisms of transferring and keeping energy into the plasma and magnetic field system. Various missions to Jovian planets observed that expansion of plasmadiscs is not adiabatic and plasma temperature is increasing with radial distance. Magnetometer measurements from Juno mission were examined to determine plausibility of turbulent fluctuations as the plasma heating mechanism. Extensive azimuthal map of magnetic activity in Jupiter’s nightside plasmadisc is presented. Observations show that magnetic activity is distributed asymmetrically, with active and quiet regions. This is similar to the asymmetrical distribution of activity observed in Saturn’s magnetosphere. However, comprehensive study of temperature measurements showed that the only systematic change of temperature in magnetospheres of giant planets is in the radial direction. Observed breakfrequency in the magnetometer time series is systematically greater than the ion cyclotron frequency. Examination of the power spectrum points to that the kinetic energy of the corotating plasma as a source of increase of plasma temperature. This study shows that turbulent fluctuations themselves are not good candidates as a plasma heating mechanism. External pressure fluctuation however, can be used to convert kinetic energy of the plasma flow into thermal.

Volcanic activity on Jupiter's moon Io expels neutral gas at the rate of about 1 ton/s 1,2 .Part of the neutral population is ionized, forming a plasma torus consisting of oxygen and sulfur ions.Heavy ions are picked up by Jupiter's strong rotating magnetic field, causing plasma to flow around the planet 3 .Moving plasma expands radially forming the plasmadisc.The tilt of Jupiter's magnetic dipole causes magnetic field and plasmadisc to wobble.Here measurements and remote sensing of emissions of wobbling plasmadisc are described with respect to the centrifugal equator 4 .
Observations on various missions to the giant planet showed that Jupiter's plasmadisc expands nonadiabatically.Temperature of plasma increases with increase in the radial distance.Magnetosphere of a giant planet is a natural laboratory that contains within a mechanism of transferring and keeping energy into the plasma fluid and magnetic field system.One suggestion for the mechanism of an increase in plasma temperature in Jovian magnetospheres is heating of plasma fluid with turbulent oscillations.
Turbulent flow consists of eddies on different scales.Turbulent power is introduced to the system at the stirring scale 5 .Energy then cascades down from larger eddies at lower frequencies to smaller eddies with higher frequencies.The cascade is local and one way.Energy transfers from larger scales to neighboring smaller scales without skipping across the range.A major marker on the power spectrum of plasma fluid is the scale of gyration of ions in the magnetic field.At low frequency large fluctuations fluid dynamics is governed by magnetohydrodynamic (MHD) principles where ion fluid motion is associated with fluctuations of the magnetic field a "frozen in" condition.In that subrange power spectrum decays at Kolmogorov rate P(f ) ∝ f −5/36 .At frequencies higher than gyrofrequency, in the kinetic (KAW) subrange, variations of ion velocities and magnetic fluctuations decouple.In this subrange spectral index changes to roughly P(f ) ∝ f −7/3 , resulting in a steeper decay.Turbulent power cascades from mixing scale, through MHD subrange to gyrofrequency and then continues to cascade through the kinetic subrange.
In this study magnetic activity in Jupiter's nightside plasmadisc and power spectrum of magnetic fluctuations are examined to determine plausibility of turbulence as a heating source.Methodology of magnetometer time series measurements and systematic power spectrum analysis are presented in Section "Magnetometer measurements".An azimuthal map of magnetic activity in Jupiter's nightside magnetosphere is presented in Section "Local time distribution of magnetic activity".Examination of power spectrum features and plasma temperature estimate is presented in Section "Examination of the power spectrum of moving plasma and estimation of plasma temperature".Discussion of results is in Section "Discussion".

Magnetometer measurements
Time series of the magnetic field in radial Jupiter magnetic (VIP4) coordinates 7 measured by magnetometer instrument on Juno spacecraft 8 were analyzed in 10 min window intervals with 1 s resolution.This allows to observe plasma with similar properties within the sample and provides a consistent high temporal resolution of magnetometer time series.The main magnetic field B 0 (t) was calculated from moving average of the magnetometer time series 9 (blue line in Fig. 1a).The fluctuation of the magnetic field is then calculated as δB(t) = B(t) − B 0 (t) .Parallel perturbation of the magnetic field is found by δB(t) � = [δB(t) • n(t)] n(t) where n(t) is a unit vector in the direction of B 0 (t) .Perpendicular perturbation was then calculated as δB(t) ⊥ = δB(t) − δB(t) � .
Power spectrum is calculated from a time series of δB ⊥ components with Continuous Wave Transform using Morlet wavelet 10 .
Where W i (t i , f ) is a Morlet wavelet, t is the time series step size, f is a sampling frequency and N is a number of frequency samples.Total power density is then a square root of the sum of squares of the power spectrum of vector components.Perpendicular perturbation of the magnetic field in frequency space is calculated from the power spectrum as (δB ⊥ ) 2 = P(f )f 11 and then averaged across frequency range to get value for the sample window (Fig. 1c).Latitude off the centrifugal equator lat cent was calculated as in 4 and then distance from the centrifugal equator is calculated as Z cent = r sin(lat sysIII − lat cent ) .Where lat sysIII is latitude in planet's rotating frame 3 and the radial distance is denoted as r in Jupiter radii.
MHD and kinetic subranges are separated by scales of gyration of ions in the magnetic field.In this study gyrofrequency and frequency of the break in the power spectrum were found independently for each 10 min sample.Gyrofrequency is calculated as f g = qB 2πm i .Frequency of the break in the power spectrum was determined with the subroutine employing slope spectral index fit to the expected power law decay (Fig. 1b).Power spectrum was examined for a break in the interval f ∈ [10 −2 , 1.5 × 10 −1 ] Hz in 10 −2 Hz steps.Slopes were fitted to the power spectrum to determine subrange indices and then compared to the expected power law decay rate in the subrange ( −5/3 for MHD and −7/3 for KAW).MHD slope was fitted in the subrange 7 × 10 −3 Hz, f b /1.5 and KAW slope was fitted in the subrange f b × 1.5, 2 × 10 −1 Hz .Distribution of spectral indices in MHD and KAW subranges separated by the breakfrequency found in the power spectrum samples is shown in Fig. 2. Some spread in slopes shows that simply not all fluctuations in the magnetosphere are due to turbulent cascade.
Differences between observed frequencies of the break in the power spectrum and ion cyclotron frequencies in 10 min samples are shown in Fig. 3.There are a few samples, where found breakfrequency was lower than the gyrofrequency.However, in the vast majority of cases breaks in the power spectrum of moving plasma in the Jupiter magnetosphere are observed at frequencies that are higher than ion cyclotron frequencies.Similar phenomenon was observed in the solar wind 12 , where breakfrequency was found to vary between f g < f b < 10f g .In cases when breakfrequency was not found by the subroutine it was estimated by the median value f b = 5f g (not included in the Fig. 3).The systematic difference between the scales of gyration of ions and scales of observed break in the power spectrum is examined in more detail in Section "Examination of the power spectrum of moving plasma and estimation of plasma temperature".
It is important to note that often in literature 13,14 (and let's not forget motivation for this study 15 ) gyrofrequency is calculated as ω g = qB m i .This however is a calculation of angular velocity which has units of rad/s and should not be plotted on the Hz axis.Interestingly, unitless 2πf g ∼ 6f g in Jupiter's magnetosphere would put this in the expectation range in Fig. 3.So this then estimates a break in the power spectrum in Jupiter's magnetosphere in the right ballpark for a wrong reason.

Local time distribution of magnetic activity
Equatorial map of magnetic activity of the night side of Jupiter's magnetosphere is shown in Fig. 4. Map was assembled from perpendicular fluctuation values measured in 10 min time series samples (Fig. 1c).Perpendicular fluctuation values were averaged over vertical distance Z cent ∈ ±10 .Figure 4 demonstrates an asymmetric azimuthal distribution of magnetic activity in the nightside side plasmadisc.Dusk is more active than dawn.Observed part of the plasmadisc was divided into two local time (LT) regions.Amplitudes of magnetic fluctuations were observed with respect to the azimuthal bend of the magnetic field.
Dusk to midnight region LT ∈ [−3, 3] hrs is more active with higher amplitude of perpendicular magnetic fluctuations with δB ⊥ ∼ 10 −9.75 T (Fig. 4b).Natural topology of the magnetic field inside the plasmadisc is an azimuthal bend back configuration.The observed bend of the magnetic field of the majority of cases in that region is close to neutral.Active magnetic field would tend to oscillate between natural bend back and bend forward configurations.So that the chance observation of the bend of the active magnetic field would probably happen somewhere between the extreme forward and backward deflections.Dawn region LT ∈ [3, 5] hrs is much quieter with majority of cases near δB ⊥ ∼ 10 −10.5 T with negative bend of the magnetic field (Fig. 4c).The spread of fluctuation magnitude values in the active region is about a half an order of magnitude, while values in the quiet region are tightly packed.This is characteristic of the turbulent environment.Turbulent formations exist on multiple scales where energy cascades from larger scales to smaller ones.So that increase in turbulence not only increases the amplitude of fluctuations but also increases the spread of amplitude ranges.This points to that the activity in dusk to midnight region is indeed turbulent.This result is similar to what was observed in Saturn's magnetosphere 15 .Survey of magnetic fluctuations quantified as turbulent heating indicated presence of a quiet LT ∈ [3, 9] hrs and an active LT ∈ [10, 20] hrs regions in Saturn plasmadisc.The dusk-dawn asymmetry of magnetic activity in magnetospheres of Jovian planets must be a result of the interaction of the spinning magnetosphere and the solar wind.This implies that this interaction is a major contributor to turbulence in magnetospheres of Jovian planets.Systematic study of magnetic fluctuations around giant planets 15,16 surprisingly point to that the magnetic activity is higher in azimuthal regions where the flow in the sheath opposes the natural backward azimuthal bend of the magnetic field rather than the flow in the plasmadisc.It is also surprising to see the effect of interaction of the magnetosphere with solar wind well within r < 35 R J .
There is no discernible dependence of magnetic activity on local time beyond r > 35 R J within azimuthal regions observed by Juno spacecraft.This study also shows that there is no System III dependence of magnetic fluctuations in the plasmadisc.This differs from results published in 17 .Additional description and figures are in the appendix of this paper.
Extensive examination of moments from the Galileo mission showed that there is no significant local time or System III variation of ion temperature in Jupiter's plasmadisc and that the only systematic variation is in the radial direction 1 .Similarly no such LT asymmetry of ion temperature in Saturn's magnetosphere was reported 18 .Which points to that asymmetrically distributed turbulent activity in plasmadisc is not a good candidate for the heating mechanism of plasma in Jovian magnetospheres.

Examination of the power spectrum of moving plasma and estimation of plasma temperature Doppler effect on observed break in the power spectrum
Observation of magnetometer time series showed that break in the power spectrum systematically occurs at frequencies that are higher than the gyrofrequency (Figs. 1 and 3).Fluctuations in a moving medium are described in terms of oscillations as observed in the frame moving with the medium modified by a Doppler shift.In Jupiter's magnetosphere Doppler shift effect can be written as.
(1) Where ⊥ implies quantity perpendicular to the main magnetic field.Here v ⊥ = v sin(θ vB ) 13,19 and θ vB is the angle between plasma bulk flow velocity and main magnetic field in Jupiter sun orbit (JSO) coordinate system 7 .
Here note that measurement of δB is a measurement of the amplitude of the oscillation and not wavelength or direction of the wave propagation.This simplified estimate of the dot product assumes that waves detected by the spacecraft are propagating along the direction of the perpendicular component of the bulk velocity, moving towards the spacecraft.It is reasonable to assume that the wavelength of the perpendicular fluctuation of the field at ion gyrofrequency is comparable to the circumference of ion gyration k ⊥ (f g )ρ i ∼ 1 (comparable cycles of motion of ions and magnetic field).In this work it is deduced that Doppler shift is responsible for the difference between observed frequency of the break and gyrofrequency (Fig. 1).Then Eq. ( 1) can be written in terms of breakfrequency, gyrofrequency and gyroradius.
The ion gyroradius is calculated in terms of the ion temperature.
Introducing Doppler shift constant equation ( 3) can be rewritten as Substituting for gyrofrequency and gyroradius it can then be solved for the ion temperature Note that in Eq. ( 2) a bulk flow in Jupiter's plasmadisc was assumed to be faster than the wave speed.So that fluctuations of waves moving towards the sensor are observed.In more general situation, the Doppler shift k • v can very well be negative.In that case break in the power spectrum will be observed at frequencies lower than the gyrofrequency f b < f g .In principle, magnitude of the angle between bulk velocity and wave propagation can be anywhere between 0 and π .If direction to the source of the wave is known then the more general expression for the ion temperature is As the angle between k and v approaches π/2 , f b will approach f g so that both c ds and nk • nv will become small with canceling effect and as the angle between k and v approaches π both c ds and nk • nv will become negative.Also note that c ds > −1∀f b and f g > 0.
In absence of a better estimate for nk • nv temperature was calculated using Eq. ( 4) with cases in the range f b > 2f g .Here bulk velocity modeled as a corotation velocity up to 20 R J , After that velocity is v = v φ = 200ê φ km/s 18 .Near the planet gyro and breakfrequencies for cases described above increase with increase in the strength of the magnetic field.In samples from magnetometer measurements 8 where break is beyond the upper bound of the range used by the algorithm 1.5 × 10 −1 Hz (Fig. 1b) or where break in power spectrum was not found, breakfrequency is estimated by the median value (Fig. 3).
Figure 5 shows distribution of temperature estimates.Temperature calculations made with breakfrequency found by the subroutine are shown with blue markers.These estimates compare fairly well to the plasma temperature fit from Galileo mission 18 with the spread that is similar to the spread in reported temperature measurements.Temperature calculations made using the median estimate for the breakfrequency are shown with red markers.These estimates land on the Galileo profile.So that if one would not bother to try to find a break in the power spectrum and just calculated it as 5f g then Eq. (4) will give the temperature profile from 18 .
An intriguing possibility is the use of external fluctuations to trigger cascade in moving plasma closer to gyrofrequency.In this case both c ds and v b in Eq. (4) will decrease.This in principle could allow to convert kinetic energy of the plasma flow into thermal.It is possible that waves generated by magnetic activity in Jupiter's magnetosphere are used to precipitate a break in the power spectrum closer to gyrofrequency.The amount of kinetic energy transferred to thermal, then depends on the difference between gyro and breakfrequencies.In Jupiter's magnetosphere this process is inefficient, the spread of the difference between break and gyrofrequencies is about an order of magnitude (Fig. 3).This then results in the natural spread of temperature values 18 and Fig. 5.

Doppler effect on general frequency spectra
Assuming that the Doppler effect is similar for fluctuations at different frequencies a more universal relation between observed wave numbers and wave numbers in the bulk flow reference frame can be made.In that case a relation between observed and original frequencies is similar to the relation between break and gyro frequencies.
Here the observed frequency of the break in the power spectrum is used as a marker for the Doppler shift of observed wave fluctuations in the rest of the frequency range.So then Eq. ( 1) can be solved for k.
Note that Eq. ( 5) will deviate from Taylor's hypothesis 20 when observed frequency is very similar to the original frequency and c ds ≪ 1 , so in the case of a very low velocity or a very long wavelength.It allows for the observation of waves with an arbitrary angle between directions of wave and bulk velocities.So that this enables an observation of incoming wave fluctuations from observed frequency, bulk velocity and the breakfrequency in the power spectrum.Alternatively, if wavenumber features k(f) in rest frame are known then using observed f ′ and f b Eq. ( 5) can be solved for direction to the source nk • nv .

Discussion
In this work Juno magnetometer measurements were used to examine magnetic activity in Jupiter's night side magnetosphere and to attempt to determine mechanism responsible for increase of plasma temperature in Jupter's plasmadisc.Magnetic field observations show an asymmetry of activity between night: dusk and dawn sectors of the magnetosphere.The dusk side is more active than the dawn side.Amplitudes of fluctuations in the active dusk region are about two orders of magnitude greater than in the dawn region.Local time asymmetry of magnetic activity must be a result of interaction of solar wind with magnetosphere of the planet.Surprisingly, observation of magnetic fluctuations show that magnetic activity is higher in azimuthal regions where the flow in the sheath opposes the natural backward azimuthal bend of the magnetic field rather than the flow in the plasmadisc.Equally surprising result is to see the effect of interaction of the solar wind with magnetosphere well within r < 35 R J .Similar azimuthally asymmetric distribution of magnetic activity was observed in Saturn's magnetosphere 15,16 .Survey of fluctuations of the magnetic field from Cassini magnetometer data, quantified as turbulent heating, showed local time dependence of activity in Saturn's magnetosphere 15 with active and quiet sites.Azimuthal locations where flow in the sheath is opposite to the flow in the plasmadisc have lower magnetic activity than areas where flow in the sheath and disc are in the same direction.This is opposite to expected Kelvin-Helmholtz instability driven turbulence process.So then Kelvin-Helmholtz vortexes had to be transported to the other side of the giant's magnetosphere 16 .Examination of magnetic activity in Jovian magnetospheres point to that it is more likely that magnetic fluctuations due to field topology in conjunction with fluid velocity variations play a crucial role in generating turbulent oscillations.
Turbulence was suggested as a possible source of energy needed to heat plasma in Jovian magnetospheres.Estimates of the total power input necessary to counteract an adiabatic cooling of expanding plasmadiscs in Saturn's and Jupiter's magnetospheres were done in 18 .However, on closer examination it is unclear what volume was used in calculation of the adiabatic expansion in 18 .There an attempt was made to scale the volume of a shell with change in the magnetic flux at different L dipole McIlwain parameters.Volume was calculated as L 3 H .This volume however does not scale as a disc where plasma fluid is contained.Indeed R 3.85 H in equation ( 13) of that publication is not a 3D volume at all, as is required for adiabatic process relation derived from thermodynamic equations for plasma fluid contained in Jupiter's plasmadisc or coffee contained in a thermos.
If turbulent power is responsible for the increase in plasma temperature and the timescale of transfer of turbulent energy in to ion temperature "heating" is less than the rotation period of the plasmadisc ∼ 10 hrs , then one would expect to see a systematic increase in plasma temperature as it passes active regions.However, extensive examination of moments from Galileo showed that there is no significant local time or System III variation of ion temperature in Jupiter's plasmadisc 1 .Similarly no such local time asymmetry of ion temperature in Saturn's magnetosphere was reported 18 .Which points to that azimuthally asymmetric turbulent activity in plasmadiscs of Jovian planets is not a good candidate for the mechanism to rationalize a radial increase in plasma temperature.
Magnetohydrodynamic and kinetic subranges are separated by scales of gyration of ions in the magnetic field.It is reasonable to assume that the wavelength of the perpendicular fluctuation of the field at ion gyrofrequency is comparable to the circumference of ion gyration.Examination of power spectrum of magnetometer time series samples shows that break in the power spectrum of moving plasma consistently occurs at frequencies that are on average half an order of magnitude higher than the ion cyclotron frequency.In this work it was deduced that this effect is due to the Doppler shift of observed oscillations.Doppler effect relation can then be used in conjunction with comparison of field wavenumber and inverse of ion gyroradius to express plasma temperature in terms of bulk flow velocity.This then implies that kinetic energy of the plasma flow is the primary source of increase of ion temperature.Turbulent fluctuations themselves are not good candidates for a plasma heating mechanism.External fluctuations however, can be used to convert kinetic energy of the plasma flow into thermal.
Difference between observed frequency of the break in the power spectrum and ion gyrofrequency can be used as a marker for the Doppler shift for the rest of the frequency spectrum.This then enables to modify Taylor's hypothesis to take into account oscillations of passing wave features.Wavenumber spectrum of the incoming wave in the bulk flow reference frame can be inferred from the observed frequency, bulk velocity and breakfrequency.

2 Figure 1 .
Figure 1.(a) Magnetometer time series of the 10 min sample window.Blue line is the 200 s centered moving average.Area outside dashed lines is a buffer for a moving average sliding window.(b) Power spectrum of perpendicular fluctuations.(c) Perpendicular magnetic field fluctuation in frequency space.

Figure 2 .
Figure 2. Powerlaw spectral index distribution in MHD and KAW subranges separated by f b using a slope fitting subroutine.Samples are taken in r ∈ [5, 40] R J , Z cent ∈ ±1 R J in all sampled local time.

Figure 3 .
Figure 3.Comparison of gyro and break frequencies in samples from Fig. 2.

Figure 4 .
Figure 4. (a) Equatorial map of perpendicular magnetic fluctuations.Values of δB ⊥ are averaged over Z cent ∈ [−10, 10] R J .Dashed line is r = 35 R J .Color bar shows an order of magnitude. (b) 2D histogram of magnetic fluctuation vs azimuthal bend of the magnetic field for the region r ∈ [5, 35] R J , LT ∈ [-3, 3] hrs .(c) 2D histogram of magnetic fluctuation vs azimuthal bend of the magnetic field for the region r ∈ [5, 35] R J , LT ∈ [3, 5] hrs .Azimuthal bend of the field was calculated as a proportion of the B φ component of the magnetic field with respect to the total strength of the field and multiplied by a sign of the B r component.Colorbar shows number of sample counts. https://doi.org/10.1038/s41598-023-41500-ywww.nature.com/scientificreports/

Figure 5 .
Figure 5. Radial profile of Temperature estimate made with breakfrequency found by the subroutine of ion temperature (blue markers).Temperature estimate made with breakfrequency calculated as f b = 5f g (red markers).Samples are taken in r ∈ [5, 40] R J , Z cent ∈ ±1 R J in all sampled local time.Temperature fit from Galileo mission 18 (dashed line).