Umeå universitets logga

Aims. We study the instantaneous asymmetry of the Martian bow shock during a change in the direction of the interplanetary magnetic field (IMF) and for steady-state conditions. Specifically, we study the asymmetry with regard to the convective electric field and to the crustal fields of Mars.

Methods. Two methods were used: First, a single-spacecraft method in which a switch in hemisphere in the Mars solar-electric (MSE) coordinate system was studied during a change in the direction of the interplanetary magnetic field. Second, we used a dual-spacecraft method wherein near simultaneous bow shock crossings on opposite hemispheres were studied. The dual bow shock crossings were then compared to a bow shock model, and the difference in the distance to the model was used as a measure of asymmetry.

Results. With the single-spacecraft method, an asymmetry with respect to the solar wind convective electric field, Esw, was found, wherein the bow shock was farther from the planet in the ZMSE <0 hemisphere, that is, the - E hemisphere. With the dual-spacecraft method, the mean of the magnitude of the asymmetries in the individual case was 0.13 RM. However, the standard deviation was as high as the mean, and no significant asymmetry could be attributed either to the solar wind convective electric field or to the Martian crustal fields. A strong asymmetry without a clear correlation to these factors was found nonetheless. Possible causes of the measured asymmetry are discussed.

Conclusions. The magnitude of the asymmetries in individual observations is larger than the average asymmetries. This indicates that the shape of the Martian bow shock is dynamic and influenced by fluctuations or wave phenomena.

Plasma flows with enhanced dynamic pressure, known as magnetosheath jets, are often found downstream of collisionless shocks. As they propagate through the magnetosheath, they interact with the surrounding plasma, shaping its properties, and potentially becoming geoeffective upon reaching the magnetopause. In recent years (since 2016), new research has produced vital results that have significantly enhanced our understanding on many aspects of jets. In this review, we summarise and discuss these findings. Spacecraft and ground-based observations, as well as global and local simulations, have contributed greatly to our understanding of the causes and effects of magnetosheath jets. First, we discuss recent findings on jet occurrence and formation, including in other planetary environments. New insights into jet properties and evolution are then examined using observations and simulations. Finally, we review the impact of jets upon interaction with the magnetopause and subsequent consequences for the magnetosphere-ionosphere system. We conclude with an outlook and assessment on future challenges. This includes an overview on future space missions that may prove crucial in tackling the outstanding open questions on jets in the terrestrial magnetosheath as well as other planetary and shock environments.

In the Earth's magnetosphere wave-particle interaction is a major ion energization process, playing an important role for the atmospheric escape. A common type of ion heating is associated with low-frequency broadband electric wave fields. For such waves the energy is not concentrated to a certain narrow frequency range and exhibits no peaks or dips in a power spectrum. If there are enough fluctuations close to the ion gyrofrequency the electric field may still come in resonance with gyrating ions and heat them perpendicular to the background magnetic field. We perform a proof-of-concept study to investigate if this heating mechanism may contibute significantly to the energization of planetary ions also in the induced magnetosphere of Venus. We assume Alfvénic fluctuations and estimate the electric field spectral density based on magnetic field observations. We find typical estimated electric spectral densities of a few (Formula presented.) /Hz close to Venus. This corresponds to a heating rate of a few eV/s. We consider an available interaction time of (Formula presented.) 300 s and conclude that this mechanism could increase the energy of an oxygen ion by about a keV. Observed thermal energies are in the range 100–1,000 eV and thus, resonant wave heating may also be important at Venus.

At a low activity comet the plasma is distributed in an asymmetric way. The hybrid simulation code Amitis is used to look at the spatial evolution of ion velocity distribution functions (VDFs), from the upstream solar wind (SW) to within the comet magnetosphere where the SW is heavily mass-loaded by the cometary plasma. We find that the spatial structures of the ions and fields form a highly asymmetric induced magnetosphere. The VDFs of SW and cometary ions vary drastically for different locations in the comet magnetosphere. The shape of the VDFs differ for different species. The SW protons show high anisotropies that occasionally resemble partial rings, in particular at small cometocentric distances. A second, decoupled, proton population is also found. Solar wind alpha particles show similar anisotropies, although less pronounced and at different spatial scales. The VDFs of cometary ions are mostly determined by the structure of the electric field. We perform supplementary dynamic particle backtracing to understand the flow patterns of SW ions that lead to these anisotropic distributions. This tracing is needed to understand the origin of cometary ions in a given part of the comet magnetosphere. The particle tracing also aids in interpreting observed VDFs and relating them to spatial features in the electric and magnetic fields of the comet environment.

