Journal of Geophysical Research: Space Physics最新文献

We extend the previous three-dimensional (3D) empirical reconstruction (ER) model for a set of ideal magnetohydrodynamics (MHD) constraints into a resistive MHD 3D ER model that includes additional resistive MHD constraints and additional measurements from NASA's Magnetospheric Multiscale (MMS) mission. The same form of a stochastic optimization algorithm is used as in the previous ideal MHD 3D ER model to directly minimize the loss function that includes a few more highly nonlinear terms characterizing the model-measurement differences and the model departures from physical constraints. The resistive MHD 3D ER model is applied to three regions of MMS measurements that correspond to direct sampling of an electron diffusion region (EDR), a region adjacent to the EDR, and one far away from the EDR. The reconstructed plasma and electromagnetic fields are of high quality in all three regions as measured by model-measurement difference indices and physics-based quality indicators. The reconstructed fields in the EDR provide us with a good view of the spatial configuration of the reconnection site. We specifically examine the effect of resistivity on energy exchange in the vicinity of the EDR. It was discovered that in the EDR, the energy exchange shows an exclusive and systematic one-channel process between the plasma thermal energy and electromagnetic energy with the conversion rate highly correlated with the strength of the turbulent electromagnetic fields. In the other two regions away from the EDR, the energy exchange between the electromagnetic energy and the plasma thermal and kinetic energies shows rapidly-varying and random characteristics.

Fundamental processes in plasmas act to convert energies into different forms, for example, electromagnetic, kinetic and thermal. Direct derivation from the Vlasov-Maxwell equation yields sets of equations that describe the temporal evolution of magnetic, kinetic and internal energies in either the monofluid or multifluid frameworks. In this work, we focus on the main terms affecting the changes in kinetic energy. These are pressure-gradient-related terms and electromagnetic terms. The former account for plasma acceleration/deceleration from a pressure gradient, while the latter from an electric field. Although limited spatial and temporal deviations are expected, a statistical balance between these terms is fundamental to ensure the overall conservation of energy and momentum. We use in-situ observations from the Magnetospheric MultiScale (MMS) mission to study the relationship between these terms. We perform a statistical analysis of those parameters in the context of magnetic reconnection by focusing on small-scale Electron Diffusion Regions and large-scale Flux Transfer Events. The analysis reveals a correlation between the two terms in the monofluid force balance, and in the ion force and energy balance. However, the expected relationship cannot be verified from electron measurements. Generally, the pressure-gradient-related terms are smaller than their electromagnetic counterparts. We perform an error analysis to quantify the expected underestimation of gradient values as a function of the spacecraft separation compared to the gradient scale. Our findings highlight that MMS is capable of capturing energy and force balance for the ion fluid, but that care should be taken for energy conversion terms based on electron pressure gradients.

Space weather forecasts are of utmost importance in safeguarding navigation, communication, and electric power system operations, satellites from orbital drag, and the astronauts in the International Space Station from hazardous space radiation during extreme space weather conditions. The finest space weather prediction requires a clear understanding of solar wind-magnetosphere coupling. The in-situ measurements of the solar wind properties give unique information about the Sun and its activity on smaller to longer timescales. The present work investigates the influence of solar activity on the coupling of solar wind and Earth's magnetosphere during 23–24 solar cycles. The geomagnetic storms with Symmetric H-component (SYMH) ≤ −85 nT during the 23–24 solar cycles are considered. We present the results of statistical analysis and relationships between the various solar wind parameters such as the total strength of interplanetary magnetic field (B) and its three-axis components (Bx, By, and Bz), solar wind proton density (Nsw), solar wind speed (Vsw), SYMH indices, the amplitude, duration, and profile of the geomagnetic storms. The integrated electric field and integrated SYMH index during storms show the highest correlation of 0.92, implying that integrated SYMH is a better proxy of the injected solar wind energy in the magnetosphere in the form of the ring current. Moreover, we do see the difference in the solar wind-magnetosphere coupling efficiency during the phases of 23–24 solar cycles which is intriguing.

