Solar Wind Convection Electric Field Research Articles

Effects of an Intrinsic Magnetic Field on Ion Escape From Ancient Mars Based on MAESTRO Multifluid MHD Simulations

AbstractIon escape has played a key role in atmospheric loss on ancient Mars due to intense solar activity. Under the existence of a strong global intrinsic magnetic field, as expected on ancient Mars, differential flows between different ion species can be important for ion escape. To assess effects of differential flows, we here developed a new global multifluid magnetohydrodynamics model with a cubed sphere grid named Multifluid Atmospheric Escape Simulations Toward Real elucidatiOn (MAESTRO). A test simulation under present‐day Mars conditions showed solar wind‐Mars interactions, for example, plasma boundaries, ionospheric profiles, and tailward outflows, consistent with observations and simulation studies. We then conducted multifluid and multispecies simulations with six different dipole field strengths under ancient Mars conditions. The multifluid cases show asymmetric outflow distributions with respect to the solar wind convection electric field, as pointed out by previous studies on present‐day Mars. Compared with multispecies cases, the multifluid representation increases the escape rates of O2+ and CO2+ by more than two orders of magnitude in the strong dipole field cases where the cusp outflow is dominant. The O+ escape rate is slightly lower in the multifluid cases with no or weak dipole field due to suppression of ion pickup in the −E hemisphere, while it is reduced by one order of magnitude in the strongest dipole field case. The existence of a dipole field can reduce the total escape rate by a factor of six. The impact on atmospheric evolution is diminished but still significant.

The Influence of Solar Irradiation and Solar Wind Conditions on Heavy Ion Escape from Mars

AbstractWe apply a recently proposed method to estimate heavy ion escape from Mars. The method combines in situ observations with a hybrid plasma model, which treats ions as particles and electrons as a fluid. With this method, we investigate how solar upstream conditions, including solar extreme ultraviolet (EUV) radiation, solar wind dynamic pressure, and interplanetary magnetic field (IMF) strength and cone angle, affect the heavy ion loss. The results indicate that the heavy ion escape rate is greater in high EUV conditions. The escape rate increases with increasing solar wind dynamic pressure, and decreases as the IMF strength increases. The ion escape rate is highest when the solar wind is parallel to the IMF and lowest when they are perpendicular. The plume escape rate decreases when the solar wind convective electric field increases.

Forces, electric fields and currents at the subsolar martian MPB: MAVEN observations and multifluid MHD simulation

The Martian MPB (Magnetic Pileup Boundary) is a key boundary in the Mars/Solar Wind interaction as it is here that part of the momentum and energy from the solar wind plasma are transferred to the planetary plasma. Since this interaction is for the most part collisionless, the transfer is mediated by electric and magnetic fields. The acceleration processes and the interaction of particles with electromagnetic fields operate at spatial scales determined by the ambient particle populations. In particular, in regions with sizes of the order of the ion inertial length (ion scales), the Hall electric field is expected to be dominant.In the present work we combine data from the MAVEN spacecraft along one orbit around Mars and multifluid MHD simulation results to study the role of electric fields, currents and forces at the MPB at ion scales. In particular, we find that the current densities deduced from MAVEN data (J ∼ 238 nA/m2) of the same order as the values obtained in the simulation (J ∼ 56-156 nA/m2) and that the Hall electric force points sunward in both cases. In addition, we find that in the subsolar MPB current layer the Hall electric field (∼3.2 mV/m) dominates over the solar wind convective electric field (∼0.4 mV/m) and electron pressure gradient (∼0.8 mV/m). These values are consistent with previous results suggesting that the MPB thickness is of the order of the solar wind proton inertial length and support the idea that non ideal terms in Ohm’s law must be considered when analysing the dynamics of particles around plasma boundaries with ion scale thicknesses.

Correlating the interplanetary factors to distinguish extreme and major geomagnetic storms