Context. The direction of the interplanetary magnetic field determines the nature of the interaction between a Solar System object and the solar wind. For comets, it affects the formation of both a bow shock and other plasma boundaries, as well as mass-loading. Around the nucleus of a comet, there is a diamagnetic cavity, where the magnetic field is negligible. Observations by the Rosetta spacecraft have shown that, most of the time, the diamagnetic cavity is located within a solar-wind ion cavity, which is devoid of solar wind ions. However, solar wind ions have been observed inside the diamagnetic cavity on several occasions. Understanding what determines whether or not the solar wind can reach the diamagnetic cavity also advances our understanding of cometsolar wind interaction in general.

Aims. We aim to determine the influence of an interplanetary magnetic field directed radially out from the Sun that is, parallel to the solar wind velocity on the cometsolar wind interaction. In particular, we explore the possibility of solar wind protons entering the diamagnetic cavity under radial field conditions.

Methods. We performed global hybrid simulations of comet 67P/Churyumov-Gerasimenko using the simulation code Amitis for two different interplanetary magnetic field configurations and compared the results to observations made by the Rosetta spacecraft.

Results. We find that, when the magnetic field is parallel to the solar wind velocity, no bow shock forms and the solar wind ions are able to enter the diamagnetic cavity. A solar wind ion wake still forms further downstream in this case.

Conclusions. The solar wind can enter the diamagnetic cavity if the interplanetary magnetic field is directed radially from the Sun, and this is in agreement with observations made by instruments on board the Rosetta spacecraft.

Diamagnetic cavities at comets were predicted already in the 1960s [1], and then observed at comet lP/Halley by the ESA/Giotto spacecraft in 1986 [2]. Recently, the ESA/Rosetta spacecraft spent two years orbiting comet 67P/Churyumov-Gerasimenko and encountered the diamagnetic cavity of comet 67P more than 700 times [3, 4]. Most encounters lasted a few minutes, with the duration varying from a few seconds up to more than 30 minutes. As the spacecraft moved very slowly (~lms^-1), it can be considered stationary with respect to the plasma. Therefore, the quick succession of detections indicates that the boundary of the diamagnetic cavity moved over the spacecraft. Figure 1 (left) shows three diamagnetic cavity signatures observed with the plasma instruments on Rosetta on 16 September 2015 when the comet was close to perihelion. Rosetta was in the diamagnetic cavity during the periods of nearly zero magnetic field (marked by the coloured regions). Outside the cavity, the plasma was often characterised by a series of asymmetric, steepened waves which are visible in the magnetic field, as well as in the plasma density [5]. Since all observations to date have been made using a single spacecraft, the shape of the diamagnetic cavity boundary cannot be well constrained by measurements. However, it has been suggested, based on wave observations, that bulges on the cavity boundary move past the spacecraft, causing the latter to quickly move in and out of the cavity [6].

Context. In past decades, several spacecraft have visited comets to investigate their plasma environments. In the coming years, Comet Interceptor will make yet another attempt. This time, the target comet and its outgassing activity are unknown and may not be known before the spacecraft has been launched into its parking orbit, where it will await a possible interception. If the approximate outgassing rate can be estimated remotely when a target has been identified, it is desirable to also be able to estimate the scale size of the plasma environment, defined here as the region bound by the bow shock.

Aims. This study aims to combine previous measurements and simulations of cometary bow shock locations to gain a better understanding of how the scale size of cometary plasma environments varies. We compare these data with models of the bow shock size, and we furthermore provide an outgassing rate-dependent shape model of the bow shock. We then use this to predict a range of times and cometocentric distances for the crossing of the bow shock by Comet Interceptor, together with expected plasma density measurements along the spacecraft track.

Methods. We used data of the location of cometary bow shocks from previous spacecraft missions, together with simulation results from previously published studies. We compared these results with an existing model of the bow shock stand-off distance and expand on this to provide a shape model of cometary bow shocks. The model in particular includes the cometary outgassing rate, but also upstream solar wind conditions, ionisation rates, and the neutral flow velocity.

Results. The agreement between the gas-dynamic model and the data and simulation results is good in terms of the stand-off distance of the bow shock as a function of the outgassing rate. For outgassing rates in the range of 1027-1031-s-1, the scale size of cometary bow shocks can vary by four orders of magnitude, from about 102 km to 106 km, for an ionisation rate, flow velocity, and upstream solar wind conditions typical of those at 1 AU. The proposed bow shock shape model shows that a comet plasma environment can range in scale size from the plasma environment of Mars to about half of that of Saturn.