An approach for creating continental-scale, multi-scale plasma convection maps in the nightside high-latitude ionosphere using the spherical elementary current systems technique has been developed and evaluated. The capability to reconstruct meso-scale flow channels improved dramatically, and the velocity errors were reduced by ∼30% compared to the spherical harmonic fitting method. Uncertainties of velocity vectors estimated by varying the model setup was also low. Convection maps for a substorm event revealed multiple flow channels in the polar cap, dominating the convection in the quiet time and early growth phase. The meso-scale flows extended toward the nightside auroral oval and had continuous flow channels over >20° of latitude, and the flow channels dynamically merged and bifurcated. The substorm onset occurred along one of the flow channels, and the azimuthal extent of the enhanced flows coincided with the initial width of the auroral breakup. During the expansion phase, the meso-scale flows repetitively crossed the oval poleward boundary, and some of them contributed to subauroral polarization streams enhancements. Increased flows extended duskward, along with the westward traveling surge. Then, flows near midnight weakened and evolved to the Harang flow shear. The meso-scale flow channels had significant (∼10%–40% on average) contributions to the total plasma transport. The meso-scale flows were highly variable on ∼10 min time scales and their individual maximum contributions reached upto 73%. These results demonstrate the capability of specifying realistic convection patterns, quantifying the contribution of meso-scale transport, and evaluating the relationship between meso-scale flows and localized auroral forms.

We present a new technique for the upcoming tri-static incoherent scatter radar system EISCAT 3D (E3D) to perform a volumetric reconstruction of the 3D ionospheric electric current density vector field, focusing on the feasibility of the E3D system. The input to our volumetric reconstruction technique are estimates of the 3D current density perpendicular to the main magnetic field, j_⊥, and its covariance, to be obtained from E3D observations based on two main assumptions: (a) Ions fully magnetized above the E region, set to 200 km here. (b) Electrons fully magnetized above the base of our domain, set to 90 km. In this way, j_⊥ estimates are obtained without assumptions about the neutral wind field, allowing it to be subsequently determined. The volumetric reconstruction of the full 3D current density is implemented as vertically coupled horizontal layers represented by Spherical Elementary Current Systems with a built-in current continuity constraint. We demonstrate that our technique is able to retrieve the three dimensional nature of the currents in our idealized setup, taken from a simulation of an active auroral ionosphere using the Geospace Environment Model of Ion-Neutral Interactions (GEMINI). The vertical current is typically less constrained than the horizontal, but we outline strategies for improvement by utilizing additional data sources in the inversion. The ability to reconstruct the neutral wind field perpendicular to the magnetic field in the E region is demonstrated to mostly be within ±50 m/s in a limited region above the radar system in our setup.

We use the full NASA Van Allen Probes mission (2012–2019) to extract the electron plasma density from the Electric and Magnetic Field Instrument Suite and Integrated Science (EMFISIS) and Electric Field and Waves (EFW) instruments and discuss the evolution of the plasmasphere. We generate new statistics including mean and standard deviations of the plasma density with respect to L-shell, magnetic local time (MLT), and various geomagnetic indices. These statistics are generated to be applied in radiation belt physics and space weather codes (with fits provided). The mean plasmasphere is circular around Earth with respect to MLT for Kp < 1. The mean 100 cm⁻³ level line is above L = 5 and mean 10 cm⁻³ level expands above the Van Allen Probes apogee for Kp < 1. The outer electron belt lies within the plasmasphere for 60% of all times. As activity increases (Kp > 2), a gradual MLT asymmetry forms with higher mean density in the afternoon sector due to plumes expanding outward. Conversely, the mean density decreases on the dawn and night sectors. The mean density is between ∼500 and ∼50 cm⁻³ between L ∼ 4 and L ∼ 6 during quiet and moderately active times (Kp < 3), representing ∼80% of all times. Statistics in regions of high density below L = 2 are underdefined for intense activity. The highest standard deviation of density represents a factor 2.5 to 3 times the mean above L = 5 and for active times. We find the percent difference between the EFW and EMFISIS densities is bounded by ±20% for quiet and moderate activity (Kp < 5) and goes up to ±100% for extreme activity.

We survey crossings of Jupiter's dawn magnetopause during the Juno prime mission to identify and characterize Jupiter's magnetopause boundary layer. Using plasma and magnetic field observations from Jovian Auroral Distributions Experiment and Juno Magnetic Field investigation, we identify 53 boundary layer events from the 62 magnetopause crossings studied here. We find that the boundary layer generally exhibits mixed properties of magnetosheath and magnetosphere electron distributions, including lower characteristic electron energies and denser ion populations than in the magnetosphere, but higher characteristic electron energies and less dense ion populations than in the magnetosheath. Boundary layer proton speeds are on average slower than both the magnetosheath and magnetosphere. Other proton parameters in the boundary layer have intermediate values between the magnetosheath and magnetosphere. Through ion composition analysis in regions adjacent to the magnetopause, we find evidence of solar wind and magnetospheric plasma in the boundary layer that suggests plasma is transported across the magnetopause in both directions. This mass and energy transport may be the result of solar wind interactions such as magnetic reconnection and Kelvin-Helmholtz instabilities. However, many boundary layer events do not exhibit local signatures of these solar wind interactions and plasma may be transported by a non-local process or diffusively transported.