We investigate the correlation between Disturbance Storm Time (Dst) characteristics and solar wind conditions for the main phase of geomagnetic storms, seeking possible factors that distinguish extreme storms (minimum Dst <−250 nT) and major storms (minimum Dst <−100 nT). In our analysis of 170 storms, there is a marked correlation between the average rate of change of Dst during a storm’s main phase (ΔDst/Δt) and the storm’s minimum Dst, indicating a faster ΔDst/Δt as storm intensity increases. Extreme events add a new regime to ΔDst/Δt, the hourly time derivative of Dst (dDst/dt), and sustained periods of large amplitudes for southward interplanetary magnetic field Bz and solar wind convection electric field Ey. We find that Ey is a less efficient driver of dDst/dt for extreme storms compared to major storms, even after incorporating the effects of solar wind pressure and ring current decay. When minimum Dst is correlated with minimum Bz, we observe a similar divergence, with extreme storms tending to have more negative Dst than the trend predicted on the basis of major storms. Our results enable further improvements in existing models for storm predictions, including extreme events, based on interplanetary measurements.

Singing Comet Waves in a Solar Wind Convective Electric Field Frame

AbstractThe cometary plasma environment at low gas production rates is dominated by highly compressional, large‐amplitude magnetic field waves in the 10–100 mHz range. They are thought to be caused by an ion‐Weibel instability due to a cross‐field current, which is caused by the cometary ions that are accelerated along the solar wind convective electric field. We devise a new method to determine the location of the wave source, the wave power, frequency, and bandwidth. It is found that the wave occurs everywhere in the coma, regardless of electric field direction. There is no correlation between the wave frequency and the measured plasma density. This is not in agreement with previous studies. A dependence of the frequency on the position of the spacecraft in a comet‐fixed frame is in agreement with the prediction from the ion‐Weibel instability. We infer a wave generation region much larger than the cometocentric distances covered by Rosetta.

Asymmetric Lunar Magnetic Perturbations Produced by Reflected Solar Wind Particles

Abstract Magnetic perturbations characterize the solar wind interaction of the Moon. The solar wind plasma absorption on the dayside surface produces large-scale field perturbations behind, i.e., the field enhancement in the central wake and reduction on the wake boundary. The solar wind repellence over local lunar magnetic anomalies (LMAs) leads to small-scale magnetic compressions ahead. In this study, the magnetic perturbations around the Moon are examined by using the observations from a near-Moon satellite mission, the Lunar Prospector, and they exhibit a clear left–right asymmetry in a coordinate system related to the solar wind convection electric field ( ). The magnetic field is observed to enhance before the left terminator that points to, while on the opposite side, it is not. The test particle simulations show that can divert the particles reflected over the LMAs to the left and then the solar wind pickup of these particles leads to the field enhancement observed before the left terminator. Behind the lunar terminator, the wake field reduction is also asymmetric. On the left, the field reduction is more remarkable and located closer to the central wake. The denser plasma, consisting of the background as well as the reflected solar wind particles, may produce a stronger diamagnetic current and thus more significant field reduction there. The asymmetric plasma and magnetic perturbations associated with the reflected particles may be a common and nonnegligible element during the solar wind interaction of a small-scale magnetic field, such as that of an asteroid or a comet.

Statistical analysis of interplanetary coronal mass ejections and their geoeffectiveness during the solar cycles 23 and 24

Interplanetary coronal mass ejections (ICMEs) are believed to be the most common and important drivers of the strongest geomagnetic storms. In this work, the geoeffective characteristics of the ICMEs occurred during the last solar cycles 23 and 24 (years 1996–2017) have been studied in detail. The maximum velocity $V$max, either of the ICME’s Sheath region or of the ICME itself, the mean velocity of the ICME, the minimum value of the southward component of the Interplanetary magnetic field Bs and the $y$-component of the solar wind convective electric field $E= -V \times B$ observed at L1 point during the pass of the ICME, were used. It was found that, in accordance to past similar studies, the most dominant characteristic of ICMEs in the generation of geomagnetic storms is the Bs component along with the Ey parameter, while the maximum velocity seems to be of less importance. Nevertheless the maximum speed is an good forecasting factor due to the fact that it is much easier to estimate the velocity of an ICME-structure many hours before it arrives at Earth compared with the observations of Bs and Ey that can only be done, for the time being, at the L1 point. That means we can use the velocity of ICME-structure to forecast the possible generation and magnitude of the geomagnetic storms. From a comparison of the ICME-generated geomagnetic storms with the total number of geomagnetic storms generated during the last two solar cycles, it seems that approximately half of the ICMEs (49% for Dst index and 53% for Kp index) produced geomagnetic storms during the solar cycles 23 and 24. Moreover the velocities of ICMEs are more in accordance with the rising and maximum phases of solar cycles 23 and 24 than the geomagnetic activity (storms) are, as well as during the first stages of the declining phases of these cycles, especially during solar cycle 23.

