1. Introduction
One of the most interesting features of the equatorial ionospheric F region is the existence of plasma density irregularities, which find the suitable conditions for their development during post-sunset hours [
1]. In the last few years, different terminologies have been used to term the equatorial ionospheric irregularities. Often, these are identified as equatorial spread-F, in virtue of the spread of the experimental trace either in frequency or in amplitude affecting the ionosonde observations, which is the marking of the existence of electron density irregularities in the reflecting layer (e.g., Booker and Wells [
2]). With the advent of in situ satellite observations, they started using terms such as “depletions”, “bite-outs,” or “plasma holes”, because the irregularities were recorded mostly at locations where the electron density is depleted with respect to the background ionosphere (e.g., McClure et al. [
3], Dyson and Benson [
4]). Then, terms such as “plumes” or “wedges” appeared to characterize the morphology of turbulent regions and the development of such irregularities, as observed by radars (e.g., Scannapieco and Ossakow [
5], Woodman and La Hoz [
6]). In the last few decades, to identify the generation process of these irregularities, the term “bubble” has taken hold more and more (e.g., Burke [
7]). It is worth noting that recent studies (e.g., Kil et al. [
8]) introduced also the expression “plasma depletion shell” to highlight the three-dimensional structure of bubbles. The formation of these irregularities is recognized to be driven by the Rayleigh–Taylor instability mechanism [
6,
9,
10] that generates in the ionosphere a situation similar to that occurring when a heavy fluid flows over a lighter one. The absence of sunlight, which leads to a different rate of recombination in the ionospheric layers, causes a density gradient between the upper and lower ionospheric F region with a plasma density in the upper F region higher than in the lower one. Such steep vertical density gradients realize the conditions for the growth of the Rayleigh–Taylor instability. These plasma density irregularities expand vertically and successively move along the magnetic field lines on each side of the geomagnetic equator. Typically, they also drifted eastward from vertical polarization electric fields associated with zonal neutral winds in the F region.
The result is the formation of magnetic field-aligned plasma-depleted regions, which are usually known as equatorial plasma bubbles. They are unstable to plasma density perturbations and polarization electric fields so that secondary irregularities can be generated, which are mainly the result of cascading processes. Nowadays, many important characteristics of these irregularities are known. They have different scale sizes, from hundreds of kilometers down to a few decameters [
11,
12], are mainly observed within a narrow band of
magnetic latitude [
9], are observable in a wide range of altitudes [
6,
9] up to around
km [
13,
14,
15], and their spatial and temporal distributions are strongly controlled by various geophysical parameters such as solar [
16] and geomagnetic activity [
17], season [
16,
18], local time, and longitude [
8,
19,
20,
21]. What drives the scientific community to investigate the plasma density irregularities is the knowledge of their role in the propagation of electromagnetic signals. They can be responsible for the disturbance of radio waves and are one of the most important sources of disturbance for the Global Navigation Satellite System (GNSS) signal. Indeed, in the worst case, plasma density irregularities can cause a loss of lock event, a condition for which a GNSS receiver can no longer track the signal sent by the satellite, with a consequent degradation of the positioning accuracy. In the last few years, the relevance of the GNSS in our society has substantially increased, as many critical infrastructures and economies are dependent on the positioning, navigation, and timing services of GNSS, so that our society is nowadays vulnerable to damages due to the malfunction of these systems. For this reason, one of the research priorities in the space weather community is to improve the knowledge of these plasma density irregularities, trying to better understand their features and generation mechanisms.
An interesting aspect of these irregularities, which requires more study, is their possible turbulent nature. Indeed, the equatorial ionosphere has proper physical conditions for the establishment of a fully developed turbulence [
22,
23]. Through numerical simulations, it has been possible to reproduce the ionospheric irregularities and study the corresponding dynamics under collisional and inertial regime flows [
24,
25]. It has been found that the plasma irregularities may be characterized by energy spectra with slopes ranging between
and
associated with the occurrence of energy and/or enstrophy cascades, respectively [
26,
27]. Over the years, these numerical models have been improved trying to describe the ionospheric two-dimensional plasma turbulence in a more realistic way. It has been found that after sunset, at apex heights between 600 and 1000 km, large-scale irregularities present a dynamics typical of inertial regime flow [
28] and that although the foundations of the Kraichnan’s theory could be undetermined in the topside equatorial ionospheric F region, the numerical simulation results seem to confirm the probable existence of large plasma density irregularities. These are characterized by velocity fluctuation spectra following power laws with spectral indices similar to those predicted for energy and enstrophy cascades in two-dimensional turbulent flows [
23]. In the last twenty years, thanks essentially to the rapid increase in the performance of computers, it has been possible to improve the models describing the evolution of plasma depletions. Thus, new models capable of reproducing the three-dimensional turbulent structures of equatorial plasma bubbles with high spatial and temporal resolutions have been developed (see for example Yokoyama et al. [
29], Yokoyama [
30]). The results of such simulations confirm observations of plasma density data recorded, for example, on board DE-2 satellite. In this case, it has been found that plasma density is characterized by an energy spectrum with a slope around
for scale sizes larger than 1 km in all regions characterized by strong plasma density irregularities [
28]. Nowadays, the numerical simulations have become useful tools both to study the plasma bubbles and their nonlinear evolution, and to understand some of their features, which cannot be fully understood from theoretical predictions. Therefore, it is extremely important to combine the numerical simulation results with real observations.
In order to detect and characterize plasma density irregularities, different methods, approaches, instruments, and data can be used such as, for example, ionosondes, airglow imager, coherent/incoherent scatter radars [
31], as well as Global Positioning System (GPS) Total Electron Content (TEC) [
32,
33] in the topside ionosphere. The ground-based observations are mostly fixed in a single location and cannot give the spatial variations of plasma density irregularities in a large area. Conversely, observations from Low Earth Orbit (LEO) satellites such as for example DMSP (Defense Meteorological Satellite Program), ROCSAT-1 (Republic of China Satellite), C/NOFS (Communications/Navigation Outage Forecasting System), and COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate) permitted climatology studies and provided the latitude-local time distributions of plasma density irregularities. The recently launched Global-scale Observation of Limb and Disk (GOLD) mission, which is located at a geo-stationary orbit, seems to be a promising mission to study bubbles, offering the opportunity to observe the Earth’s complete disk continuously from the geostationary orbit [
34,
35]. A different method for the detection of ionospheric irregularities based on the diamagnetic effect was suggested a few years ago. Plasma density irregularities are indeed associated with magnetic field perturbations [
36,
37] that can be used as a proxy for plasma irregularities. There are different ways to describe the so-called diamagnetic effect but, perhaps, the easiest one is to consider that the sum of magnetic and plasma pressures has to be constant in a quasi-stationary state. This means that when the plasma pressure is reduced, the magnetic pressure and consequently the magnetic field strength must increase. The equatorial plasma bubbles identify regions of reduced plasma density; consequently, they have to be characterized by a magnetic field strength higher than the ambient one. That was confirmed using simultaneous observations of magnetic field and electron density recorded by the CHAllenging Minisatellite Payload (CHAMP) satellite. These measurements revealed the increase of a few nT in the main field intensity inside the plasma density depletions [
38] Park et al. [
39] quantitatively investigated the actual balance between plasma and magnetic pressures within the plasma bubbles, showing that a dominant part of the magnetic pressure change could be explained in term of plasma density change under specific assumptions. Based on the high correlation between the variations of plasma density and those associated with magnetic field strength, an extensive survey of magnetic signatures related to plasma density irregularities was used for statistical studies [
15,
40,
41] devoted to the analysis of the global distribution of equatorial plasma irregularities in the topside ionosphere. These studies revealed features of the magnetic signatures that closely reflected those of the plasma bubbles previously obtained using different methods as GPS scintillation [
19], in situ plasma density measurements [
42,
43], and radio wave propagation [
44]. Furthermore, magnetic field measurements recorded on board the CHAMP satellite with different time resolution (50 Hz and 1 Hz) permitted addressing some properties of the plasma density irregularities related to their different spatial sizes. Evidently, these results do not reflect directly the features of plasma depletions but their associated magnetic signatures. However, the features obtained using the magnetic signatures of the plasma depletions are consistent with those obtained by directly using the electron density measurements. For this reason, the diamagnetic effect, as an indirect way of sampling plasma density structures, can be extremely useful in studying plasma irregularities.
While it is correct to assume that some plasma bubbles features can be inferred using the corresponding magnetic signatures, the assumption that plasma bubbles scaling properties can be inferred using the associated magnetic field scaling properties may be not necessarily correct as well. The purpose of this study is to understand if the hypothesis of a linear relation between the scaling properties of the electron density and of the magnetic field is valid. So, needless to say, the clarification of this point is of crucial importance.
2. Data Description
To investigate the scaling features of the electron density irregularities and of their magnetic signatures at low latitude, we use measurements on board the Swarm A satellite during its near-polar orbit around the Earth at an altitude of approximately 460 km [
45]. Swarm A is one of the satellites constellation launched by ESA for Earth observation. Specifically, the ESA mission consists of three identical satellites named Alpha, Bravo, and Charlie (A, B, and C), which were launched on 22 November 2013 into a near-polar orbit. Swarm A and C fly side-by-side at an altitude of 462 km (initial altitude) and at an inclination angle of 87.35°, whereas Swarm B flies at a higher orbit of 511 km (initial altitude) and at an inclination angle of 87.75°.
We consider measurements collected from 1 April 2014 to 31 January 2016 at low and mid latitudes (
) under conditions of high solar activity. It is known indeed that the equatorial plasma bubbles occur mainly during periods characterized by high extreme ultraviolet solar flux [
41]. In the selected time interval, the value of the solar radio flux index (F10.7) is equal to
sfu. This is not an extremely high value for the solar activity, but the selected period, which covers a part of the descending phase of Solar Cycle 24, is the only period of moderate/high solar activity during which Swarm is flying. Electron density and magnetic field vector data at a rate of 1 Hz are selected from the ESA dissemination server (
http://earth.esa.int/swarm accessed on 12 February 2022). To highlight the magnetic signatures associated with plasma bubbles, we remove from the magnetic field measurements the contribution due to the magnetic field of internal origin, as modeled by CHAOS-6 [
46]. CHAOS-6 is the latest generation of the CHAOS series of global high-resolution geomagnetic field models, which spans a period between 1997 and 2016. It is derived primarily from magnetic satellite data (Swarm, CHAMP, Ørsted, and SAC-C), although ground-based activity indices and observatory monthly means are also used. It is able to model the dominating core field, the crustal field, and the magnetic fields from large-scale magnetospheric currents, but there is no explicit representation of the ionospheric field or fields due to magnetosphere—ionosphere coupling currents. One of the applications of this model is the removal from the observed magnetic field of the large contributions from the prime sources when the weaker effects such as ionospheric currents, plasma irregularities, or even ocean tidal motions want to be investigated. We use the CHAOS-6 package to compute the geomagnetic field of internal origin, which is successively subtracted to the magnetic field observed by Swarm. Lastly, we focus on the 18:00–24:00 local time (LT) sector and select periods characterized by low/moderate geomagnetic activity. To select these periods, we used the
index, which is designed to measure the solar particle radiation through its magnetic effects, and today, it is considered a proxy for the energy input from the solar wind to Earth. It is a global geomagnetic activity index based on 3-hour measurements from ground-based magnetometers around the world. The
index ranges from 0 to 9, where a value of 0 means that the geomagnetic activity is absent, and a value of 9 means that an extreme geomagnetic storming is occurring. Data are distributed by GFZ (Helmholtz Centre Potsdam) [
47] and they are redistributed by various data centers and databases. To select periods characterized by low/moderate geomagnetic activity, we consider all those time intervals characterized by
. According to previous results [
40], the probability of observing plasma bubbles is maximum under these conditions.
To select the plasma bubbles events, we use the Swarm Level-2 Ionospheric Bubble Index (IBI) product, which was evaluated using plasma density and magnetic field data with a cadence of 1 Hz. A full description of the data file and explanations regarding the IBI data can be found in the L2-IBI product description file (
https://earth.esa.int/eogateway/documents/20142/37627/Swarm-Level-2-IBI-product-description.pdf/3e9f6c3a-1ffc-ea53-0161-a18b63f90c6f, accessed on 13 February 2022). This product is composed of three parameters: Bubble Index, Bubble Probability, and Bubble Flag. The Bubble Index can assume three different values: 1, when the analyzed data are affected by a plasma bubble, 0 when there is not a plasma bubble, and −1 when conditions do not allow unequivocally identifying a plasma bubble in the data [
48]. Furthermore, when IBI = 1, there is an additional flag, the Bubble Flag (BF), indicating the quality of plasma bubbles detection that is related to exceeding a certain probability threshold. When BF = 1, we are in the presence of a high correlation between plasma density depletion and magnetic field signature. Here, we chose the condition IBI = 1 and BF = 1 to identify the plasma bubbles. However, we note that this automatic detection algorithm cannot distinguish equatorial plasma bubbles from plasma blobs, which are localized plasma density enhancements. This means that the events identified by the L2-IBI processor can contain also blob events. However, the occurrence probability of equatorial plasma bubbles is significantly higher than that of blobs mainly during solar maximum years.
3. Method and Data Analysis
The scaling features of the electron density and magnetic field at different temporal scales are evaluated using the so-called generalized
-order structure function, namely
. For a signal
,
is given by the following equation:
where
t is the time,
is the time delay, and
stands for a statistical average. When
has a power-law behavior as a function of
:
the analyzed signal is characterized by a scale invariance, and
is the so-called
order scaling exponent, which defines the scaling nature of the increments of the signal
. A signal possesses a global scale-invariance if
; that is, it is a linear function of the order
q. Conversely, a signal is defined as multifractal, that is, it has a local scale-invariance, if
is a nonlinear convex function of
q. Here, we investigate the second-order scaling exponent,
, which is related to the power spectral density (PSD) exponent,
(PSD
), via the following relation:
The use of the structure function analysis, rather than the standard Fourier analysis, as a method to detect scaling features is motivated by the fact that the analysis of scaling/spectral features in the real time space, using signal increments, is less affected by the choice of a specific decomposition basis, as it is instead the case for the Fourier analysis. Furthermore, being geophysical signals, such as magnetic field and electron density fluctuations in the ionosphere, generally characterized by power-law spectra, i.e.,
, with spectral exponents
ranging between 1 and 3, one deals with nonstationary signals with stationary increments [
49]. In this case, the Fourier analysis, which requires stationarity, might produce unreliable results especially when evaluating spectral properties over short time intervals. Conversely, the structure function analysis works with signal increments that are expected to be nearly stationary, and thus, the corresponding results over short time windows/intervals are more reliable. Furthermore, the structure function analysis is a companion method that is widely used in the framework of turbulence-related studies.
The measurements collected by instruments on board a satellite require some precautions to be correctly used. The satellite crosses regions characterized by different physical processes, which change in space and time. Consequently, the collected measurements, such as in our case electron density and magnetic field strength, present properties with a local character. The scaling properties of these real-time series also acquire a local character, and it is necessary to apply methods of analysis capable of evaluating their local scaling features describing the different regions explored by the satellite. The structure function approach allows doing this. Among the different methods to evaluate the structure functions, here, we apply the local detrended structure function analysis (DSFA) method discussed in De Michelis et al. [
50,
51] and already successfully applied in a series of previous works showing results that are well in agreement with the present literature on the scaling and/or spectral features of turbulent fluctuations in the ionosphere. In detail, we apply the analysis on a moving window of
points using a time series obtained from the original one after removing its linear trend, if present. This gives us the opportunity to remove from the initial signal all the possible large-scale variations that can affect a correct estimation of the scaling features. The choice of a moving window of 301 points is the result of an accurate series of test over a wide set of simulated multifractional Brownian motions [
52,
53,
54], which are characterized by a local dependence (nonstationarity) of the scaling/spectral features. Indeed, this number of points is the right compromise to get an estimation of local scaling exponents sufficiently precise (with a typical error
%), being 1 order of magnitude larger than the longest investigated timescale (
s). In each moving window, we evaluate the second-order structure functions associated with the electron density and magnetic field strength time series, respectively. These are evaluated for time delay values (
) ranging from 1 to 40 s, and the estimated second-order scaling exponents are associated with the position of the satellite at the center of the fixed window. We assume that the results obtained in the time domain are valid in the space domain too (see e.g., [
55,
56] and references therein). Thus, considering that the Swarm A orbital velocity is around 8 km/s, the results obtained in the time domain between 1 and 30 s are valid in the spatial domain between ∼8 and ∼250 km, i.e.,
where
is the satellite velocity. We remark that the previous conversion between the timescale and the corresponding spatial scale is based on the assumption of the Taylor’s hypothesis according to which instantaneous spatial averages coincide with temporal averages calculated from the recorded signal. Given this hypothesis as valid, which can be applied for frozen turbulence, we assume that the temporal evolution of turbulent structures is longer than that required for the satellite to cross them [
22,
55,
56,
57]. Hence, by evaluating the scaling properties of the magnetic field strength and electronic density associated with the plasma bubbles, we can investigate the properties of these plasma density irregularities from both a magnetic and plasma point of view in a range of scales between some kilometers and a few hundred kilometers.
The structure function analysis is initially done on the entire time series of around points (22 months with a time resolution of 1 Hz) and only later, the scaling exponent values of the electron density () and magnetic field strength () associated with the occurrence of plasma bubbles are selected.
4. Results
Figure 1 reports a map of the values of the second-order structure function exponent,
, which was obtained analyzing
and
inside the plasma bubbles. These two different parameters are characterized by similar values of
, supporting the previous findings by Lühr et al. [
15]. However, to look for any difference, we start by plotting in
Figure 2 the corresponding histograms. The second-order structure function exponents range from 0.4 to 2 with the highest probability occurring at
in both cases.
Taking into account the relation between the second-order scaling exponent and the Fourier power spectral density exponent (
), we notice that the spectral features obtained in the case of electron density (
) are in agreement with previous results. For example, Livingston et al. [
58] analyzing the spectra of electron density data recorded from the Atmospheric Explorer-E satellite (AE-E) in the interval 20:00–02:00 LT found a distribution of spectral slopes in the range between 1 and 3 centered around 1.9, which is a mean value that agrees with our results and is close to those of Dyson et al. [
59], who used data recorded by the OGO satellite, Basu et al. [
60], who analyzed AE-E measurements, and Le et al. [
61], who analyzed ROCSAT-1 burst mode measurements at 1024 Hz sampling. Concerning the spectral scaling exponents obtained for the magnetic field strength, there are not many previous analyses with which to compare them. Lühr et al. [
15] used magnetic field data recorded by CHAMP to investigate the spectral properties of the signal and found values of the spectral indices ranging between 1.4 and 2.6 with a peak around 1.9. Again, these results agree with ours. Although the mean values of the two distributions reported in
Figure 2 are equal, their shapes are different, in which the
distribution for
is skewed. Indeed, the skewness
values of
distributions of
and
are
and
, respectively.
Figure 3 reports the values of the second-order structure function exponent obtained analyzing
and
inside the plasma bubbles according to QD magnetic latitude. For
, the highest occurrence of
is found for
at a latitudinal distance from the equator of around
in magnetic latitude in both hemispheres, while in the case of
, the highest occurrence of
is found for a value around 1 occurring always in the same regions. Thus, the peak of the highest occurrence of plasma bubbles is associated with a range of
values that is smaller for
than
.
To go into the observed differences between the two distributions, we repeat our analysis considering two-hour time intervals from 18:00 LT to 24:00 LT.
Figure 4 reports the distribution of
values relative to
and
in the three different selected time intervals. A double peak distribution of
clearly appears analyzing the scaling properties of
in the range 20:00 ≤ LT < 22:00. The
distribution has one peak in the next local time interval (22:00 ≤ LT ≤ 24:00) and a rather flat shape in the previous one (18:00 ≤ LT < 20:00). In the case of
values associated with
, the data distribution has only one peak in the second and third local time intervals while a distribution with a rather flat shape appears in the first local time interval. The emergence of a bimodal distribution for the magnetic field
exponent could be an indication of a non-unique relation between the scaling features of electron density and magnetic field fluctuations in plasma bubbles. Furthermore, as clearly shown in several works (see, e.g., [
62]) (and references therein), in the second LT interval here considered, plasma bubbles are still evolving at Swarm altitude, so that the bimodality could be a reflection of the nonstationarity of the plasma bubble structures.
To get more insights on the bimodal character of the
distribution that appears, analyzing the scaling properties of
inside the plasma bubbles between 20:00 LT and 22:00 LT, we check how the
values of the two peaks are distributed in latitude and longitude. Taking into account that the two different populations are characterized by values of
and
, respectively,
Figure 5 shows their distributions in the latitude–longitude plane. The two geographic distributions seem similar, although the values
seem to prefer the areas closest to the equator. This feature is well highlighted on the bottom of
Figure 5, where the
values are sorted according to magnetic latitude. Here, the red line clearly shows how the plasma bubbles around the magnetic equator (
) are often characterized by a magnetic field strength with a second-order scaling exponent
. The values of the scaling exponents of the magnetic field strength seem to depend on magnetic latitude at least during the hours following their formation. To better investigate this point, we report in the left panels of
Figure 6 the joint probability distributions of
values relative to
(
) and
(
) along with the mean values of
exponent (white circles) for fixed
at the three different selected two-hour time intervals. There is a large spreading of the observed distributions with most of the values concentrated around
, which indicates that the relation between electron density and magnetic field scaling exponents can be assumed to be linear only in a zero-order approximation. However, looking at the trend of the mean values of
fixed
, we can realize how the relation between the two scaling indices significantly departs from a linear trend for
at LT
. This result suggests that a more complex relation might exist between plasma density and magnetic field variations. Furthermore, the departure from a linear relation increases with the midnight approaching. We can conjecture that the observed departure from linearity occurs when the fluctuations of the electron density show scaling features similar to those expected for convective turbulence, i.e.,
[
63]. In the right panels of
Figure 6, we report the conditioned probably distribution density functions of
values for fixed
. The observed distributions are multimodal, thus supporting the hypothesis that a linear relation between magnetic field intensity and plasma density spectral features has to be considered valid only on average being a very crude approximation. The multimodal structure of the observed distribution could also be the counterpart for different physical processes at the origin of the observed fluctuations, or it could be due to the role played by other physical quantities, such as for example the plasma temperature fluctuations, that could modify the relation between electron density and magnetic field second-order scaling exponents. Moreover, a very interesting result stands in the clear bimodal character of the distribution of the magnetic field
exponent in the first LT sector. Indeed, there is not a one-to-one correspondence between the magnetic field
exponent and that of the electron density fluctuations, i.e.,
. A possible explanation could be that, as shown in Gruzinov et al. [
63], we are still in the early stage of a rising bubble so that the spectral features are still not stationary.
We conclude our analysis reporting in
Figure 7 the probability density functions of the location of the plasma bubbles as a function of both the magnetic latitude and LT. These distributions suggest the tendency of plasma bubbles to move away from the equator as the time goes by. This could indicate that the plasma bubbles observed near the midnight sector have not been generated in that LT sector but could have been generated near sunset, have moved at high altitudes, and successively have started to descend in altitude, moving along the magnetic field lines. Indeed, the flux tube intersecting the 460 km Swarm A orbit at ∼10
magnetic latitude has an apex height around 600 km at the equator where the plasma bubbles are generated. This is a crucial point in the interpretation of the LT evolution of the scaling features of magnetic field and density fluctuations inside the plasma bubbles. It makes you understand how, moving away from the sunset, Swarm A crosses older plasma bubbles at locations where the turbulence had time to fully develop.
5. Discussion
Our findings seem to suggest that a non-trivial linear relation exists between the spectral features of the plasma density and magnetic field in the plasma irregularities which depends on the local time and latitude. Actually, the fact that the magnetic field and electron density scaling properties are not simply linearly related is somehow expected. Usually, when we study the plasma bubbles dynamics, we start with a few assumptions.
First, following the derivation of the diamagnetic current in Alken et al. [
64], i.e., considering the standard magnetohydrodynamic equations for a plasma in a steady-state configuration, that is,
, and neglecting gravity, we can write that the plasma pressure force (
) is counter-balanced by the Lorentz force (
), i.e.,
where
is the ambient magnetic field. Combining Equation (
4) with the Maxwell’s equation describing the Ampère’s circuital law in absence of displacement currents
where
is the permeability of free space, we can write
By using the vector calculus identity
we can rewrite Equation (
4) as
The term on the left of Equation (
8) is the sum of plasma pressure force and magnetic pressure one, while on the right, the term represents magnetic tension due to the curvature of the field lines. Taking into account that the plasma structures are much smaller than the geomagnetic field curvature radius, it is possible to assume a linear field line geometry. In this case, we can neglect the magnetic curvature term,
, in the momentum balance so that a change in plasma pressure is, in first approximation, compensated by a change in the magnetic pressure (second assumption). Thus, according to Lühr et al. [
37] and taking into account that
, we can write that across the irregularity walls, the following relation holds:
where
denotes a finite scale difference (i.e., the gradient in Equation (
4) computed at the bubble intensity scales),
is the Boltzmann constant,
is the electron density, and
and
are electron and ion temperatures, respectively.
Third, considering that the strength of the ambient magnetic field is around four orders of magnitude greater than the magnetic field generated by the equatorial plasma bubbles, Equation (
9) can be linearized as:
Fourth, assuming a constant plasma temperature, i.e., assuming a local thermodynamic equilibrium
, according to the latter expression, a high correlation must exist between
and
across plasma bubble walls, and plasma density irregularities can be detected reliably through magnetic field data, being valid the following relation
where
is a constant. Thus, the scaling properties of the electron density and magnetic field data should be equal, leading to a universal character of the statistical properties of fluctuations as
that can be generalized as
Conversely, our analysis find a discrepancy between the theoretical prediction by Lühr et al. [
37] and the experimental results.
As observed by Alken et al. [
64], studying the global flow of the gravity-driven and diamagnetic currents, the diamagnetic effect formula proposed by Lühr et al. [
37] is strictly valid only when the plasma has reached equilibrium and the ambient field has no curvature. The Earth’s magnetic field has a large dipole moment and the assumption to ignore the magnetic field line curvature could be one of the possible reasons responsible for the observed discrepancies. Furthermore, it is however true that the plasma bubbles are unstable systems that evolve rapidly in time. This means that the assumption of a plasma in equilibrium that, in first approximation, may be only valid in the early hours after sunset, is no longer true in the following hours when the turbulence inside plasma bubbles is fully developed [
51]. This point is supported by the discrepancies from a linear relation between magnetic field and plasma density scaling exponents for increasing LT, as shown in
Figure 6, and by the geographic location of plasma bubbles in LT, as shown in
Figure 7. We remark that Equation (
10), being a linear approximation of Equation (
4), is expected to be valid in the first phases of the plasma bubbles rising, i.e., for LT near the dusk–night sector, as highlighted in
Figure 6. There is also an additional issue related to the assumption of a constant plasma (electron and ion) temperature that has to be considered. This might be true in the first instants of plasma bubbles formation, i.e., when plasma bubbles rise just after the sunset and the instability growth is still quasi-linear. Conversely, when these structures start evolving, the emerging of turbulent fluctuations can strongly affect both plasma density and temperature fluctuations which could passively respond to turbulent fluctuations. This means that the assumption of negligible temperature fluctuations could be no longer valid for hours approaching midnight. This is a crucial point that will be investigated in a successive work.
Last but not least, we note that Equation (
10) neglects the role of other terms that might affect the evolution of plasma bubbles, such as for instance the emergence of
convective turbulence and/or Hall MHD compressive turbulence, which might be at the basis of some discrepancies between density and magnetic field intensity fluctuations. In particular, we observe that some of the features of the observed fluctuations of the magnetic field, such as the spectral slope near
are in good agreement with the emergence of 2D
turbulence of electron density, which is indeed isomorphic to the viscous convection of an ordinary fluid driven by temperature gradients [
63]. These extra elements require the investigation of the dynamics of the magnetic field components with a further detail than that discussed here in the case of the magnetic field intensity.
6. Summary and Conclusions
In order to understand whether it is possible to study the dynamical features of plasma density irregularities by using independently either the magnetic field or the electron density measurements, this study focused on determining the relationship between the spectral features of electron density and magnetic field strength inside plasma bubbles. In the past, many important features of the equatorial plasma bubbles have been obtained indirectly from the analysis of their magnetic signatures, using the so-called diamagnetic effect. For example, it has been possible to investigate their dependence on latitude, longitude [
8,
19,
20,
21], and solar and geomagnetic activity [
16,
17] by using magnetic field data recorded by LEO satellites. However, the study of the dynamical features of these plasma density irregularities by using only the magnetic field might not be correct, since it is necessary to assume that the scaling properties of the electron density and magnetic field fluctuations are equal.
In order to address this point, we studied the scaling properties of both the electron density and the magnetic field inside plasma bubbles using measurements on board the Swarm A satellite from 1 April 2014 to 31 January 2016. These were obtained using the so-called generalized
-order structure function analysis evaluated by applying the local detrended structure function analysis (DSFA) method [
50,
51]. The considered method, along with the time resolution of the electron density and magnetic field data, permitted investigating spatial scales ranging from some kilometers to a few hundred kilometers. At this range of spatial scales, our findings suggest that it is not correct to hypothesize that the spectral features of the plasma irregularities can be inferred from their magnetic properties. For example, this type of assumption was done in the past by Lühr et al. [
15], who using 50 Hz magnetic field measurements from CHAMP obtained the spectral properties of the magnetic data and considered them also valid to investigate the spectral features of the electron density irregularities.
Our findings suggest that a more complex relation may exist between the spectral features of the electron density and magnetic field in the plasma irregularities, which depends on the local time and latitude. Indeed, the results of this research support the idea that there are some relevant discrepancies between the spectral features of electron density fluctuations and magnetic field intensity ones. The discrepancies depend on the evolution of plasma bubbles with local time and on its turbulent nature. In the early stage of plasma bubble formation, i.e., when plasma bubbles rise just after the sunset and the instability growth is still quasi-linear, we find a quasi-linear relationship on average between the scaling properties of electron density and those of magnetic field. Conversely, when these structures develop, the evolution of turbulent fluctuations can strongly affect also other physical quantities such as the plasma temperature fluctuations, which could modify the spectral features. This means that the assumption of negligible temperature fluctuations could no longer be valid for hours approaching midnight. In other words, the electron temperature fluctuations could play an important role when turbulence is fully developed.
A better comprehension of the plasma bubbles dynamics and of the turbulence processes that characterize their time evolution may benefit from the use of very high-resolution vector magnetic field and plasma density measurements such as those that should be available from the future NanoMagSat mission [
65]. This future satellite, which should fly at an altitude of 500 km and an inclination angle of
, should allow a quick full local time coverage and for the first time joint measurements of electron density and magnetic signals at very high frequencies. The mission should give the opportunity to investigate plasma variability at very short spatial scales ranging from a few meters to some kilometers.