Observations using Arecibo Observatory's highly sensitive Incoherent Scattering Radar (AO-ISR) show ionospheric descending layers from as high as ∼400 km, much higher than earlier studies, with continuity down to 90 km. The AO-ISR was operated to observe the ion-line and plasma-line with coded-long-pulse for high temporal and spatial resolution of 35/10 s and 300 m, respectively, during 01–06 February 2019. We found multiple layering structures descending from 400 to 90 km in all these six days. These layers are traditionally called intermediate descending layers (IDLs) (>130 km and below F-peak), upper semi-diurnal daytime and nighttime layers (110–130 km), and lower diurnal layers (<110 km). We have denoted the new daytime descending layers above the hmF2 as top-side descending layers (TDLs). All these layers are collectively named ionospheric descending layers (IonDLs) since all of them are connected with some discontinuity at the F1-peak (i.e., 170 km), except for the daytime lower-diurnal layer. The most pronounced IonDLs occur in the twilight times. IonDLs mainly occur in shear zones of the vertical ion drifts and are favored by downward ion drifts, and their descent speeds increase with increasing altitude. The estimated phase velocities of the waves in the F-region are comparable with the descending speed of the IonDLs. Furthermore, IonDLs/IDLs occur with and without spread-F events but intensified spread-F events raise their beginning altitude. The TDLs and IDLs are driven by gravity waves with periods of 1.5–4 hr.

Data products relating to auroral arc systems are often sparse and distributed while ionospheric simulations generally require spatially continuous maps as boundary conditions at the topside ionosphere. Fortunately, all-sky auroral imagery can provide information to fill in the gaps. This paper describes three methods for creating electrostatic plasma convection maps from multi-spectral imagery combined with plasma flow data tracks from heterogeneous sources. These methods are tailored to discrete arc structures with coherent morphologies. The first method, “reconstruction,” builds the electric potential map (from which the flow field is derived) out of numerous arc-like ridges that are then optimized against the plasma flow data. This method is designed for data from localized swarms of spacecraft distributed in both latitude and longitude. The second method, “replication,” uses a 1D across-arc flow data track and replicates these data along a determined primary and secondary arc boundary while simultaneously scaling and rotating to keep the flow direction parallel to the arc and the flow shear localized at the arc boundaries. The third, “weighted replication,” performs a replication on two data tracks and calculates a weighted average between them, where the weighting is based on data track proximity. This paper shows the use of these boundary conditions in driving and assessing 3D auroral ionospheric, multi-fluid simulations.

Pc5 ULF waves play an important role in transporting energy and particles in the coupled magnetospheric and ionospheric system. They are known to be initiated by dynamic pressure fluctuations upstream of the magnetopause, including those induced by hot flow anomalies (HFAs). However, the role of HFAs in generating magnetospheric and ground magnetic Pc5 ULF oscillations has not been investigated statistically yet. Thus, in this paper, we investigate the contribution of HFAs to ground magnetic Pc5 ULF oscillations and analyze how the characteristics of HFAs influence these oscillations, based on the coordinated observations between the THEMIS probes and the ground magnetometers at high latitudes during the years 2008, 2009 and 2019. We find that HFAs can serve as a notable source of ground magnetic Pc5 ULF oscillations, with about 18.9% of Interplanetary Magnetic Field (IMF) discontinuity-induced HFAs associated with discernible enhancements in Pc5 ULF wave power, whereas spontaneous HFAs play a comparatively minor role in generating these oscillations. Furthermore, we observe that the cores of HFAs are likely to contribute more significantly to modulating the induced ground magnetic Pc5 ULF oscillations than their compressed boundaries. More dynamic pressure reductions within HFA cores correspond to stronger ground magnetic Pc5 ULF oscillations. Additionally, HFAs can propagate with the IMF discontinuity along the bow shock, continuously generating ground magnetic Pc5 ULF oscillations during their propagation. This research sheds light on the mechanisms underlying Pc5 ULF wave generation and underscores the role of HFAs in driving magnetospheric-ionospheric interactions.