Diamagnetic Plasma Clouds in the Near Lunar Wake Observed by ARTEMIS

Abstract Because of the plasma pressure gradient across the lunar wake boundary, the interplanetary magnetic field is enhanced in the lunar wake. In previous works, the solar wind ions that enter into the near lunar wake are found to have an unimportant influence on the magnetic field within the lunar wake. In this study, two cases of a 22%–70% reduction of the magnetic field in the near lunar wake are first observed by the Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon’s Interaction with the Sun mission (ARTEMIS). The magnetic field depressions are caused by the refilling of dense plasma clouds (with densities of 0.20–0.47 cm−3) into the near lunar wake. The ions of the plasma clouds that originate from the reflected solar wind ions in the lunar dayside are accelerated into the near lunar wake by the solar wind convection electric field. The source regions of the reflected ions can be traced back to the lunar magnetic anomalies over the sunlit surface. Pressure (magnetic pressure + thermal pressure) balances are roughly maintained for both cases at the boundary between the plasma cloud and the ambient plasma. Our results imply that the interaction between the solar wind and lunar magnetic anomalies drastically disturbs the near-Moon electromagnetic environment.

The Convective Electric Field Influence on the Cold Plasma and Diamagnetic Cavity of Comet 67P

Abstract We studied the distribution of cold electrons (&lt;1 eV) around comet 67P/Churyumov–Gerasimenko with respect to the solar wind convective electric field direction. The cold plasma was measured by the Langmuir Probe instrument and the direction of the convective electric field conv = − × was determined from magnetic field ( ) measurements inside the coma combined with an assumption of a purely radial solar wind velocity . We found that the cold plasma is twice as likely to be observed when the convective electric field at Rosetta’s position is directed toward the nucleus (in the − conv hemisphere) compared to when it is away from the nucleus (in the + conv hemisphere). Similarly, the diamagnetic cavity, in which previous studies have shown that cold plasma is always present, was also found to be observed twice as often when in the − conv hemisphere, linking its existence circumstantially to the presence of cold electrons. The results are consistent with hybrid and Hall magnetohydrodynamic simulations as well as measurements of the ion distribution around the diamagnetic cavity.

Statistical Study of Heavy Ion Outflows From Mars Observed in the Martian‐Induced Magnetotail by MAVEN

AbstractIt is important to include the effects of cold ions when we consider heavy ion outflows from Mars. We here report on statistical properties of heavy ion outflows (including cold ion outflows) observed by Mars Atmosphere and Volatile EvolutioN in the optical wake region. Using data from July 2015 to December 2017, we statistically investigate the effects of solar wind convection electric fields and crustal magnetic fields on the heavy ion outflows. Results show that the average density ratio of O+:O2+:CO2+ is ~29:68:04. In the southern hemisphere where the strong crustal magnetic fields are located, the heavy ion outflow flux becomes smaller and the relative contribution of molecular ions to heavy ion outflows is larger than the northern hemisphere. The solar wind convection electric field strongly affects the heavy ion outflows. Heavy ion density is larger in the −E (electric field) hemisphere than in the +E hemisphere, while the dependence of velocity is opposite. Acceleration by the solar wind convection electric field in the +E hemisphere is expected to cause these dependencies. The heavy ion flux is larger in the −E (total O: 8.6 × 106 cm–2/s) than in the +E hemisphere (2.9 × 106 cm–2/s) due to the large density in the −E hemisphere. Velocity ratios of O+ to O2+ suggest that heavy ion outflows with large velocities tend to have the same energy to each other, while the O+ to O2+ ions are more likely to have the same velocity in outflows with small velocities.