Conclusions. The model-data agreement allows for the planning of upcoming spacecraft comet encounters, such as that of Comet Interceptor, when a target has been identified and its outgassing rate is determined. We conclude that the time a spacecraft can spend within the plasma environment during a flyby can range from minutes to days, depending on the comet that is visited and on the flyby speed. However, to capture most of the comet plasma environment, including pick-up ions and upstream plasma waves, and to ensure the highest possible scientific return, measurements should still start well upstream of the expected bow shock location. From the plasma perspective, the selected target should preferably be an active comet with the lowest possible flyby velocity.

Here we describe the novel, multi-point Comet Interceptor mission. It is dedicated to theexploration of a little-processed long-period comet, possibly entering the inner Solar System for the first time, or to encounter an interstellar object originating at another star. Theobjectives of the mission are to address the following questions: What are the surface composition, shape, morphology, and structure of the target object? What is the composition ofthe gas and dust in the coma, its connection to the nucleus, and the nature of its interactionwith the solar wind? The mission was proposed to the European Space Agency in 2018, andformally adopted by the agency in June 2022, for launch in 2029 together with the Arielmission. Comet Interceptor will take advantage of the opportunity presented by ESA’s FClass call for fast, flexible, low-cost missions to which it was proposed. The call requireda launch to a halo orbit around the Sun-Earth L2 point. The mission can take advantage ofthis placement to wait for the discovery of a suitable comet reachable with its minimum ΔV capability of 600 ms⁻¹. Comet Interceptor will be unique in encountering and studying, at anominal closest approach distance of 1000 km, a comet that represents a near-pristine sample of material from the formation of the Solar System. It will also add a capability that noprevious cometary mission has had, which is to deploy two sub-probes – B1, provided by theJapanese space agency, JAXA, and B2 – that will follow different trajectories through thecoma. While the main probe passes at a nominal 1000 km distance, probes B1 and B2 willfollow different chords through the coma at distances of 850 km and 400 km, respectively.The result will be unique, simultaneous, spatially resolved information of the 3-dimensionalproperties of the target comet and its interaction with the space environment. We presentthe mission’s science background leading to these objectives, as well as an overview of thescientific instruments, mission design, and schedule.

Magnetosheath jets represent localized enhancements in dynamic pressure observed within the magnetosheath. These energetic entities, carrying excess energy and momentum, can impact the magnetopause and disrupt the magnetosphere. Therefore, they play a vital role in coupling the solar wind and terrestrial magnetosphere. However, our understanding of the morphology and formation of these complex, transient events remains incomplete over two decades after their initial observation. Previous studies have relied on oversimplified assumptions, considering jets as elongated cylinders with dimensions ranging from $0.1\, R_{\rm E}$ to $5\, R_{\rm E}$ (Earth radii). In this study, we present simulation results obtained from Amitis, a high-performance hybrid-kinetic plasma framework (particle ions and fluid electrons) running in parallel on graphics processing units (GPUs) for fast and more environmentally friendly computation compared to CPU-based models. Considering realistic scales, we present the first global, three-dimensional (3D in both configuration and velocity spaces) hybrid-kinetic simulation results of the interaction between solar wind plasma and the Earth. Our high-resolution kinetic simulations reveal the 3D structure of magnetosheath jets, showing that jets are far from being simple cylinders. Instead, they exhibit intricate and highly interconnected structures with dynamic 3D characteristics. As they move through the magnetosheath, they wrinkle, fold, merge, and split in complex ways before a subset reaches the magnetopause.

Aims: We aim to quantify the width of the quasi-perpendicular Martian bow shock region to deepen the understanding of why the width is variable and which factors affect it, and to explore the implications on thermalization.

Methods: To quantify the width, 2074 quasi-perpendicular bow shock crossings from a database were studied. Upstream conditions, such as Mach numbers, dynamic pressure, ion densities, and other factors, were considered. Furthermore, the difference between the downstream and upstream temperature was measured.

Results: We found that the shock region width is correlated with the magnetosonic Mach number, the critical ratio, and the overshoot amplitude. The region was found to be anticorrelated with dynamic pressure. The width is not affected by the upstream ion density of the investigated species or by the upstream temperature. The difference between the downstream and upstream temperature is not affected by the shock region width.

Conclusions: We found that the factors that affect the stand-off distance of the bow shock, such as the magnetosonic Mach number and dynamic pressure, also affect the width. The width is also positively correlated with the overshoot amplitude, indicating that the structures are coupled or that they are affected by largely the same conditions. The lack of a correlation with the ion temperature difference indicates that the shock region width does not affect the ion thermalization.