Oxygen Ion Energization at Mars: Comparison of MAVEN and Mars Express Observations to Global Hybrid Simulation

AbstractWe study oxygen ion energization in the Mars‐solar wind interaction by comparing particle and magnetic field observations on the Mars Atmosphere and Volatile EvolutioN (MAVEN) and Mars Express missions to a global hybrid simulation. We find that large‐scale structures of the Martian‐induced magnetosphere and plasma environment as well as the Mars heavy ion plume as seen by multispacecraft observations are reproduced by the model. Using the simulation, we estimate the dynamics of escaping oxygen ions by analyzing their distance and time of flight as a function of the gained kinetic energy along spacecraft trajectories. In the upstream region the heavy ion energization resembles single‐particle solar wind ion pickup acceleration as expected, while within the induced magnetosphere the energization displays other features including the heavy ion plume from the ionosphere. Oxygen ions take up to 80 s and travel the distance of 20,000 km after their emission from the ionosphere to the induced magnetosphere or photoionization from the neutral exosphere before they have reached energies of 10 keV in the plume along the analyzed spacecraft orbits. Lower oxygen ion energies of 100 eV are reached faster in 10–20 s over the distance of 100–200 km in the plume. Our finding suggests that oxygen ions are typically observed within the first half of their gyrophase if the spacecraft periapsis is on the hemisphere where the solar wind convection electric field points away from Mars.

Influence of the Interplanetary Convective Electric Field on the Distribution of Heavy Pickup Ions Around Mars

AbstractThis study obtains a statistical representation of 2–15 keV heavy ions outside of the Martian‐induced magnetosphere and depicts their organization by the solar wind convective electric field (ESW). The overlap in the lifetime of Mars Global Surveyor (MGS) and Mars Express (MEX) provides a period of nearly three years during which magnetometer data from MGS can be used to estimate the direction of ESW in order to better interpret MEX ion data. In this paper we use MGS estimates of ESW to express MEX ion measurements in Mars‐Sun‐Electric field (MSE) coordinates. A new methodological technique used in this study is the limitation of the analysis to a particular instrument mode for which the overlap between proton contamination and plume observations is rare. This allows for confident energetic heavy ion identification outside the induced magnetosphere boundary. On the dayside, we observe high count rates of 2–15 keV heavy ions more frequently in the +ESW hemisphere (+ZMSE) than in the −ESW hemisphere, but on the nightside the reverse asymmetry was found. The results are consistent with planetary origin ions being picked up by the solar wind convective electric field. Though a field of view hole hinders quantification of plume fluxes and velocity space, this new energetic heavy ion identification technique means that Mars Express should prove useful in expanding the time period available to assess general plume loss variation with drivers.

MAVEN observations of tail current sheet flapping at Mars

AbstractThe Martian magnetotail is a complex regime through which atmospheric particles are lost to space. Our current understanding of Mars' tail continues to develop with the comprehensive particle and field data collected by Mars Atmosphere and Volatile EvolutioN (MAVEN). In this work, we identify periods when MAVEN encounters multiple current sheet crossings through a single tail traversal in order to understand tail dynamics. We apply an analysis technique that has been developed and validated by using multipoint measurements in order to separate the spatial and temporal properties associated with current sheet flapping. Events are classified into periods of steady flapping, due to a global motion of the current sheet, and kink‐like flapping, resulting from localized wave propagation along the tail current sheet. Out of 106 periods during which multiple current sheet crossings were observed, 20 were due to steady flapping and 10 from kink‐like flapping. A majority of the kink‐like events resulted from waves propagating in the opposite direction of the solar wind convection electric field, regardless of their location in the tail, unlike at Earth and Venus. This finding suggests that possible magnetosphere energy sources, whereby plasma is accelerated and removed from the Martian environment, are not located in the central magnetotail; rather, these waves may be driven by a source located at the tail flank based on the direction of the solar wind electric field. Therefore, by identifying potential sources of impulsive energy release in the tail, we may better understand mechanisms that drive atmospheric loss at Mars.

Mass loading at 67P/Churyumov‐Gerasimenko: A case study

AbstractWe study the dynamics of the interaction between the solar wind ions and a partially ionized atmosphere around a comet, at a distance of 2.88 AU from the Sun during a period of low nucleus activity. Comparing particle data and magnetic field data for a case study, we highlight the prime role of the solar wind electric field in the cometary ion dynamics. Cometary ion and solar wind proton flow directions evolve in a correlated manner, as expected from the theory of mass loading. We find that the main component of the accelerated cometary ion flow direction is along the antisunward direction and not along the convective electric field direction. This is interpreted as the effect of an antisunward polarization electric field adding up to the solar wind convective electric field.

Emission of hydrogen energetic neutral atoms from the Martian subsolar magnetosheath

AbstractWe have simulated the hydrogen energetic neutral atom (ENA) emissions from the subsolar magnetosheath of Mars using a hybrid model of the proton plasma charge exchanging with the Martian exosphere to study statistical features revealed from the observations of the Neutral Particle Detectors on Mars Express. The simulations reproduce well the observed enhancement of the hydrogen ENA emissions from the dayside magnetosheath in directions perpendicular to the Sun‐Mars line. Our results show that the neutralized protons from the shocked solar wind are the dominant ENA population rather than those originating from the pickup planetary ions. The simulation also suggests that the observed stronger ENA emissions in the direction opposite to the solar wind convective electric field result from a stronger proton flux in the same direction at the lower magnetosheath; i.e., the proton fluxes in the magnetosheath are not cylindrically symmetric. We also confirm the observed increasing of the ENA fluxes with the solar wind dynamical pressure in the simulations. This feature is associated with a low altitude of the induced magnetic boundary when the dynamic pressure is high and the magnetosheath protons can reach to a denser exosphere, and thus, the charge exchange rate becomes higher. Overall, the analysis suggests that kinetic effects play an important and pronounced role in the morphology of the hydrogen ENA distribution and the plasma environment at Mars, in general.

Strong plume fluxes at Mars observed by MAVEN: An important planetary ion escape channel

AbstractWe present observations by the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission of a substantial plume‐like distribution of escaping ions from the Martian atmosphere, organized by the upstream solar wind convection electric field. From a case study of MAVEN particle‐and‐field data during one spacecraft orbit, we identified three escaping planetary ion populations: plume fluxes mainly along the upstream electric field over the north pole region of the Mars‐Sun‐Electric field (MSE) coordinate system, antisunward ion fluxes in the tail region, and much weaker upstream pickup ion fluxes. A statistical study of O+ fluxes using 3 month MAVEN data shows that the plume is a constant structure with strong fluxes widely distributed in the MSE northern hemisphere, which constitutes an important planetary ion escape channel. The escape rate through the plume is estimated to be ~30% of the tailward escape and ~23% of the total escape for &gt; 25 eV O+ ions.

The spatial distribution of planetary ion fluxes near Mars observed by MAVEN

AbstractWe present the results of an initial effort to statistically map the fluxes of planetary ions on a closed surface around Mars. Choosing a spherical shell ~1000 km above the planet, we map both outgoing and incoming ion fluxes (with energies &gt;25 eV) over a 4 month period. The results show net escape of planetary ions behind Mars and strong fluxes of escaping ions from the northern hemisphere with respect to the solar wind convection electric field. Planetary ions also travel toward the planet, and return fluxes are particularly strong in the southern electric field hemisphere. We obtain a lower bound estimate for planetary ion escape of ~3 × 1024 s−1, accounting for the ~10% of ions that return toward the planet and assuming that the ~70% of the surface covered so far is representative of the regions not yet visited by Mars Atmosphere and Volatile EvolutioN (MAVEN).

Dynamics of planetary ions in the induced magnetospheres of Venus and Mars

We compare dynamics of planetary ions in the induced magnetospheres of Venus and Mars in a global hybrid simulation to study factors controlling the ion escape at unmagnetized planets. In the simulation we find that the finite Larmor radius (FLR) effects of escaping heavy ions are stronger at Mars than Venus under nominal solar wind conditions. But, varying upstream conditions, especially the IMF, affects the strength of the FLR effects. We classify three basic types of planetary ion dynamics in an induced magnetosphere. First, light ions such as hydrogen follow the E×B drift, and escape in the wake in the hemisphere where the solar wind convection electric field is pointing towards the planet. Second, heavy ions like oxygen undergo FLR effects, and escape mainly outside of the wake in the hemisphere where the solar wind convection electric field is pointing away from the planet. Third, ion species between light and heavy ions can have both the E×B and FLR type dynamics in the same time.

Lack of relationship between geoeffectiveness and orientations of magnetic clouds with bipolar Bz and unipolar southward Bz

In this study, 38 magnetic clouds (MCs) that caused significant geomagnetic storms (the minimum SYM-H, SHmin, ≤−50nT) are examined, in which 17 MCs were unipolar Bz in south (S-type) and 21 MCs were bipolar Bz (north-to-south, NS-type, or south-to-north, SN-type). For S-type MC, inclination angle of the axis of the MC, |θ|, is ≥45°, while |θ|<45° for bipolar MC. This paper aims to address a question: is the intensity of a MC-driven storm correlated with the orientations of bipolar and S-type MCs? Our results demonstrate that there is no direct and significant relationship between geoeffectiveness and orientations of bipolar and S-type MCs. In other words, there is no MC preference (bipolar or S-type MC) to regulate the SHmin of the storm. On the whole, the SHmin is found to strongly correlate with southward field Bz (cc=0.96) and with the y component of the solar wind convective electric field (cc=−0.91) but to weakly correlate with solar wind speed (cc=−0.65). This result is consistent with previous studies by Wu and Lepping (2002), J. Geophys. Res. 107 (A10), 1314. doi:10.1029/2001JA000161. By separating MC-driven storms by size into moderate (−100nT<SHmin≤−50nT), intense (−200nT<SHmin≤−100nT), and super-intense (SHmin≤−200nT) categories, we find that (1) there was a high-speed stream at the trailing of NS-type MCs for moderate and intense storms, and (2) for super-intense storms driven by bipolar or S-type MCs, both solar wind speed and dynamic pressure had a significant jump at the same time before the onset of the storm, as compared to the other two categories. These jumps are the signatures of a shock. We conclude that the SHmin of a MC-driven storm is regulated by the cloud’s southward field Bz and the y component of the electric field, regardless of whether the MC is bipolar or S-type. Also, the ambient solar wind structure (e.g., shock) ahead of MC may play a role in regulating the storm’s intensity.

Hemispheric asymmetries of the Venus plasma environment

We study the Venus‐solar wind interaction and the hemispheric asymmetries of the Venus plasma environment in the global HYB‐Venus hybrid simulation. We concentrate especially on the role of the flow‐aligned interplanetary magnetic field (IMF) component (i.e., the Parker spiral angle or the IMF cone angle) and analyze the dawn‐dusk and Esw asymmetries between four magnetic quadrants around Venus. Using the simulation model, we study two upstream condition cases in detail: the perpendicular IMF to the solar wind flow case and the nominal Parker spiral case (dominant flow‐aligned IMF component). Several differences and similarities were found in these two simulation runs. Common features of the Venus plasma environment between the two cases include asymmetric magnetic barrier and tail lobes and asymmetric planetary ion escape in the direction of the solar wind convection electric field. Further, protons of planetary origin and of solar wind origin were found to follow similar velocity patterns in the Venus plasma wake in both cases. The differences when the IMF flow‐aligned component is dominating compared to the perpendicular IMF case, the so‐called (magnetic) dawn‐dusk asymmetries, include the parallel bow shock and the foreshock region, the asymmetric magnetic barrier, the asymmetric tail current system, and the asymmetric central tail current sheet. Further, the escaping planetary H+ and O+ ion fluxes are concentrated more on the hemisphere of the parallel bow shock. When interpreting in situ plasma and magnetic observations from Venus, the features of at least these two basic IMF configurations should be considered.