# BepiColombo Science Investigations During Cruise and Flybys at the Earth, Venus and Mercury

Valeria Mangano, Melinda Dósa, Markus Fränz, Anna Milillo, Joana S. Oliveira, Yeon Joo Lee, Susan Mckenna-Lawlor, Davide Grassi, Daniel Heyner, Alexander S. Kozyrev, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Valeria Mangano, Melinda Dósa, Markus Fränz, Anna Milillo, Joana S. Oliveira, et al. BepiColombo Science Investigations During Cruise and Flybys at the Earth, Venus and Mercury. Space Science Reviews, 2021, 217, pp.23. 10.1007/s11214-021-00797-9. insu-03139759

HAL Id: insu-03139759

https://insu.hal.science/insu-03139759

Submitted on 12 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# BepiColombo Science Investigations During Cruise and Flybys at the Earth, Venus and Mercury

Valeria Mangano<fnref n="1" />, Melinda Dósa<fnref n="2" />, Markus Fränz<fnref n="3" />, Anna Milillo<fnref n="1" />,

Joana S. Oliveira<fnref n="4" />,<fnref n="5" />, Yeon Joo Lee<fnref n="6" />, Susan McKenna-Lawlor<fnref n="7" />, Davide Grassi<fnref n="1" />,

Daniel Heyner<fnref n="8" />, Alexander S. Kozyrev⁹, Roberto Peron<fnref n="1" />, Jörn Helbert¹⁰,

Sebastien Besse¹¹, Sara de la Fuente¹², Elsa Montagnon¹³, Joe Zender<fnref n="4" />,

Martin Volwerk¹⁴, Jean-Yves Chaufray¹⁵, James A. Slavin¹⁶, Harald Krüger<fnref n="3" />,

Alessandro Maturilli¹⁰, Thomas Cornet¹¹, Kazumasa Iwai¹⁷, Yoshizumi Miyoshi¹⁷,

Marco Lucente<fnref n="1" />, Stefano Massetti<fnref n="1" />, Carl A. Schmidt¹⁸, Chuanfei Dong¹⁹,

Francesco Quarati²⁰,²¹, Takayuki Hirai²², Ali Varsani¹⁴, Denis Belyaev²³,

Jun Zhong²⁴, Emilia K.J. Kilpua²⁵, Bernard V. Jackson²⁶, Dusan Odstrcil²⁷,²⁸,

Ferdinand Plaschke¹⁴, Rami Vainio²⁹, Riku Jarvinen³⁰,³¹,

Stavro Lambrov Ivanovski³², Ákos Madár<fnref n="2" />, Géza Erdős<fnref n="2" />, Christina Plainaki³³,

Tommaso Alberti<fnref n="1" />, Sae Aizawa³⁴, Johannes Benkhoff<fnref n="4" />, Go Murakami³⁵,

Eric Quemerais¹⁵, Harald Hiesinger³⁶, Igor G. Mitrofanov⁹, Luciano Iess³⁷,

Francesco Santoli<fnref n="1" />, Stefano Orsini<fnref n="1" />, Herbert Lichtenegger¹⁴, Gunther Laky¹⁴,

Stas Barabash³⁸, Richard Moissl³⁹, Juhani Huovelin²⁵, Yasumasa Kasaba⁴⁰,

Yoshifumi Saito³⁵, Masanori Kobayashi⁴¹, Wolfgang Baumjohann¹⁴

Received: 13 June 2020 / Accepted: 13 January 2021
© The Author(s) 2021

**Abstract** The dual spacecraft mission BepiColombo is the first joint mission between the European Space Agency (ESA) and the Japanese Aerospace Exploration Agency (JAXA) to explore the planet Mercury. BepiColombo was launched from Kourou (French Guiana) on

The BepiColombo mission to Mercury

Edited by Johannes Benkhoff, Go Murakami and Ayako Matsuoka

<figure><img src="image_3.png" /></figure>

V. Mangano

<footnote n="1">¹ INAF-IAPS (Institute for Astrophysics and Planetology from Space), via del Fosso del Cavaliere 100, 00133, Rome, Italy</footnote>

<footnote n="2">² Department of Space Research and Space Technology, Wigner Research Centre for Physics, Budapest, Hungary</footnote>

<footnote n="3">³ Max Planck Institute for Solar System Research, Goettingen, Germany</footnote>

<footnote n="4">⁴ Science Division, European Space Agency (ESA), European Space Research and Technology Centre (ESTEC), Noordwijk, The Netherlands</footnote>

<footnote n="5">⁵ Space Magnetism Area, Payloads and Space Sciences Department, INTA, Torrejón de Ardoz, Spain</footnote>

<footnote n="6">⁶ Zentrum für Astronomie und Astrophysik, Technische Universität, Berlin, Germany</footnote>

<footnote n="7">⁷ Space Technology Ireland, Ltd., Maynooth, Co. Kildare, Ireland</footnote>

<footnote n="8">⁸ Institut für Geophysik und extraterrestrische Physik, Technische Universität, Braunschweig, Germany</footnote>
October 20th, 2018, in its packed configuration including two spacecraft, a transfer module, and a sunshield. BepiColombo cruise trajectory is a long journey into the inner heliosphere, and it includes one flyby of the Earth (in April 2020), two of Venus (in October 2020 and August 2021), and six of Mercury (starting from 2021), before orbit insertion in December 2025. A big part of the mission instruments will be fully operational during the mission cruise phase, allowing unprecedented investigation of the different environments that will encounter during the 7-years long cruise. The present paper reviews all the planetary flybys and some interesting cruise configurations. Additional scientific research that will emerge in the coming years is also discussed, including the instruments that can contribute.

**Keywords** BepiColombo · Mercury · Venus · Earth · Cruise · Flyby

# 1 Introduction

The dual spacecraft mission BepiColombo is the first joint mission between the European Space Agency (ESA) and the Japanese Aerospace Exploration Agency (JAXA) to explore the planet Mercury. BepiColombo was launched from Kourou (French Guiana) on October 20th, 2018 in its packed configuration named Mercury Composite Spacecraft, MCS) including the two spacecraft: the ESA Mercury Planetary Orbiter (MPO) and the JAXA Mio (formerly named Mercury Magnetospheric Orbiter) (see Benkhoff et al., this journal). Hence, BepiColombo packed configuration (usually named Mercury Composite Spacecraft,

<footnote n="9">⁹ Institute for Space Research, Moscow, Russian Federation</footnote>

<footnote n="10">¹⁰ Deutsches Zentrum für Luft- und Raumfahrt (DLR), German Aerospace Center Institute for Planetary Research, Berlin, Germany</footnote>

<footnote n="11">¹¹ Aurora Technology B.V. for ESA, European Space Agency, Madrid, Spain</footnote>

<footnote n="12">¹² ESA/ESAC, Villanueva de la Canada, Spain</footnote>

<footnote n="13">¹³ ESA/ESOC, Darmstadt, Germany</footnote>

<footnote n="14">¹⁴ Space Research Institute, Austrian Academy of Sciences, Graz, Austria</footnote>

<footnote n="15">¹⁵ LATMOS, Université Versailles Saint-Quentin, Guyancourt, France</footnote>

<footnote n="16">¹⁶ Climate and Space Sciences and Engineering, University of Michigan, Ann Arbor, MI, USA</footnote>

<footnote n="17">¹⁷ Institute for Space-Earth Environmental Research, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, 464-8601, Japan</footnote>

<footnote n="18">¹⁸ Boston University, Boston, USA</footnote>

<footnote n="19">¹⁹ Department of Astrophysical Sciences and Princeton Plasma Physics Laboratory, Princeton University, Princeton, NJ, USA</footnote>

<footnote n="20">²⁰ LumMat/RST/AP, Delft University of Technology, Delft, The Netherlands</footnote>

<footnote n="21">²¹ Gonitec BV, Den Haag, The Netherlands</footnote>

<footnote n="22">²² Planetary Exploration Research Center (PERC), Chiba Institute of Technology, Chiba, Japan</footnote>

<footnote n="23">²³ Space Research Institute (IKI), Moscow, Russia</footnote>

<footnote n="24">²⁴ Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, China</footnote>

<footnote n="25">²⁵ Department of Physics, University of Helsinki, Helsinki, Finland</footnote>
MCS) is composed by ESA spacecraft MPO, the JAXA spacecraft Mio, the Mercury Transfer Module MTM, and the sunshade cone MOSIF (acronym for Magnetospheric Orbiter Sunshield and Interface Structure) designed to protect Mio during the cruise phase.

BepiColombo cruise trajectory is a long journey into the inner heliosphere, and it includes one flyby of the Earth (in April 2020), two of Venus (in October 2020 and August 2021), and six of Mercury itself (starting from 2021), before the orbit insertion in December 2025.

Figure 1 shows a schematic of the complete cruise trajectory (from launch to orbit insertion), and Fig. 2 shows the composite spacecraft MCS with its reference axes: +Y is directed into the radiator panel of MPO, +Z is along the main axis of MCS, from Mio towards the MTM; +X completes the Cartesian 3-axis system.

The MPO payload comprises 11 experiments and instrument suites, and Mio carries 5 experiments and instrument suites. The former will focus on the planet itself and close interactions with its surface, while the latter will study the environment around the planet including the planet's exosphere and magnetosphere, and their interaction processes with the solar wind. However, the MCS configuration will not allow full operability of all instruments onboard. In fact, Mio will be partly shielded by MOSIF, thus allowing instruments to detect signals only within a conical field-of-view around the MCS's -Z axis. On-board the MPO, the instruments with apertures in the +Z axis are obstructed by the MTM and, hence, are not able to acquire scientific measurements. Nevertheless, all the instruments not requiring pointing or with apertures in the other directions will be, in principle, able to operate.

Hence, BepiColombo cruise phase has become an important opportunity to investigate different environments of the inner Solar System, even if it was not originally designed for that purpose.

<footnote n="26">26 Center for Astrophysics and Space Sciences, University of California at San Diego, La Jolla, CA, USA</footnote>

<footnote n="27">27 Department of Physics and Astronomy, George Mason University, Fairfax, VA, USA</footnote>

<footnote n="28">28 NASA/GSFC Heliophysics Science Division, Space Weather Laboratory, Greenbelt, MD, USA</footnote>

<footnote n="29">29 Department of Physics and Astronomy, University of Turku, Turku, Finland</footnote>

<footnote n="30">30 Finnish Meteorological Institute, Helsinki, Finland</footnote>

<footnote n="31">31 Department of Electronics and Nanoengineering, School of Electrical Engineering, Aalto University, Espoo, Finland</footnote>

<footnote n="32">32 INAF-OATs (Osservatorio Astronomico di Trieste), Trieste, Italy</footnote>

<footnote n="33">33 Italian Space Agency, Rome, Italy</footnote>

<footnote n="34">34 Institut de Recherche en Astrophysique et Planétologie, CNRS, CNES, Université de Toulouse, Toulouse, France</footnote>

<footnote n="35">35 Japan Aerospace Exploration Agency, Institute of Space and Astronautical Science, Sagamihara, Kanagawa 252-5210, Japan</footnote>

<footnote n="36">36 Institut für Planetologie, Westfälische Wilhelms-Universität Münster Münster, Germany</footnote>

<footnote n="37">37 Dipartimento di ingegneria meccanica e aerospaziale, Università La Sapienza, Roma, Italy</footnote>

<footnote n="38">38 Swedish Institute of Space Physics, Kiruna, Sweden</footnote>

<footnote n="39">39 Research and Scientific Support Department, ESTEC, Noordwijk, The Netherlands</footnote>

<footnote n="40">40 Planetary Plasma and Atmospheric Research Center, Graduate School of Science, Tohoku University, Sendai, Miyagi 980-8578, Japan</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 1** BepiColombo spacecraft sketch of different cruise phases: launch in 2018, MCS orbit around the Earth, Earth's flyby, cruise, Venus flybys, Mercury's flybys, separation, and orbit insertion of the MPO and Mio spacecraft in 2026 after separation of the 4 components of MCS (courtesy of Airbus Defense and Space)</figcaption></figure>

The present paper aims to review all the additional science investigations that may take place during the planetary flybys and during the cruise in specific configurations with other spacecraft.

Section 2 will discuss the cruise phase, including attitude details and limitations of the MCS configuration, and focus on the possible scientific objectives to be addressed during the cruise phase in general. It will include discussions on multi-spacecraft coordinated observations, and a brief description of how the present space weather modeling tools may help in predicting and interpreting the data.

Section 3 will focus on the radio science experiments which are the only ones originally planned during the cruise of BepiColombo. Section 4 will report on the Earth flyby investigations that occurred on April 10th 2020. Section 5 will give an overview of the planet Venus characteristics in terms of atmosphere and induced magnetosphere, and will describe in detail the possible Venus investigations during the two flybys with the instrumentation onboard MCS in operation. Finally, Section 6 will focus on the Mercury flybys, from October 2021 to January 2025, being with the first approach to the planet. The primary goal of the BepiColombo mission starting from 2026 is to study Mercury and its environment, when the two spacecraft will be in their own orbits, with the instrumentation onboard fully operative. Nevertheless, the peculiar geometry of the flybys will allow first investigations of Mercury exosphere, intrinsic magnetic field and magnetosphere, and also surface and interior starting from 2021, i.e. 5 years in advance with respect to the nominal mission.

<footnote n="41">41 Planetary Exploration Research Center, Chiba Institute of Technology, 2-17-1, Tsudanuma, Narashino, Chiba 275-0016, Japan</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig. 2** The composite spacecraft MCS during cruise and its reference frame *XYZ*: *Y*-axis points towards the Sun (and MPO radiator – not visible in this image – points towards –*Y* axis); +*Z*-axis is along the main axis of the spacecraft pointing at the anti-MTM direction (i.e., Mio – or MMO – axis is pointing towards –*Z*), and +*X* axis completes the right hand coordinate system. Note that the magnetometer boom of MAG and the antenna used by MORE are both fully deployed during cruise (courtesy of Airbus Defense and Space)</figcaption></figure>

Section 7 will summarize the different scientific tasks that BepiColombo may contribute to investigate during its cruise to Mercury, together with our conclusions. An Appendix will complement the present paper, to provide the reader with other useful information related to the wide field of investigations discussed.

## 2 Cruise

During the interplanetary cruise phase, the radiator panel of MPO (−*Y* axis) is pointing anti-sunwards (see Fig. 2), thus +*Y* axis is directed towards the Sun. Only small deviations from this pointing described by solar aspect angle (SAA), that is the angle between +*Y* axis and the Sun direction projected onto the *XZ* plane are allowed during the cruise phase. The allowed deviations are between +47° and +17° depending on the spacecraft distance to the Sun (see Table 9 in the Appendix).

Orbit changes are achieved by the planetary flyby maneuvers and by firing of the Solar Electric Propulsion System into specific periods (see Table 10 in the Appendix). The spacecraft attitude is controlled around the Sun line such that the angular momentum loading is minimized while ground contact is maximized during the ground station passes. This results in tight attitude constraints, as well as limited communication windows with the High Gain Antenna, typically only up to a maximum of 10 hours per day.

Hence, during the cruise phase, unobstructed instrument operations are mainly limited by the following constraints:

1) as a baseline, no instrument operations are planned during solar electric propulsion, both for MPO and Mio;

2) as BepiColombo moves away from the Earth, flexibility in attitude and the downlink capability will further limit possible instrument operations.

Since the launch, and outside the solar electric propulsion periods, onboard the MPO, the magnetometer MAG (Glassmeier et al. 2010), the Italian Spring Accelerometer ISA (Santoli
et al. 2020), the radio science experiment MORE (Iess et al. 2021, this journal), the Mercurian Gamma-ray and Neutron Spectrometer MGNS (Mitrofanov et al. 2010; Mitrofanov et al. 2021, this journal), and the environmental radiation sensor BERM (Moissl et al., this journal) are fully operative.

The X-ray and particle spectrometers of SIXS (Huovelin et al. 2020, this journal) and the UV-spectrometer PHEBUS (Quémerais et al. 2020, this journal) will be operated throughout the cruise phase depending on dedicated feasibility analysis.

The magnetometer MAG onboard MPO is mounted on the fully deployed 2.8 m boom at the edge of the radiator (−Y axis) and the X axis side, and the antenna used by MORE is fully deployed too (see Fig. 2). The MGNS spectrometers are located inside the MPO spacecraft structure and measures gamma-rays and neutrons, however during cruise these measurements are impacted by the surrounding material of the composite stack. The ISA accelerometer is also inside the spacecraft structure but does not need pointing and is fully operative.

The BERM sensor is mounted behind the radiator panel and looks to the −Y direction, thus always anti-sunward. The three SIXS-X X-ray detectors are oriented towards the +Y direction (unfortunately with an offset such that none of the detectors has the Sun in its field-of-view during the cruise phase). The SIXS-P particle detector is operational and views in five orthogonal directions (three of them not obstructed during cruise). The PHEBUS baffle is located at the radiator panel, and can rotate 360° around the −Y-axis. The infrared spectrometer MERTIS (Hiesinger et al. 2020, this journal) can be operated by its calibration telescope also pointing in −Y-direction. Finally, the two ion sensors of SERENA are on the edge of the −X axis, and can be operated too (the other two neutral sensors in the +Z face being blocked by the MTM). Mio instrumentation is only partially obstructed by the sunshield, and have a free field of view of about 30 degrees in the −Z direction (Murakami et al. 2020). MPPE sensors (HEP-e, ENA, MSA, MIA, and MEA) (Saito et al. 2021, this journal) and MDM (Kobayashi et al. 2020) have some open field of views around −Z axis above MOSIF. The MGF (Baumjohann et al. 2020) and PWI (Kasaba et al. 2020; Yagitani et al. 2020; Karlsson et al. 2020) instruments are not deployed yet; hence, noise levels are higher during the cruise phase than during the Mercury orbit phase. Table 1 summarizes the instruments viewing directions, field of views and constraints valid during the cruise phase.

## 2.1 Scientific Objectives

During the long cruise of BepiColombo in the inner Solar System, the only official science goal originally planned was the superior solar conjunction experiment to be performed by MORE (see Sect. 3 and Iess et al. 2021, this journal). Nevertheless, in recent years it was clear that the opportunity to operate the other available instruments should be exploited.

During the BepiColombo cruise phase (2018–2025), the solar cycle will be in its ascending phase with predicted solar maximum predicted in July 2025. This fact will provide the opportunity to study interplanetary physics under different conditions of solar activity (Fig. 3).

Even limited instrument operations of BepiColombo can contribute to a wide range of scientific investigations, as for example:

- In-situ solar wind observations on small scales and large scales, with special regard to the study of transient events such as Coronal Mass Ejections (CMEs), and Solar Energetic Particles (SEPs), but also background solar wind and scale-invariant properties (turbulence, intermittency, MHD instabilities), high speed streams, Co-rotating Interaction Regions (CIRs), and Heliospheric Current Sheet (HSC).
<table><caption>**Table 1** Viewing directions of the instruments and their constraints during the cruise phase. Viewing directions are given in the BepiColombo spacecraft stacked configuration, as given in Fig. 2. (*) for more details refer to Fig. 4 in Murakami et al. (2020, this journal)</caption><thead><tr><th>Host</th><th>Instrument</th><th>Location</th><th>Viewing direction (x, y, z)</th><th>Viewing constraints</th></tr></thead><tbody><tr><td rowspan="10">MIO (*)</td><td>MDM</td><td>on Mio side panel</td><td>(0.0000, 0.0000, -1.0000)</td><td>FOV partially obstructed by MOSIF, 30° FOV free</td></tr><tr><td>MGF</td><td>on Mio side panel</td><td>N/A</td><td>not deployed yet</td></tr><tr><td>MPPE-HEP-e</td><td>on Mio side panel</td><td>(0.2615, -0.4556, -0.8509)</td><td>FOV partially obstructed by MOSIF, 15° FOV free</td></tr><tr><td>MPPE-ENA</td><td>on Mio side panel</td><td>(-0.4099, -0.2411, -0.8789)</td><td>FOV partially obstructed by MOSIF, 20° FOV free</td></tr><tr><td>MPPE-MSA</td><td>on Mio side panel</td><td>(-0.0135, 0.0931, -0.9956)</td><td>FOV partially obstructed by MOSIF, 45° FOV free</td></tr><tr><td>MPPE-MIA</td><td>on Mio side panel</td><td>(0.0942, -0.1312, -0.9869)</td><td>FOV partially obstructed by MOSIF, 45° FOV free</td></tr><tr><td>MPPE-MEA1</td><td>on Mio side panel</td><td>(-0.0150, -0.0474, -0.9988)</td><td>FOV partially obstructed by MOSIF, 45° FOV free</td></tr><tr><td>MPPE-MEA2</td><td>on Mio side panel</td><td>(-0.0236, -0.1411, -0.9897)</td><td>FOV partially obstructed by MOSIF, 45° FOV free</td></tr><tr><td>PWI</td><td>on Mio side panel</td><td>N/A</td><td>not deployed yet</td></tr><tr><td rowspan="5">MPO</td><td>BERM</td><td>behind radiator panel</td><td>N/A</td><td>mounted in radiator panel with boresight is in -Y</td></tr><tr><td>ISA</td><td>inside the MPO</td><td>N/A</td><td>mounted inside the MPO, no constraints</td></tr><tr><td>MPO-MAG</td><td>top-edge of radiator panel</td><td>(1.0000, 0.0000, 0.0000)</td><td>deployed in October 2018, no constraints</td></tr><tr><td>MERTIS-SPACE</td><td>inside MPO on +Z direction</td><td>(0.0000, -1.0000, 0.0000)</td><td>Planet view port is in Z-direction, thus obstructed</td></tr><tr><td>MGNS</td><td>inside the MPO</td><td>N/A</td><td>no preferred viewing direction.</td></tr><tr><td rowspan="5">MTM</td><td>MORE</td><td>N/A</td><td>N/A</td><td>Measurements are impacted by surrounded material</td></tr><tr><td>PHEBUS</td><td>on radiator panel</td><td>x-z-plane</td><td>high-gain, medium-gain, and low-gain antenna deployed</td></tr><tr><td>SERENA-PICAM</td><td>-X-edge of radiator panel</td><td>(-1.0000, 0.0000, 0.0000)</td><td>rotates 360°</td></tr><tr><td>SERENA-MIPA</td><td>-X-edge of radiator panel</td><td>(-1.0000, 0.0000, 0.0000)</td><td>1.5 π steradian FOV with boresight in -X</td></tr><tr><td>SIXS-X</td><td>top-edge of radiator panel</td><td>(0.0000, 1.0000, 0.0000)</td><td>2 π steradian FOV with boresight in -X</td></tr><tr><td rowspan="3">MCAM1</td><td>SIXS-P</td><td>top-edge of radiator panel</td><td>(0.0000, 1.0000, 0.0000)</td><td>Three X-ray detectors oriented in the +Y direction, with an offset such that none of the detectors has the Sun in its field-of-view during the cruise phase.</td></tr><tr><td>MCAM2</td><td>-X panel</td><td>(-0.9811, 0.0858, -0.1736)</td><td>Four of the five particle detectors are oriented in the y-x-plane, with 90° offsets to each other</td></tr><tr><td>MCAM3</td><td>-X panel</td><td>(-0.5826, -0.2717, -0.7660)</td><td>The fifth detector is oriented towards -Z</td></tr><tr><td></td><td></td><td></td><td>(0.4545, -0.4545, -0.7660)</td><td>no constraints<br>no constraints<br>no constraints</td></tr></tbody></table>

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 3** Solar cycle 25 forecast update compared to previous ones as released on December 2019 by NOAA. The rising phase of cycle 25 is expected to peak in July 2025 (courtesy of NOAA)</figcaption></figure>

- Measurements of the plasma composition, ion flux and density, and magnetic fields in the vicinity of Venus and Mercury during flybys (see devoted Sects. 5 and 6).

It has to be emphasized, that while MPO and Mio will have polar orbits around Mercury, the flybys will occur at low-latitudes. These are the best opportunities to study low-latitude solar wind interaction with the magnetosphere, Kelvin-Helmholtz boundary waves and important internal dynamics like plasma sheet convection and reconnection.

- coordinated science observations together with other spacecraft in Earth orbit and in the inner Solar System, e.g., Akatsuki, Parker Solar Probe, SDO, Proba-2, Hisaki, SOHO, Solar Orbiter and JUICE, to provide measurements from additional vantage points;

- observations of InterPlanetary Scintillations (IPS) coordinated with Earth-based measurements to provide information on solar wind density;

- analysis of cometary composition, and detection of dust particles of different origins;

- monitoring the local radiation background due to bombardment by energetic particles of Galactic Cosmic Rays;

- Gamma Ray Burst (GRB) detection and localization;

- and, of course, the above mentioned Superior Solar Conjunction measurements (SSC) to test general relativity.

### 2.1.1 Solar Wind and Interplanetary Magnetic Field

The solar wind is a continuous flow of charged particles that originates from the Sun's corona. The magnetic energy contained at lower regions rises to the higher corona and dissipates via processes not yet understood in detail, causing coronal heating and solar wind acceleration (e.g. Shoda et al. 2019). As solar wind propagates radially outward, it fills up the heliosphere and interacts with various obstacles (e.g. planetary bodies). The parameters of solar wind plasma and its dynamic changes determine all plasma processes in space. Also the Interplanetary Magnetic Field (IMF) carried by the solar wind is prone to various kinds of structures and discontinuities on a broad range of scales (see e.g. Burlaga et al. 1969) from the large scale ones, like CMEs and CIRs, to small scale structures like Magnetic Holes (MHs) and Mirror Mode (MM) waves, via middle scale reconnection structures and flux tubes.

In particular, BepiColombo can provide a huge dataset of the interplanetary magnetic field through the almost continuous measurements of MPO-MAG, plus localized (in time and space) measurements of the solar wind parameters like density and velocity by other ion sensors, and radiation at other wavelengths through the X- and gamma-sensors onboard.
**Small Scale Heliospheric Structures** Magnetic holes (MHs) were first detected in the Explorer 34 (IMP I) data by Turner et al. (1977), and identified as depressions in the solar wind magnetic field strength with $|B| < 1$ nT in a background field around 6 nT and a very small rotation of the magnetic field over the structure ($\theta < 10^\circ$). MHs were assumed to be diamagnetic, created by inhomogeneities in the solar wind, but the plasma data for these structures were insufficient to actually prove this assumption. Later studies (e.g. Stevens and Kasper 2007) have shown that MHs mainly occur in high solar wind plasma-$\beta$ regions, which makes them similar to Mirror Mode (MM) waves (Hasegawa 1969). MMs occur in high-$\beta$ plasmas when there is a temperature asymmetry. This MM instability usually leads to non-propagating wave trains of structures where the magnetic field and plasma density oscillate in anti-phase, which are transported by the plasma flow. The MHs seemed to occur mainly in MM stable regions of the solar wind (Stevens and Kasper 2007) and wave trains mainly appear in MM unstable plasma regions (Winterhalter et al. 1995). This may lead to the hypothesis that MHs are the final stages of MMs. In addition, Hasegawa and Tsurutani (2011) proposed a Bohm-like diffusion mechanism (Bohm and Gross 1949), in which the higher frequencies of the structure diffuse faster than the lower ones, which results in larger magnetic structures. This mechanism was shown to happen also in Venus's magnetosheath and in comet 1P/Halley's coma (Schmid et al. 2014).

Recent investigations of linear MH occurrence rates using MESSENGER magnetometer data (Anderson et al. 2007) during the cruise from Venus to Mercury in 2007–2011 were performed (Volwerk et al. 2020). The results show a gradual decrease in the occurrence rate of linear MHs from Mercury (~14%) to Venus (~10%). This agrees with the observation that MHs in the solar wind occur between ~0.1 to ~7.5 per day at different locations in the solar system. Hence, an analysis throughout the cruise, in particular before and after Venus flybys, would be of particular interest to better study the magnetic holes characteristics. With the MPO-MAG boom deployed and the instrument operational, the detection and measurements of MHs can be repeated through BepiColombo's cruise phase, as done during the MESSENGER cruise phase (years 2004–2011).

In particular, this type of measurements by BepiColombo may fill the gap of magnetic data at heliospheric distances between Earth and Venus that MESSENGER did not acquire. In addition, magnetic field measurements will enable the study of, e.g., the solar wind turbulence (such as its evolution with the heliocentric distance), the energy transfer from large to small scales (by looking at the behavior of magnetic field spectrum, e.g. Korth et al. 2011), the role of stochastic fluctuations, magnetic reconnection, intermittency, self-organized structures and waves on spectral and scaling properties (e.g. Dong et al. 2018b). Concerning the radial evolution of solar wind turbulence we will understand more deeply the role of non-linearities during the expansion phase of the solar wind (Tu and Marsch 1995), i.e., a key feature for properly characterizing energy transfer across different scales and at different heliocentric distances (Alberti et al. 2020; Chen et al. 2020). Moreover, we will also have the opportunity to investigate small-scale processes (e.g., sub-ion and ion-kinetic effects) as well as some characteristic phenomena of turbulent flows as coherent structures, i.e., plasma regions of concentrated vorticity, intermittency, i.e., the manifestation of sudden field changes manifesting in field gradients and anomalous scalings of the field increments, and the validity of an exact relation for the scaling of the third-order moment of fluctuations, e.g., the Yaglom law, at various heliocentric distances (e.g., Bruno and Carbone 2016).

**Large Scale Heliospheric Structures: Transients** The understanding of the properties and evolution of large scale interplanetary structures is one of the major challenges in current heliospheric research. This is largely due to the fact that only limited observations have been
available with varying heliospheric distances as most measurements come at or near the Earth orbit, plus the Helios missions in the 1970 and occasional measurements by planetary missions, such as MESSENGER and Venus Express (Schwenn 1990; Bothmer and Schwenn 1998; Winslow et al. 2015; Good and Forsyth 2016).

Particular attention should be put on:

- Interplanetary Coronal Mass Ejections (ICMEs). ICMEs are the interplanetary counterparts of coronal mass ejections (CMEs) and they are strongly distorted by their interaction with the background solar wind which determines significantly their evolution. The overall structure of the resulting interplanetary disturbance is very complex. ICMEs are important drivers of space weather events not only in the circumterrestrial space but also at planetary environments (see discussions in Lilensten et al. 2014; Plainaki et al. 2016; Kilpua et al. 2017, and references therein).

The integral part of the ICME is a magnetic flux rope that consists of bundles of helical magnetic field lines that wind about the common axis (e.g., Klein and Burlaga 1982; Burlaga 1988; Exner et al. 2018). When sufficiently faster than the ambient solar wind, ICMEs drive fast forward leading shocks through the solar wind and a turbulent sheath forms between the shock and the ICME leading edge. ICME-driven sheaths are compressive structures, characterized by large amplitude magnetic field variations and large dynamic pressure. ICME flux ropes are associated with smoothly changing magnetic field direction and lower dynamic pressure. The global morphology of ICMEs and the nature of interactions between multiple ICMEs and between ICMEs and SIRs (Stream Interaction Regions) is still unresolved (e.g., whether they are generally inelastic, elastic or super-elastic and how the structure and kinematics of the ICMEs involved are affected, see Lugaz et al. 2017). Such issues could benefit from multi-point observations (see Sect. 2.2.1). Sometimes interactions can occur relatively quickly e.g., as featured by Kilpua et al. (2019) who investigated a case where two ICMEs were just about to start interacting at Venus, but by reaching the Earth had coalesced together, resembling a single coherent ICME flux rope, and the shock of the second ICME had also propagated past the preceding ICME.

Another open question is how CMEs that start as more or less coherent flux ropes at the Sun, often transform into quite complex ICMEs in interplanetary space (e.g., Dasso et al. 2007; Kilpua et al. 2013). On average, only about one-third of ICMEs near Earth orbit exhibit clear flux rope signatures (e.g., Cane and Richardson 2003), while the rest have more disorganized field characteristics.

The lack of flux rope signatures can be explained either by interplanetary evolution, e.g. magnetic flux erosion (e.g., Ruffenach et al. 2015), or by the spacecraft's crossing the ICME far from the center. Another possible reason can be that interactions occur with other large scale heliospheric structures (e.g., Manchester et al. 2017; Winslow et al. 2016). This is an additional reason why more measurements are needed.

In fact, observations made at the varying heliospheric distances may help clarify the physical mechanisms of the evolution/interactions of ICMEs and how they progressively change their structure. Since the BepiColombo cruise phase coincides with the ascending and maximum phases of Solar Cycle 25, a high number of ICMEs (though not predictable) are expected to be encountered.

- Solar Energetic Particles (SEPs). SEPs are mainly protons of solar origin and, to a lesser extent, heavier ions; they are accelerated to energies up to hundreds of MeV (and in some rare cases also GeVs) (Miroshnichenko 2018), either at the flare reconnection sites or at the shock regions driven by CMEs propagating through the solar corona in the interplanetary space (e.g., Desai and Giacalone 2016, and references therein). SEPs have been
observed at Mercury, Earth, and Mars and are important manifestations of space weather processes.

In the most extreme cases, they can fill the whole inner heliosphere with greatly enhanced fluxes of energetic electrons, protons and heavy ions (e.g., Dresing et al. 2014), that can last from days up to weeks. There are still several key open questions on their acceleration and transport that the inner heliospheric probes Solar Orbiter (ESA) and Solar Parker Probe (NASA) can help to answer, as it was done recently, when the combination of ACE, MESSENGER, STEREO-A, B, MAVEN, and Juno data at widely separated heliospheric longitudes was used to study:

(a) the cross-field diffusion processes in both the solar corona and interplanetary space,

(b) broad particle sources associated with coronal and interplanetary shocks initially driven by CMEs, and

(c) complex magnetic field configurations in the corona or interplanetary space that allow SEPs injected from a narrow source to reach distant locations (Lario et al. 2016).

In particular, comparison of SEP flux time-profiles and variations of the energy spectrum at different locations can provide crucial knowledge on whether the shocks are able to accelerate and inject SEPs into different longitudes, and whether cross-field diffusion causes wide spreading SEPs at different heliospheric longitudes. In addition, the combined studies of in-situ data, type-II radio burst data (e.g., Iwai et al. 2020), and global MHD simulations may provide clues on the structures of the background interplanetary magnetic fields and propagation characteristics of CMEs. Such clues are important to specify the acceleration regions of SEPs, and the 3D field line connection close to the solar surface from any specified in-situ locations including planets or a spacecraft.

## Large Scale Heliospheric Structures: Background Solar Wind

Solar wind outflow is continuous, and plasma velocity is usually defined either as “slow” (~400 km/sec) or “fast” (~700 km/s), even though showing wide range of values. In regions where slow and fast solar wind meet, solar wind SIRs (Richardson 2018) are formed. SIRs develop with radial distance, either the fast wind overtaking slow wind (compression region), or slow wind lagging behind fast wind (rarefaction region).

SIRs are also often named Co-rotating Interaction Regions (CIRs) to reflect their repeating nature with solar rotation on 27 days as well as on short-term harmonics (7 days, 13.5 days), especially during solar minimum conditions.

A typical SIR fingerprint shows a sudden and well-marked increase in density, temperature, and magnetic field strength and a constant change of velocity over the interface. The azimuthal component of the magnetic field has an opposite sign on each side of the interface, such as the non-radial components of the velocity vector.

Long-term observations of SIRs as they co-rotate with the Sun give us insight into coronal properties as well. The source of CIRs lies in the corona on the border of slow and fast solar wind sources. The study of typical CIR features by Dósa and Erdős (2017), e.g., alternating deviation of the velocity vector, magnetic flux density enhancements, etc., during the 23rd solar cycle, evidenced that during the calmer decreasing phase the source of the emerging CIRs were located on constant heliographic longitudes. Two exceptionally strong features during 2007–2008 are visible in Fig. 4. They appear in the magnetic field dataset as inclined red stripes around Carrington longitudes –180° and –60°, and also in the velocity deviation dataset around 180° and 300°. The coronal source longitude on the other hand rotated with a different rotation rate than the Carrington rotation (faster during the declining phase), but showed a constant trend, giving us some hints of coronal mechanisms. This
**Fig. 4** Source longitude of CIRs over the long term: unsigned magnetic flux density and solar wind azimuth angle (the angle between the velocity vector and the radial direction) measured by ACE during cycle 23 mapped back to the longitude of origin. Each horizontal line corresponds to 25.3802 days of data (time period corresponding to the Carrington rotation). Year to the right, Carrington rotation number to the left (y axis). Inclined red stripes are visible, implying that magnetic field enhancements (occurring at the same time as the deviation changes from positive to negative) recurred later or earlier than the Carrington rotation rate (from Dósa and Erdős 2017)

<figure><img src="image_8.png" /></figure>

finding should be further investigated in relation to active longitudes identified on the photosphere. Propagation properties of CIRs can be best studied when observed at different radial distances: especially the relationship between the magnetic sector boundaries and SIRs, but also the 3-dimensional evolution of the SIR feature.

Interfaces of slow and fast solar wind streams are often, but not necessarily located at magnetic sector boundaries. Dósa and Erdős (2017) also showed that during the years 2002–2006 the Heliospheric Current Sheet (identified as magnetic sector boundaries at 1 AU near the Earth) rotated also somewhat faster than the Carrington rate (25.3802 days per sidereal rotation) but with a different rate than the CIRs. Examining the rotation speed of the different features at different radial distances could enhance our understanding of how these features evolve in time and space, and what is their exact relationship with each other.

### 2.1.2 Dust and Comets

The study of dust distribution in the inner Solar System is important to understand the evolution of interplanetary dust travelling towards the Sun due to solar radiation drag (Poynting-Robertson effect), as well as a better constraint on the inner Solar System environment, both chemically and dynamically. Very few observations of dust in the inner Solar System are currently available. In the 1970s the Helios spacecraft observed the interplanetary dust distribution between 0.3 AU and 1.0 AU from the Sun for the first time (Grün et al. 1980). Analysis of these data led to the discovery that an interstellar dust stream originated from outside the Solar System and penetrates into the inner heliospheric region (Altobelli et al. 2006). Only recently, the JAXA’s solar power sail demonstrator IKAROS provided new data by its large-area dust impact sensor, ALADDIN (Yano and Hirai 2016), and revealed a 10 µm-sized interplanetary dust distribution from 1.1 AU to 0.72 AU (at Venus orbit). It is
<figure><img src="image_9.png" /><figcaption>**Fig. 5** Field of view of MDM during the cruise phase. The unobstructed field of view will be of 30° at maximum, and it derives from the relative positions of MDM, Mio, and MOSIF</figcaption></figure>

noteworthy that on 2nd September 2019, a comet-like object 322P/SOHO passed its perihe-
lion flying to a heliocentric distance of 0.12 AU, and swept by Parker Solar Probe (PSP) at a
relative distance as close as 0.025 AU, where images photographed by WISPR (Wide-field
Imager for Solar Probe Plus) onboard PSP shows messy trails caused by dust bombardment
(He et al. 2020). The 7-year cruise phase of BepiColombo, hence, is a precious opportunity
for in-situ measurement of cosmic dust including cometary dust trail particles in the inner
Solar System. The dust impact sensor MDM on a side panel of the Mio spacecraft will be
able to provide measurements of dust particles with diameter of ~ 1 µm (at an impact ve-
locity of 10 km/s and dust material density of 2.0 g/cm³) (Kobayashi et al. 2020). In fact,
though it is surrounded by the MOSIF (see Figs. 2 and 5), there is a possibility to detect dust
impacts through the MOSIF’s topside aperture which partially provides the MDM sensor
with ~ 30° field of view (see Fig. 5).

In addition, the MDM sensor, which is made of piezoelectric ceramics, can proba-
bly detect also acoustic vibrations propagating through the MOSIF wall and possibly the
whole BepiColombo spacecraft, induced by impacts of relatively large dust particles, e.g.
> 100 µm in diameter. For instance, the spacecraft passage through a comet dust trail (con-
sisting mainly of large millimeter-sized dust particles emitted from the parent comet), will
allow MDM to measure the small vibration of the spacecraft structure generated by the im-
pacts. Capitalising on its high sensitivity, also the ISA accelerometer can likely contribute
to this research topic. Although not directly exposed to outer space, ISA working principle
would allow detecting eventual structural vibrations induced by particle/dust impacts prop-
agating along the spacecraft. In this perspective, a synergy with MDM instrument can be
envisaged to support the identification of eventual impact events (see Sect. 3.2).

Simulations by the new universal model IMEX (Interplanetary Meteoroid Environment for eXploration) of the recently created cometary dust trails of different comets in the inner Solar System (Soja et al. 2015a,b) may help to study cometary dust trails, as measured by in-situ detectors and observed in infrared images, and program in advance devoted observations in particular periods of the BepiColombo cruise. Large variations in the predicted dust fluxes from comet to comet have to be expected because the ejection velocity, mass distribution, and dust production rate – all parameters of the IMEX model – likely vary for each comet and are not well constrained for many comets yet. This may be improved in the future for the comets found in the present analysis to yield more reliable flux predictions. Figure 6 shows the expected dust flux simulated with IMEX for cometary trails crossed by BepiColombo during the cruise. By far the largest flux is encountered in the trail of comet 2P/Encke. Fluxes predicted for the crossings of other cometary trails are at least two orders of magnitude lower.
<figure><img src="image_10.png" /><figcaption>**Fig. 6** Left: Simulated fluxes of particles larger than 100 µm for cometary trails crossed by the BepiColombo spacecraft during its interplanetary cruise to Mercury. The comet 2P/Encke (red square) shows the largest flux after 2022 when the spacecraft orbits the Sun inside Venus' orbit. The simulations were performed with a step size of two days. Right: BepiColombo interplanetary trajectory from 2021 to 2025 (blue). The orbit of comet 2P/Encke is shown in red; diamonds indicate BepiColombo's crossings of comet Encke's dust trail as derived from the IMEX model, and red lines attached to the diamonds show the approach direction (speed vector) of trail particles on to the spacecraft. The X–Y plane is the ecliptic plane with vernal equinox oriented towards the +X direction</figcaption></figure>

The flux peaks are rather narrow in space/time with a typical peak width of only several days (simulations in Fig. 6 have a two-day step size). On the right side of Fig. 6, the geometry of BepiColombo's orbit and the orbits of Encke are shown. The impact directions of Encke trail particles on to the BepiColombo spacecraft are indicated for a few crossings of this comet's trail. Maximum fluxes of 1 event/m2/day are reached for trail crossings in October 2023 and February 2024. These fluxes should be detectable with the MDM dust detector. In addition, during the period of maximum flux contemporary observation for the comet's hydrogen coma can be possibly performed by the PHEBUS instrument (see below and in Fig. 5). If this will be confirmed, it would be the first coordinated observation of in-situ cometary dust trail and coma observation in the inner Solar System.

The study of cometary composition is an important goal that could be addressed by Bepi-
Colombo during its cruise phase. One or two comets per year are visible in the inner helio-
sphere, on average. Hence, a total of about 10 comets could be observable during the 7-year
cruise to Mercury. This is an observation goal during the cruise of the UV spectrometer PHE-
BUS. For each target, it will be possible to compute the time when it intersects the PHEBUS
instrument FOV without a change of the spacecraft attitude. Depending on the brightness of
the object and its distance to the spacecraft, observations of the hydrogen coma at Lyman-α
(EUV) may be done, as done regularly by SWAN/SOHO (Bertaux et al. 1995), to derive the
water production rates (e.g., Combi et al. 2018). If possible, OH (FUV line at 308 nm) and
other weaker emissions could also be observed, as done in the past by SPICAV on Venus
Express (Chaufray and Bertaux 2015). Comet 2P/Encke should reach its perihelion on 26
June 2020 and 23 Oct 2023. While the geometry is not favourable for observations near
June 2020 by PHEBUS, observations will be possible at the end of November 2023. At the
end of November 2023, slightly after the perihelion, the comet-sun distance will be ~ 1 AU
and the comet–BepiColombo distance equal to 0.6 AU (see Fig. 7). This will allow optimal
observation of the comet tail, while the MDM measurements 1 month earlier and 3 months
later will be able to detect the dust amount originated from the comet.
<figure><img src="image_11.png" /><figcaption>**Fig. 7** Geometry of possible PHEBUS observation of comet 2P/Encke in November 2023. Blue cross is for BepiColombo at that date, and black dashed line indicates the direction of PHEBUS observation crossing the trajectory of cope 2P/Encke (in red)</figcaption></figure>

### 2.1.3 Other Science Goals: MGNS

The Curiosity Gamma-ray and Neutron Spectrometer (MGNS) is a scientific instrument developed to study the elementary composition of Mercury's sub-surface by measurements of neutron and gamma-ray emission of the planet. MGNS measures neutron fluxes in a wide energy range from thermal energy up to 10 MeV and gamma-rays in the energy range of 300 keV up to 10 MeV with energy resolution of 5% at 662 keV and of 2% at 8 Mev (Mitro-fanov et al. 2021, this journal) thanks to the innovative crystal of CeBr (Kozyrev et al. 2016). During the cruise to Mercury, it is planned that the MGNS instrument will operate continuously (excluding solar electric propulsions periods) to perform measurements of neutrons and gamma-ray fluxes for achieving two main goals of investigations:

1. monitoring of the radiation background of the prompt spacecraft emission due to bombardment by energetic particles of Galactic Cosmic Rays. This data will be taken into account at the mapping phase of the mission on the orbit around Mercury. Detailed knowledge of the spacecraft background radiation during the cruise will help to derive the data for neutron and gamma-ray emission of the planet during the nominal phase at Mercury starting from 2026. In fact, many elements like Si, O, Al, C, Mg, and others – the abundance of which at the uppermost layer of the planet is studied – are also present in the material of the spacecraft (Evans et al. 2017). Figure 8 shows the spectrum of the gamma-ray background of the spacecraft measured during the initial part of the cruise. Indeed, the nuclear lines of Al, Mg, and O are well-pronounced in the spectrum, which are also expected to be detectable in the gamma-ray spectrum of the Mercury emission.

2. participation in the InterPlanetary Network (IPN) program for the localization of sources of Gamma-Ray Bursts in the sky (Hurley et al. 2013). The localization accuracy by the interplanetary triangulation technique is inversely proportional to the distance between the spacecraft that jointly detected a GRB. Before the launch of BepiColombo, the IPN network included a group of spacecraft in the near-to-Earth orbit (e.g., Konus-Wind, Fermi-GBM, INTEGRAL, Insight-HXMT) and the Mars Odyssey spacecraft on orbit around Mars. Now, MGNS provides an additional interplanetary location, potentially increasing the accuracy of GRBs localization. During the first 13 months of continuous operation, MGNS detected 24 GRBs. Since November 2019 the pre-set time resolution of 20 seconds for GRB profiles measurements was increased into 1 second, and downlink resources allocated. Since then, the corresponding GRB detection rate is increased to
<figure><img src="image_12.png" /><figcaption>**Fig. 9** Time profile of gamma-ray burst GRB191125A detected by MGNS gamma-ray spectrometer (1 s acquisition integrated counts)</figcaption></figure>

<figure><img src="image_13.png" /></figure>

about 2–3 per month. An example of MGNS GRB measured with time resolution of 1 sec is presented in Fig. 9.

Gamma-rays originating from solar flares are also detectable by MGNS. Solar flares are non-stationary and anisotropic processes, and the ability to observe them from different directions in the Solar System is crucial for further understanding of their development and propagation, as it has been demonstrated in the case of HEND instrument onboard Mars Odyssey (Livshits et al. 2017). Though MGNS has not detected any solar events during its first 13 months, as the solar activity will increase during the present Solar Cycle 25, many flares are expected to be detected.

The MGNS instrument will also perform special sessions of measurements during flybys of Earth, Venus, and Mercury with the objective to measure neutron and gamma-ray albedo of the upper atmosphere of Earth and Venus and of the surface of Mercury. Another objective is to test the computational model of the local background of the spacecraft using the data measured at different orbital phases of flyby trajectories. The low altitude flybys (such as the 551 km for 2nd Venus flyby and three 200 km flybys for Mercury) would be the
<figure><img src="image_14.png" /><figcaption>**Fig. 10** Case of BepiColombo-Earth opposition. Trajectories of different spacecraft (BepiColombo, Solar Orbiter, Akatsuki (Venus), Parker Solar Probe) for 18 March 2021, plus/minus three months. The Sun is at the center, the location of Earth is fixed on the +X axis (black asterisk)</figcaption></figure>

most useful for such tests, because the spacecraft would be shadowed from cosmic radiation at very different distances from the planet. Neutron and gamma-ray measurements during Earth flybys should also enable investigation of the interaction between solar wind and Earth environments as well as studies of spacecraft neutron and gamma-ray background upon its passage through the Earth's radiation belts.

## 2.2 Methods

### 2.2.1 Multi-spacecraft Coordination

The heliospheric investigations can greatly benefit of BepiColombo related observations during cruise, and the data achieved can improve the accuracy of the space weather predictions beyond the location of BepiColombo (e.g. at the Earth). As already pointed out, the journey of BepiColombo to Mercury and the operation around the Mercury orbit cover different portions of the solar cycle, starting from the late declining phase of cycle 24, and up to all the rising phase of cycle 25 (see Fig. 3). In the intervening period between the April 2020 Earth flyby and Bepi Colombo's Mercury orbit insertion, several opportunities for heliospheric multi-point observations are possible via coordinated activity between Bepi-Colombo and other active spacecraft, as well as with ground-based IPS observations. The opportunity to obtain measurements simultaneously at many different heliospheric locations inside 1 AU is unprecedented (see Figs. 10, 11, and 12). Potential coordinated observations that may be planned would involve the following observation platforms:

- NASA's Parker Solar Probe (PSP) and ESA's Solar Orbiter (SolO) (launched in August 2018 and in February 2020, respectively, and that will be operating in the inner heliosphere for 7 years);

- JAXA's Akatsuki spacecraft orbiting around Venus since December 2015 (mission extension to 2024 presently under review);

- ESA's JUICE spacecraft (launch date June 2022), whose flight path includes an orbit around the Sun in the inner heliosphere prior to embarking on a direct trajectory to Jupiter in November 2026;
**Fig. 11** Plasma flow studies. Trajectories of different spacecraft (BepiColombo, Solar Orbiter, Akatsuki (Venus), Parker Solar Probe) for 10 August 2021, plus/minus three months. The Sun is in the center, the location of the Earth is fixed on the +X axis (black asterisk), the arrow shows the direction of plasma flow

<figure><img src="image_15.png" /></figure>

**Fig. 12** A magnetic flux rope constellation of BepiColombo, Parker Solar Probe and Solar Orbiter on 18 March 2023. The Sun is in the center, the location of the Earth is fixed on the +X axis (black asterisk). The black curved line shows the calculated Parker spiral of solar wind plasma with a velocity of 400 km/s

<figure><img src="image_16.png" /></figure>

- NASA's Polarimeter to Unify the Corona and Heliosphere (PUNCH) spacecraft, to be launched in March 2023;

- NASA's STEREO A, which is an Earth orbiting, heliophysics observatory;

- Other Earth-orbiting spacecraft located at the L1 point such as SOHO, DSCOVR, ACE and Wind.

Below are the main interesting geometries of multi-spacecraft constellations to be analyzed:

- radial alignment of two or more spacecraft;

<figure><img src="image_17.png" /></figure>
- occultation, when two spacecraft lie on the opposite sides of the Sun, thus the heliographic longitude they face differ in 180°;

- magnetic alignment, when the spacecraft lie on the same Parker spiral (where the curve of the computed Parker spiral depends on the assumed, modelled or measured solar wind velocity).

For radial alignment studies, a distinction has to be done between two cases: (a) the alignment of several spacecraft at the same time (when each spacecraft is located on the same heliospheric longitude at the same time and facing the same heliographic longitude); and (b) the alignment that takes into account the plasma travel time (with retardation): in this latter case, the different spacecraft will be on the same heliospheric longitude at different times, hence they will have a chance to observe the same plasma parcel (provided that the solar wind velocity used to define the plasma travel time is accurate).

In order to understand solar wind acceleration, we need simultaneous observations of the solar surface, the corona, and the interplanetary solar wind. In addition to the in situ observations, the magnetic field line at the solar photosphere can be derived from high spatial resolution observations (i.e. by using Hinode/SOT Tsuneta et al. 2008). The global structure of the solar wind density and velocity in the inner heliosphere (0.2–1.0 AU) can be derived by the study of ground-based interplanetary scintillation (IPS) observations (e.g., Iwai et al. 2019), combined with global MHD simulations such as SUSANOO (Shiota et al. 2014), IPS-ENLIL (Jackson et al. 2015), or EUHFORIA (Pomoell and Poedts 2018) (see next Sect. 2.2.2). The solar wind structures obtained with these models, are combined with the in situ observations at the corresponding trajectory, in order to distinguish the time and spatial variations in the in situ data, and to understand background solar wind propagation characteristics (Fujiki et al. 2003). The ratio between Fe and H, the so-called First Ionization Potential (FIP) bias is used to determine the location of solar wind source based on the abundance of the elements. It can be measured by the PICAM and MIPA ion sensors of the SERENA package (Orsini et al., this journal) onboard MPO. The FIP bias at the low solar corona can be observed by EUV spectrometer onboard Hinode/EIS (Culhane et al. 2007) by its campaign observation modes. Therefore, we can find the origin of solar wind at the BepiColombo location by comparing the full disc mosaic observation of the FIP bias by Hinode/EIS and in-situ measurements by SERENA.

Heliospheric transients studies will greatly benefit from multi-spacecraft perspective. As discussed in detail above (Sect. 2.2.1), magnetic field observations from several space probes inside 1 AU will be important to shed light on the propagation and evolution of ICMEs, their shocks and sheaths, SIRs and fast solar wind streams.

Magnetic field and plasma parameters observations will also be of help to characterize the evolution of plasma across the heliosphere and across different scales, thus advancing the investigating dynamical properties of solar wind evolution. Specifically, thanks to multi-spacecraft locations and observations (BepiColombo, SolO, PSP, ACE, Wind) several solar wind properties will be highlighted: the changes in the solar wind composition at different locations, how turbulent features evolve due to both solar sources and in situ processes, the role of thermal and suprathermal particle populations and their source mechanisms, and so on.

In addition to the local plasma and field measurements, STEREO A, PSP and SolO carry heliospheric imagers (Eyles et al. 2009; Vourlidas et al. 2016; Howard et al. 2019) that can be used to probe the density structures at the location of the inner heliospheric probes, including BepiColombo. Then, the ground-based IPS observations can be also used to derive the global density and velocity distributions of the inner heliosphere including CMEs (Iwai
et al. 2019). MGNS, if operated during solar flares, can deliver important information on the gamma and neutron spectrum related to solar flares (see also Sect. 2.1.3).

The study of Solar Energetic Particle events can also benefit a lot from multi-point obser-
vations: one of the main advantages is the greatly enhanced capability to separate transport
effects from the dynamics of the particle accelerator at the source, which is smeared to a
great extent when particles are observed from 1-AU locations only (Desai and Giacalone
2016). Both the number of observing spacecraft and their varying radial and longitudinal
separation will be an asset in disentangling the effects of source dynamics and transport.
The SIXS-P instrument measuring energetic particles, if operated during the cruise phase
and combined with energetic particle observations of PSP, SOLO, and spacecraft at 1 AU,
may help understanding of the transport and acceleration processes.

Figures 10, 11 and 12 show some of the above mentioned configurations for Bepi-
Colombo, Solar Orbiter, Akatsuki (orbiting around Venus) and Parker Solar Probe:

1) on 15–18 March 2021 BepiColombo and Akatsuki will be behind the solar corona (hid-
den from Earth) close to each other (Fig. 10). Radio waves transmitted from Bepi-
Colombo and Akatsuki can cross the solar corona almost simultaneously to allow multi-
point measurements along the same diagonal. Parker Solar Probe will be close on the
Eastern side and Solar Orbiter will be on the Western side of the two occulting space-
craft. In March 2021 solar activity will be low, therefore conditions in the corona can be
expected to be more or less stable for a couple of days. Having Parker Solar Probe and
Solar Orbiter at different longitudes makes it possible to investigate eventual fast streams
or CMEs as they sweep across the different spacecraft. On 15 March BepiColombo and
Akatsuki will be radially aligned, on the 16th they will be radially aligned with retarda-
tion, i.e. Akatsuki will observe the same plasma package that left BepiColombo on the
15th. On 18 March BepiColombo will be in opposition with Earth (Akatsuki very close)
making it ideal for radio science.

2) The period around 10 August 2021 offers remarkable possibilities to coordinated mea-
surements, for BepiColombo flyby Venus, 1 day after Solar Orbiter flyby at Venus too.
The event is ideal for cross-calibration and eventual detection of magnetic holes around
Venus (see Fig. 11 and Sect. 5). Having Parker Solar Probe also radially aligned makes
it possible to have a 3-dimensional picture of solar wind (eventually fast streams, CMEs,
SEP) propagation. The plasma parcel observed by Parker Solar Probe on the 10th of Au-
gust at the indicated time reaches BepiColombo on the 14th of August at 2 a.m. if its
velocity is 400 km/s. The latitudinal difference between Parker Solar Probe and Bepi-
Colombo is small (~ 0.4°): thus, there is a good chance to observe exactly the same
plasma and investigate its propagation. Furthermore, the latitudinal distance between the
three spacecraft at Venus at 0.7 AU is ~ 1°. At a distance of 0.7 AU, this is enough to
observe plasma expansion. Before the “triple meeting”, Parker Solar Probe will be on
the opposite, eastern side of the Sun on the 7th of August 2021. There is a chance for
solar corona analysis with radio observations using the Parker Solar Probe and all three
probes along the same radial.

3) One of the several opportunities for coordinated measurements of magnetically aligned spacecraft will occur on 18 March 2023, when Parker Solar Probe, BepiColombo and Solar Orbiter will be very close to the same Parker field line (Fig. 12). Calculations were made by assuming solar wind velocity to be 400 km/s, but a time window can be defined using different velocity values. If the probes are also close in latitude, they can study the changes of the solar wind source region. In the case shown, Solar Orbiter, Parker Solar Probe, and BepiColombo are wider apart in latitude: they will be located at 0, 6 and 12° from the Ecliptic. With this constellation one can study the 3-D variability of solar wind
plasma. This will also offer the unique opportunity to investigate the spatial evolution of solar wind turbulence properties as well as to provide a deeper understanding of the nature of the intermittency, continuously debated between a temporal and a spatial phenomenon. Moreover, this will be helpful for testing Taylor frozen-turbulence hypothesis and investigating the role of non-stationarities (Taylor 1938; Alberti et al. 2020; Chen et al. 2020).

The possible multi-spacecraft constellations are not limited to the above cited cases; several further special geometries (e.g. quadrature, widely spaced) can be considered where the different latitudinal position of the different spacecraft may offer a 3-dimensional investigation of solar wind properties. On the other hand, before proposing a coordinated observation campaign, a series of constraints have to be further checked, e.g. the operability of the instrument, its field of view, position, etc.

### 2.2.2 Support of Space Weather Modeling

Important support to multi-spacecraft measurements may come from the space weather simulation tools presently available in the scientific community.

In particular, the model predictions of solar wind propagation and interaction with the planetary environments may be used to support the interpretation of the measurements performed during the cruise (both far from the planetary environments, and close to planetary flybys). On the other hand, the BepiColombo measurements will provide precious data to verify and improve current space weather models and prediction tools.

Subsequent detailed analysis on the ground of the data provided by active MPO and Mio instrumentation will overall allow a database of snapshots of the solar wind environment at different distances from the planet to be built up and to, thereby, deduce how Venus and Mercury interacts over time with the extant solar wind at a range of distances from its surface, thus yielding important insights into how different kinds of solar wind stimulate different planetary responses.

Several state-of-the-art models/tools are available to provide space weather predictions at the BepiColombo spacecraft:

1. the Interplanetary Scintillation 3D-reconstruction Technique (IPS analyses) which provides precise tomographic 3-D reconstructions of the time-varying global heliosphere (Jackson et al. 2011, 2015). This methodology incorporates both the background solar wind and ICMEs and iterates to provide the best boundary values of velocity and density that are present globally, as well as measured in-situ over the viewed volume. Also, it extrapolates magnetic fields from the solar surface to this same inner boundary.

2. ENLIL (Odstrcil and Pizzo 1999; Odstrcil et al. 2005; Odstrcil 2003), is a time-dependent 3-D MHD model of the heliosphere which solves equations for plasma mass, momentum, energy density, and magnetic field using the Total-Variation-Diminishing Lax-Friedrichs (TVDLF) algorithm. The standard way ENLIL is operated using magnetic field measurements from the solar surface to provide a quasi-stationary solar wind model with velocity and density boundary input parameters derived from the Wang-Sheely-Arge/WSA Potential Field Source Surface (PFSS) model (Arge and Pizzo 2000). Injected into this background are “cone” CME inputs (e.g., Luhmann et al. 2010) from LASCO coronagraph data. ENLIL supports mass, and mass flux conservation, heating terms over solar distance, and non-radial transport of structures from its inner boundary, which is usually set at 0.1 AU.
An upgraded analysis system to be utilized at Venus during the flybys melds the first two systems described above together, thereby allowing the IPS data to iteratively update and fit ENLIL modelling so as to ultimately provide a rapid forecast of Coronal Mass Ejections and shocks, as well as of CIRs at inner heliospheric planets using the ENLIL 3-D MHD model as a kernel. This can provide all the accoutrements of the ENLIL system, and/or IPS tomography and imagery. In addition, the combined system updates this modelling as the solar wind flows outward from the ENLIL inner boundary. The combined system can trace the trajectories of interplanetary magnetic field lines in 3-D, thereby enabling simulations of the magnetic connections from locations on the Sun to BepiColombo and to other contemporaneously flying spacecraft (such as SOHO, ACE, and STEREO-A), as well as the Parker Solar Probe which monitors ambient energetic solar particles.

Other available space weather simulation tools could be put in use in a similar way to
provide a useful forecast of relevant space weather effects at BepiColombo position:

- The European heliospheric forecasting information asset (EUHFORIA – Pomoell and Poedts 2018) is a physics-based simulation model of the inner heliosphere driven by boundary conditions based on empirical models. This simulation tool has been specifically designed for space weather forecasting purposes. EUHFORIA consists of two components: a coronal model and a heliosphere model including coronal mass ejections. The coronal model reconstructs a large scale model of the coronal magnetic field and make use of empirical relations to determine the plasma state (solar wind speed and mass density). These quantities are then used as boundary conditions to drive a 3-D time-dependent magnetohydrodynamics model of the inner heliosphere. CMEs are injected into the ambient solar wind modeled by using the cone model, with their parameters obtained from fits to imaging observations. Upcoming improvement of EUHFORIA will take into account both the CME internal magnetic field and a time-evolving solar wind.

- In the last years, lightweight, fast semi-empirical or physics-based models have been used to implement an ensemble modeling of the CME and solar wind transient propagation. Ensemble modeling incorporates the intrinsic limitation of information due to measure errors and lack of knowledge in the form of probability distributions. In practice, instead of a single run to forecast an ICME propagation, a set of runs, driven with input parameters extracted from suitable distributions are used to retrieve a distribution of output parameters. Among those models, the P-DBM (Napoletano et al. 2018) and the DBEM (Dumbović et al. 2018) can run thousands of single simulations in seconds and thus explore thoroughly the parameter space. The model outputs are the most probable ICME travel time and velocity at a heliospheric position, as well as the associated prediction uncertainties.

# 3 Radio Science

## 3.1 MORE: The Mercury Orbiter Radio Science Experiment

Radio science experiment, as already cited, is the only originally planned science activity during the cruise of BepiColombo. Superior Solar Conjunctions (SSC) can be used to carry out tests of general relativity, similar to those previously performed with the Viking and Cassini missions (Reasenberg et al. 1979; Bertotti et al. 2003). In the Viking experiment the spacetime curvature generated by the mass of the Sun (controlled by the post-Newtonian parameter $\gamma$, equal to 1 in General Relativity) was derived from measurements
of the time delay of radio signals sent by a ground station and coherently returned to Earth by means of an onboard transponder. The observable quantity in the Cassini determination was the frequency shift of the carrier (proportional to the spacecraft range rate). The use of a multi-frequency link at X and Ka band (7.2–8.4 GHz and 32.5–34.0 GHz) allowed a nearly complete cancellation of the noise due to the interplanetary and coronal plasma and an improvement by a factor of 50 of the Viking results (Bertotti et al. 1993, 2003).

The Viking experiment was limited by the lower frequency of the radio link, at 2.1–2.3 GHz. Now the cruise tests to be performed by MORE on BepiColombo will combine for the first time range and range rate observables, and a plasma noise cancellation system based on the use of multiple frequencies. The Mercury Orbiter Radio Science Experiment (MORE) is based on the use of two onboard transponders: a dedicated radio science Ka band transponder (KaT), supporting a coherent link at Ka band (both uplink and downlink), and a TT&C deep space transponder (DST) supporting an X band uplink (7.2 GHz) and a dual frequency, coherent X and Ka band downlink at 8.4 and 32.5 GHz, respectively (for a complete description of MORE, see Iess et al., this journal).

MORE has two additional science goals for the cruise phase:

- the advanced radio system, complemented by data from the ISA accelerometer, can potentially provide significant improvements in the determination of the spacecraft trajectory.

- a comparison between the standard navigation system and the augmented system available for BepiColombo's geodesy and relativity experiments (especially the Ka/Ka link) can pave the way for its adoption in spacecraft operations.

Quantifying these improvements in spacecraft navigation is a primary goal of the MORE team. The spacecraft was tracked also during the Earth flyby (April 10th, 2020), and collected data that will be used to possibly better model the spacecraft trajectory orbit and improve the orbit determination codes.

In addition, MORE may contribute to studies on solar corona heating and on acceleration of the solar wind that are important topics of solar and heliospheric physics. The key processes are propagation and dissipation of the magnetic energy between the solar surface and the outer corona. MORE will produce a wealth of plasma calibration data during the SSC used for testing relativistic gravity. These data include the uplink and downlink Total Electron Content (TEC) and their variation with time, as shown by Bertotti et al. (1993, 2003). Uplink and downlink TEC and its time derivative are inferred from a linear combination of, respectively, range and range rate measurements in the three radio links (X/X, X/Ka and Ka/Ka). The ability to separate the uplink and downlink TEC along with its variations was never possible before and it is peculiar to BepiColombo. In the geometric optics limit, TEC is proportional to the integral of the refractive index along the line of sight, which, in turn, is proportional to the plasma density at microwave frequencies. (Magnetic corrections to the refractive index will likely not produce detectable effects after 2–3 solar radii.) The data could therefore be exploited not only to probe the density of the solar corona down to a few solar radii (a method often exploited in the past, see, for example, Miyamoto et al. 2014, with earlier references therein) but also for correlative analyses with the plasma instruments onboard BepiColombo and other spacecraft (such as ESA's Solar Orbiter and JAXA/Akatsuki), as well as observations from ground and near-Earth satellites. The separate determination of the uplink and downlink plasma lends itself to the space-time localization of large plasma events (such as CME, or CIRs) along the line of sight, thus complementing the information provided by imaging of the solar corona (Richie-Halford et al. 2009). In addition, the tracking data acquired during SSC may contain a contribution from the coronal magnetohydrodynamic waves and, in principle, provide also the solar wind velocity (via a time delay in the phase fluctuations at X and Ka band due to the differential bending of the two radio waves).
## 3.2 ISA: Measurements of Non-gravitational Accelerations

The measurements of MORE will be complemented by the Italian Spring Accelerometer (ISA) (Iafolla and Nozzoli 2001; Iafolla et al. 2007, 2010, 2016; Lucchesi and Iafolla 2006, and Santoli et al., this journal) by providing measurements of the non-gravitational perturbations (NGP) acting on the MPO spacecraft. These perturbations are due to surface forces, like direct solar radiation pressure, that carry away by a small amount the spacecraft trajectory from a purely gravitational (geodetic) one. ISA is a three-axis instrument, i.e. it provides (once its data are properly calibrated and reduced) the full three-dimensional vector representing the overall non-gravitational acceleration acting on the spacecraft. The instrument has been designed to operate on a wide signal frequency band ($3 \cdot 10^{-5}–10^{-1}$ Hz).

BepiColombo marks the first time that a high-sensitivity accelerometer – fully dedicated to scientific measurements – is embarked on a deep-space mission. Accelerometers with similar performance have been employed just on Earth geophysics missions, such as CHAMP (Reigber et al. 2006), GRACE (Tapley et al. 2013), GRACE-FO (Kornfeld et al. 2019), and GOCE (Drinkwater et al. 2006). It has to be noticed that, by construction, the accelerometer is able to sense both NGP due to external surface forces and internally generated signals (e.g. micro-vibrations). ISA data can be used as well to monitor the platform behaviour in terms of vibrations and rotations produced by antennas, mechanisms, solar panels, reaction wheels, etc. A non-exhaustive list of measurements that could be performed at selected times in the cruise includes:

- NGP acting on the MCS;

- planet-induced gravitational gradients during flybys;

- NGP during superior solar conjunctions;

- MCS accelerations along the thrust-no thrust transitions;

- density changes within the magnetosphere (boundary regions) during flybys;

- detection of micrometeoroid and dust impacts.

The measurement of non-gravitational effects in various cruise phases would be an interesting result by itself, besides being a direct verification of the instrument capabilities before the nominal in-orbit phase. Furthermore, the SSC will feature the full tracking capabilities of BepiColombo, opening an interesting calibration opportunity for ISA. Providing an independent measurement of the transition from thrust to no thrust orbit arcs could be useful in order to better assess the solar electric propulsion performance.

The possibility of measuring gravitational gradients on the spacecraft in the flyby phases deserves special attention. Due to design constraints, ISA has not been placed in the MPO center of mass; therefore a (small) gravitational gradient by the planet arises between the center of mass and the position of ISA proof masses. Hence, it is expected that during flyby with relatively low altitudes at closest approach, where the contribution is higher, the gravitational gradient of the planet will be part of the signals measured by the accelerometer. While this could have been an issue in the nominal orbital phase (expected to be solved by the use of the so-called Schulte vector), it will provide a potential calibration signal during the flybys. With the ever-increasing interest in gravitational gradiometers as instruments for the direct measurement of the gravitational field of Solar System bodies, where accelerometers are coupled together to detect the gravitational gradient, such a measurement can be considered as a direct test of ISA potentiality in being a basic element of a space gradiometer.

A measurement to be potentially performed by ISA is the identification and possible characterisation of the transition regions within the planet magnetospheres during flybys.
Crossing the different magnetospheric regions (like bow shock, magnetosheath, magnetopause, etc.) related to the interaction between solar wind and planetary magnetic field, means to travel through areas characterized by fluctuations in the density of particles and ions. Possible changes due to drag from charged particles (Milani et al. 1987) could be identified during flyby phases in the magnetosphere, specifically during the crossing of the magnetopause, i.e. the outer boundary of the magnetosphere. Such an opportunity is strongly dependent on the flyby geometry with respect to the magnetosphere structure and on the encountered density changes. Moreover, solar activity (solar wind, CME, 11-year cycle) deeply affects the magnetosphere shape and the particle densities.

To provide effectiveness to this approach, ISA measurements will be correlated with data from other switched-on instruments during the flyby and providing information on magnetic field variation, plasma measurements, particles energy/density/composition derived from other onboard instrumentation, such as MPO/MAG, SIXS-P, BERM, MIPA, PICAM, PHEBUS, and all Mio sensors. Collection of data from analogous instruments onboard other spacecraft monitoring the Earth environment during the flyby would be beneficial to the observations.

A further possibility for ISA science is to corroborate possible detection of micro-meteoroids and dust impacts on the spacecraft. This method has been already proved with MMS (Magnetospheric Multiscale Mission) spacecraft (Williams et al. 2016); in the present case, ISA observations should be carried out together with the Mio/MDM instrument, which is devoted to dust particles detection but suffers limited field of view (see Sect. 2.1.2).

Finally, the cruise phase is an additional opportunity to directly test and validate the instrument performance of ISA in a relatively quiet environment. Indeed, it is expected to be much quieter than the orbital phase for a number of reasons (i.e., greater distance from the Sun, reduced on-board activity, and very stable attitude). Measuring the background noise in selected cruise periods will be useful to fine tune the instrument error model and to characterize its long-term behavior.

# 4 Earth Flyby

The first BepiColombo flyby at the Earth was needed to deflect the spacecraft into the inner Solar System, and towards the orbit of Venus. The Earth flyby occurred on the 10th of April 2020 with a closest approach (CA) at 04:25 UTC, and an altitude of 12684 km.

The maximum apparent size was 37° for the Earth, and ~0.5° for the Moon. The geometry of the flyby and of the relative positions of MCS, Earth and Moon are depicted in Fig. 13 (upper panels), together with the average position of the Earth's plasma boundaries (bow shock and magnetopause) crossed by the spacecraft, with relative time in UTC; in the bottom panel the variation of altitude and of the Sun-MCS-Earth and Sun-Earth-MCS angles are also shown.

As discussed in Sect. 1, during interplanetary cruise and excluding electric propulsion phases, the default spacecraft attitude has the +*Y* axis directed towards the Sun. The need for angular momentum load minimization was relaxed for closest approach ±1 day, corresponding to the period of scientific interest. The general attitude constraints shown in Table 9 (in the Appendix) apply here: at a Sun distance of 1 AU, it is possible to offset the Sun direction in the spacecraft composite +*YZ* plan in the range between +47 and −9.1°. A roll phase around the Sun direction is also possible (360° rotation). According to attitude constraints and instrument requests, MCS operated in a quasi-inertial attitude with the Sun along +*Y* and with a phase of 250°.
<figure><img src="image_18.png" /><figcaption>**Fig. 13** Upper panels: BepiColombo trajectory on a Geocentric Solar Ecliptic System (GSE) X–Y plane (to the left) and GSE X–Y plane (to the right). Bow shock and magnetopause average positions are drawn in green and cyan respectively. In red the trajectory. Moon position is also shown in yellow. Lower panel: Linear and angular distances (Sun-MCS-Earth and Sun-Earth-MCS) of BepiColombo to the Earth’s center in the 48 hours around closest approach</figcaption></figure>
Moreover, power flux density constraints implied a switch from HGA to MGA to LGA respectively 7, 2, and 1 days before (and after) closest approach. 34 minutes of solar eclipse occurred at CA+36 minutes (the maximum duration allowed, as driven by the capacity of the MTM battery). During eclipses, the link between MTM and MPO was turned off, so each module is provided with its own power battery to be charged to 100%.

## 4.1 Scientific Objectives

The Earth flyby was the first opportunity to operate several instruments of BepiColombo at the same time. The instrument operations were mainly for calibration purposes, to observe the well known environments of the Earth, and the Earth-Moon system. The geometry of the flyby (see Fig. 13) offered a very close approach to Earth surface (< 2 $R_E$), while the Moon was farther away from the spacecraft in the opposite direction (> 300.000 km). The Moon is often used as a calibration target due to its well-known flux. Those calibration observations may sometimes result in major scientific discoveries, as happened in the past for the temporal and spatial variability of Moon surface hydration rate as observed by Deep Impact (Sunshine et al. 2009), or the detection of adsorbed water and hydroxyl on the Moon by Cassini (Clark 2009), and Chandrayaan-1 (Pieters et al. 2009).

Both MERTIS and PHEBUS observed the Moon, before and after the CA respectively. Both instruments required pointing. MERTIS obtained the first hyperspectral data ever of the Moon in the thermal infrared from space with a spatial scale of 500 km. Apart from using this type of observation for validating the calibration of the instrument, MERTIS measurements provided a new science dataset of the Moon surface composition and a baseline for comparison with the future Mercury datasets (i.e. the two celestial bodies, Moon and Mercury, being very similar in terms of surface appearance and exosphere). PHEBUS also observed the Moon, in two time-slots after the CA with the purpose of checking the absolute calibration of both EUV and FUV channels and characterize the pixel to pixel sensitivity of the two detectors (e.g. flatfield correction). These data will be compared to previous observations of the Moon in the UV range.

Regarding the measurements of the Earth environment, interesting data were acquired related to the crossing of the different regions of the planetary magnetic field (see Fig. 14).

The MPO-MAG started high rate observations (128 Hz) 33 hours before CA which allowed observing the interplanetary magnetic field (IMF) while outside the Earth's bow shock and the terrestrial magnetosphere after crossing the magnetopause. In total, the spacecraft spent 14.5 hours within the magnetospheric cavity. These measurements are of great importance to check the absolute sensor orientation. Simultaneous measurements of the magnetometer MGF onboard Mio, as well of other magnetometers onboard other spacecraft orbiting around the Earth (Cluster, Themis, MMS) have been acquired. In particular, MGF used the data of the solar wind for instrument offset determination, and the Earth magnetic field for inter-sensor alignment corrections.

PICAM and MIPA sensors from the SERENA package acquired measurements consistent with the ones observed by MPO-MAG. In particular, different ion populations have been identified with different energy levels and counts when crossing the magnetic boundaries, the outer radiation belts, the plasmasheet, and the low latitude boundary layers. The two SIXS detectors measured the proton and electron profiles well before and after the closest approach, including the range CA±7 hours, the near-Earth solar wind, the Earth's foreshock region, and the magnetosheath; and used them to cross-calibrate the detector elements with each other. In addition, by simultaneously observing with the X-ray detection system, SIXS calibrated the X-ray sensor with the particle background.
<figure><img src="image_19.png" /><figcaption>**Fig. 14** Projection of the BepiColombo trajectory on GSE frame, X-Y plane. The blue area highlights the Earth's magnetosphere, and the violet area the magnetosheath. The trajectory of the BepiColombo (red line) crosses the magnetopause and bow shock (Shue et al. 1998; Jeřáb et al. 2005). The two mesh-circles around the Earth represent the average location of inner and outer radiation belts. The Moon is not shown because outside of the figure box</figcaption></figure>

Measurements by the accelerometer ISA attempted to detect the variation of the charged drag on the MCS due to the variation of particle density (see Sect. 3). Also instruments on-board Mio had a unique opportunity to monitor the Earth magnetosphere, plasmasphere, and radiation belts, with relevant observations for calibration and sampling purposes. During the Earth flyby period, the plasma particles by MPPE-MSA, MEA, MIA, HEP-ele, and ENA, the magnetic fields by MGF, and the magnetic waves by PWI have been observed.

In addition, the MGNS detectors were on, with the goal to measure neutron and gamma-ray albedos of the upper atmosphere of the Earth, as well as background signal during the passage through the radiation belts near the phase of closest approach (see Sect. 2.1.3). The data from Earth flyby will help to properly interpret MGNS data at Mercury, and the effects of Mercury's magnetosphere over the MGNS measurements.

Finally, the geometry of the flyby provided favorable conditions for exploring spacecraft orbit determination of ingress and egress arcs. In fact, current orbit determination algorithms do not treat ingress and egress arcs with sufficient precision to reconcile the predicted and measured delta-v: this difference is often mentioned in the literature as the “Earth flyby anomaly” (Turyshev et al. 2009).

The joint measurements of MORE and ISA may reduce the uncertainty on the measured delta-v and, help to better understand the origin of this discrepancy.
<table><caption>**Table 2** Venus flybys dates and their main characteristics. No eclipse and no occultation occur in both cases. Electric propulsion will not be used close to first flyby (then, classified as N/A in the table)</caption><thead><tr><th></th><th>Venus flyby 1</th><th>Venus flyby 2</th></tr></thead><tbody><tr><td>Date</td><td>15 October 2020</td><td>10 August 2021</td></tr><tr><td>Time</td><td>03:58:31 UTC</td><td>14:00 UTC</td></tr><tr><td>Closest Approach altitude (km)</td><td>10721.6</td><td>~ 550</td></tr><tr><td>Venus-Sun distance (AU)</td><td>0.72</td><td>0.72</td></tr><tr><td>Venus-Earth distance (AU)</td><td>1.17</td><td>1.25</td></tr><tr><td>Sun-Earth-MCS angle (deg)</td><td>38</td><td>25</td></tr><tr><td>Electric propulsion end before CA</td><td>N/A</td><td>CA-30d</td></tr><tr><td>Electric propulsion start after CA</td><td>N/A</td><td>CA+7d</td></tr></tbody></table>

# 5 Venus Flybys

After the Earth flyby, the cruise to Mercury includes two consecutive flybys at Venus to reduce the perihelion to nearly Mercury distance and place the spacecraft in the same orbital plane as Mercury. The flybys take place on October 15th, 2020 and August 10th, 2021, and their main characteristics are summarized in Table 2. Figure 15 shows the flyby's trajectories in the Venus Solar Orbital (VSO) frame where the +X-axis points towards the Sun (left), the +Y-axis against the Venus orbital velocity vector and the +Z-axis northward. The map in the background shows expected ion mass fluxes and green and cyan lines show average positions of the bow shock and ion composition boundary layer (after Martinecz et al. 2009). Figure 16 provides the flybys linear and angular distances of the MCS, during the ingress and egress phases.

The two flybys are quite different both in terms of altitude and trajectory. During the first flyby the MCS approaches Venus coming from the dayside evening flank, moving towards the tail of the induced magnetic field, and crossing both the bow shock and the ionopause, with the closest approach altitude being above 10000 km. Conversely, during the second flyby the MCS will approach Venus down to about 550 km above Venus surface, and will come from the opposite direction, from night to dayside, but still passing close to the evening side (see Fig. 15).

In particular, the lower minimum altitude of the second Venus flyby will allow the MCS to fly over the sunlit part of Venus, providing Venus visibility to instruments with a field of view on $-Y$. The general attitude constraints, at a Sun distance of 0.7 AU, allow an offset of the Sun direction in the range of $+27^\circ$ to $-9.1^\circ$ in the spacecraft composite $+YZ$ plane (see Table 9 in the Appendix). Present configurations of both Venus flybys are eclipse free, therefore no power constraints will be necessary in this case. In addition, no Earth occultation by Venus will occur in both cases though, depending on the selected attitude profile, possible blockage of the antennae coverage by the spacecraft may occur.

## 5.1 Scientific Objectives

At the beginning of the space travel era, Venus was one of the most visited planets by space probes. In the last two decades, only two missions were launched to orbit around Venus, though long lived and fully successful: the ESA mission Venus Express (2006–2014), and the JAXA-Akatsuki spacecraft (2015-active). Nevertheless, there are still many open issues related to this planet that require dedicated investigations, and regular monitoring. In addition, Venus is becoming increasingly important in the rapidly growing field of exoplanet science (Kane et al. 2019; Dong et al. 2020) because exo-Venus analogs may be common
<figure><img src="image_20.png" /><figcaption>**Fig. 15** The flyby trajectories of BepiColombo at Venus on October 15th, 2020 (Bepi I, red) and August 10th, 2021 (Bepi II, blue) projected onto the VSO X-Y (top) and X-Z (bottom) planes. Numbers plotted along the orbits give the expected universal time. Green and cyan lines show average positions of the bow shock and ion composition boundary after Martinecz et al. (2009) The color map in the background shows the expected ion mass fluxes in [amu/(cm²s)] from a 3D MHD model with spherical grids based on average solar wind and ionospheric parameters using the PLUTO code (Mignone et al. 2012)</figcaption></figure>

types of planets and the near-future exoplanet detection and testing characterization methods at Venus is important to help future detections.

Hence, while waiting for a devoted mission to Venus, the frequent flybys of other spacecraft at Venus, necessary for the optimization of the orbits, are good opportunities to take additional data of the Venus environment (in terms of atmospheric measurements of composition, temperature, and wind vertical profiles, as well as of neutral and ionized particles, magnetic field, dust particles, etc.).

### 5.1.1 Atmosphere

**Composition and Thermal Structure** Key facts about the atmosphere of Venus were consolidated after the first exploratory phase carried out by Venera and Pioneer programs.
<figure><img src="image_21.png" /><figcaption>**Fig. 16** First and second Venus flybys linear and angular distances of the spacecraft: to the planet (green), to Sun and Venus (blue), to Sun and Earth (orange)</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Fig. 17** Example of air temperature profile in the Venus atmosphere, as derived from radio occultation experiment onboard Venus Express (from Limaye et al. 2018, after Tellmann et al. 2009). The tropopause is the approximate location of the main cloud decks</figcaption></figure>

The thick atmosphere of Venus is $10^{-4}$ the mass of the planet, much higher than in the case of Earth and Mars. The surface temperature is 735 K, and the pressure is 92 bars. The most abundant atmospheric gas is CO$_2$ plus several percent N$_2$ and other gases (notably, SO$_2$ and CO) as minor constituents. The water abundance is very low (between 1 and 15 ppm at the cloud top, Cottini et al. 2015). The vertical structure of the Venus atmosphere is presented in Fig. 17: within the height range of 0–55 km, the temperature lapse rate is nearly adiabatic. Between approximately 30 and 50 km, the decrease of temperature along altitude is slower than expected by the adiabatic gradient (up to 4 K/km), resulting in the formation of a region very stable against convective overturning (Seiff et al. 1985).
The entire planet is covered by layers of aerosols at heights between 49–70 km. The occurrence of clouds allows us to track a number of dynamical phenomena. A main feature of atmospheric general circulation is the global wind super-rotation at intermediate latitudes, with velocities between 60–120 m/s (Hueso et al. 2015). Upwelling of absorbing materials from deeper cloud regions associated with a Kelvin-like trapped wave was proposed to justify the 90°-tilted “Y pattern” observed in UV images extending symmetrically on the two hemispheres up to 45° (Peralta et al. 2015). Despite thick clouds, a small but non-negligible part of the solar flux penetrates to the surface and heats the atmosphere due to the greenhouse effect. The thermal structure of Venus atmosphere above the 55 km level was extensively studied by VIRTIS (Visible and Infrared Thermal Imaging Spectrometer) onboard Venus Express (see Limaye et al. 2017, 2018, and references therein), with measurements on the night-side, preferentially over the southern polar regions. VIRTIS data confirmed and provided further details about the mesospheric thermal fields already known from the data of the Venera 15 Fourier Spectrometer (Zasova et al. 1999).

At night, the temperature field exhibits a high degree of symmetry between the two hemispheres, as expected from the very small axial tilt. At levels around 55 km, air temperature tends to decrease monotonically from the equator to the poles (from 290 to 240 K), as expected by intense absorption of UV radiation by clouds at the subsolar points. At the 65 km level, dynamical effects begin to dominate, with two symmetric “cold collars” showing a minimum temperature of 220 K at 65S and 65N. Above 70 km and at least up to 90 km, equatorial regions are colder than polar regions, with monotonically increase at fixed altitude. This behaviour is consistent with a global scale Hadley circulation, where polar heating is caused by adiabatic compression of descending air. Mapping by VIRTIS unambiguously reveals patterns in local-time/latitude air temperature averages on the night side of the planet, tentatively associated with wave activities. Day-to-day variability was also observed, being particularly evident in the cold-collar region. On the dayside of Venus a few hundreds of high vertical resolution temperature profiles were obtained from radio-occultation measurements of the VeRA experiment (Tellmann et al. 2009). The overall thermal structure seems consistent to the night hemisphere, but the resulting coverage remains sparse.

Air temperatures at higher altitudes (between 90 and 110 km) at all local times are still poorly constrained by available data, but suggest very strong variations. This variability is possibly related to the overall change in global circulation pattern (from Hadley to solar/antisolar) occurring at this altitude. Above 120 km, the dayside atmosphere is dominated by absorption of UV radiation by molecules, with locations of peaks values of temperature strictly following the subsolar point. The solar/antisolar temperature gradient around 120 km was found to be around 35 K.

**Aerosols and Clouds** The most striking feature of the Venusian atmosphere is the thick coverage of clouds that permanently preclude the visual observation of the surface, with a total opacity at visible and near-IR wavelengths that varies from 20 to 40 (for a review see Titov et al. 2018). Data from Pioneer probe nephelometer is the primary standard for inferring Venus' cloud characteristics and still the best available today (Knollenberg and Hunten 1980). Aerosol populations encountered by the nephelometer in its descent from the space towards the surface corresponded first to a population of sub-micron haze (so-called “mode 1” particles), with particle detection beginning at an altitude of ~100 km and lower. A second, larger component (modal radius of ~1.2 µm) is evident below 70 km. This component (“mode 2”) is the main constituent (in terms of mass) of the upper clouds of Venus, that is typically observed from space in visible and infrared. A local minimum in cloud opacity at about 57 km marks the transition at the middle/lower cloud deck, where larger particles
$(r_m = 3–5 \mu m, “mode 3”)$ are found. At 48 km of altitude a sharp decrease in aerosol opacity is observed, and below only much optically thinner diffuse haze and possibly discrete clouds of uncertain nature can be found (Grieger et al. 2003). In situ analysis (Hoffman et al. 1980), as well as remote IR (Zasova et al. 2007) and polarimetric measurements (Hansen and Hovnenier 1974) identified a liquid mixture of sulfuric acid and water as the main constituent of haze and upper clouds: consistently, the clearing observed at the altitude of 48 km occurs where the temperature in the Venus environment allows the sulphuric acid to evaporate. Other constituents must be present as well: UV observations show high contrasting details, demonstrating the existence of a still unidentified UV absorber, strongly variable in space and time (Markiewicz et al. 2007); chlorine, phosphorus, and iron in the deeper clouds form still unidentified components as inferred by the measurements by instruments onboard the VEGA balloon (Andreichikov 1987).

The large suite of Venus Express instruments operating from thermal IR (5 µm) to UV enabled a series of studies, mostly focused on upper clouds and hazes. Among the notable results, the decrease of upper cloud heights and scale heights poleward from about 50° in longitude (Ignatiev et al. 2009; Haus et al. 2014), the detection of detached haze layers and larger particles (1 µm above) above the cloud top at 75 km (Wilquet et al. 2012; Luginin et al. 2016), variations in size and refractive indices with local time and latitude (Shalygina et al. 2015; Wilson et al. 2008). More spare are the studies about deeper aerosols: Magurno et al. (2017) found an increase of “mode 3” particle size with latitude, that is qualitatively consistent with the increase of total opacity of deep clouds towards the pole reported by Cardesín Moinelo et al. (2008) from VIRTIS data.

**The BepiColombo Flybys** Near and during the time of flybys, Venus's day and night side atmospheres will be observed by MERTIS (radiometer and spectrometer) and PHEBUS spectrometer that will make complementary measurements, and will work together with the Akatsuki spacecraft presently orbiting around Venus, and with the Earth-based telescopes (see Sect. 5.2). MERTIS (Hiesinger et al. 2020, this journal) will allow sensing the atmosphere from the top to an altitude of about 60 km above the surface, just within the upper cloud deck. It will also be able to detect mid-infrared thermal emissions from the cloud top level at daytime, and from below the cloud top at nighttime. In particular, the infrared spectrometer TIR in the range 7-40 µm allows sounding the middle atmosphere of Venus in the altitude range 55-100 km and its cloud layer (specifically CO₂, SO₂, and H₂SO₄ have absorption bands in the 7-14 µm range covered by the spectral channel TIS of MERTIS). Up to now, only a very limited set of observations of Venus in the thermal infrared spectral range is available.

PHEBUS spectrometer (Quémerais et al. 2020), designed to observe the weak emissions of the Hermean exosphere, is not optimized to observe the bright dayside atmosphere of Venus with the FUV channel (145-315 nm), that may damage the detector. In addition, due to the power limit, only one detector (EUV at 55-155 nm, or FUV) can be switched on at a time during the flybys (but they can be used both, alternatively). However, both NUV channels (at 404 and 422 nm) can be switched on at the same time. PHEBUS plans to perform two different observation modes during the Venus flyby closest approach periods: (a) star occultation at CA in FUV range at nightside (in both flybys), and (b) EUV airglow observations at high altitudes (during the second flyby). Stellar occultation is done in inertial mode and could be performed at EUV or at FUV if pointing at the nightside, to avoid scatter light from the Venusian atmosphere (Montmessin et al. 2011). It will allow deriving CO₂ maps; while O₃ and SO₂ could be derived from the FUV channel.

MGNS neutron and gamma spectrometers (Mitrofanov et al., this journal) will also contribute to atmospheric investigations by detecting elemental composition in terms of C, O,
N, and H₂O, indicative for the complex processes of diffusion, vertical mixing and chemical processing occurring in the Venus atmosphere.

**First Venus Flyby** The spacecraft approached the planet from the solar direction, over the dayside. The closest approach occurred above the evening terminator of the planet (Fig. 15), and then the spacecraft moved away from the planet to the anti-solar direction, over the night side. During the approaching phase, MERTIS started to resolve Venus from a distance of 1.4 million km when Venus size was about the same as seen from Akatsuki spacecraft, and went on until 4 hours before CA. At the closest distance, MERTIS gained a spatial resolution of about 6 km at the cloud top, i.e. higher than all the currently existing datasets. The contemporary radiometric data obtained provided insights into temperature profiles both on the illuminated and the non-illuminated hemisphere. Between 7 and 4 hours before CA, MERTIS performed close-up dayside observations from late morning to late afternoon via noontime on Venus at low latitudes of the southern hemisphere. Due to the 100% cloud cover on the planet, most of thermal emissions are corresponding to temperatures at the upper cloud level atmosphere (60–70 km) (Zasova et al. 2007). The shorter wavelength edge of the 15-µm CO₂ band covered by TIS/MERTIS will help to retrieve temperature profiles from the cloud tops to slightly above using the inversion methods of radiative transfer (Rodgers 2000). Simultaneous cloud top structures and SO₂ gas abundances above the clouds can be estimated from the observed spectra over the 7–14 µm range of TIS/MERTIS. This spectral range observation from space will be the first measurement of this kind after Venera missions in 1980s (Zasova et al. 2007).

PHEBUS planned to perform star occultation in the night side of Venus, shortly after CA time, by observing perpendicular to the dawn and dusk limbs but unfortunately no night side observations were obtained. Stellar occultations have been performed in the past in a systematic way by SPICAV-UV onboard Venus Express at the nightside. Vertical profile of aerosols of the haze layer, CO₂ density and temperature have been derived from stellar occultations between 80–140 km (Bertaux et al. 2007; Piccialli et al. 2015). An ozone layer has also been detected between 90–105 km from detailed inspection of these stellar observations at the strongest Hartley band around 250 nm (Montmessin et al. 2011). However, such an identification has been obtained only sporadically. Night side distribution of sulfur dioxide (SO₂) content at the same altitude range was retrieved from SPICAV-UV star occultations as well (Belyaev et al. 2017). Systematic SPICAV measurements in SO₂ absorption bands at 190–220 nm and 270–290 nm provided a study of the annual variability for 2006–2014. The star occultation by PHEBUS, possessing a twice better spectral resolution, should have complemented such a long-term trend. Synergetic observations of the Venusian atmosphere from MERTIS on the dayside, and from PHEBUS on the nightside, may help also in the study of the atmosphere from the chemical and dynamical point of view, for they should be able to observe how the chemistry changes before terminator, soon after, and shortly before entering again into light. Additionally, PHEBUS performed dayside disk integrated albedo in FUV when faraway from Venus (in late August–early September 2020) with the goal to monitor albedo variability over time; the visible channel was used too to measure global average SO₂ content.

**Second Venus Flyby** For the second flyby at Venus, the spacecraft will approach the planet from its night side, and CA will occur again near the evening terminator, but slightly shifted towards the late afternoon side with respect to the first flyby (Fig. 15). Then the spacecraft will move away from the planet towards the morning direction. The minimum altitude will be now only about 550 km. The close-up observations by MERTIS will be
<figure><img src="image_23.png" /></figure>

done from early night side to evening terminator at low latitudes, with a similar strategy as 1st flyby. At the CA, spatial resolution at the cloud top will be about 360 m providing an unprecedented detailed look at the distribution of minor species and the cloud structure in the atmosphere of Venus. PHEBUS will observe the exosphere of Venus with the EUV channel. The Venusian exosphere is mainly composed of atomic H, He and O. Exospheric H and O are composed of a population at thermal equilibrium with the atmosphere and a suprathermal population. While the hot H population is observed systematically by UV spectrometers (Chaufray et al. 2012 and references therein), the hot O population has only been observed by Pioneer Venus Orbiter UV Spectrometer (Nagy et al. 1981) and by Venera 11 (Bertaux et al. 1981). Several emission lines could be observed, depending on the altitude of the line of sight of PHEBUS. The hydrogen Lyman-α brightness, should vary between ~2 to ~30 kR for a line of sight tangent altitude from 500 to 2000 km based on SPI-CAV/Venus Express observations (Chaufray et al. 2012; Chaufray and Bertaux 2015) and then it could be detected by PHEBUS at all altitude during the flyby; the Lyman-β (~300 times weaker) could be detected from 1000 km altitude. Extrapolated from disk observations by Cassini/UVIS during the Venus flyby (Gérard et al. 2011a), the He brightness should be of only a few Rayleighs at 1000 km, due to its fast vertical decrease (scale height ~70 km) and then difficult to detect; nonetheless, it will surely be worth to try. The O emission at 130.4 nm, based on Pioneer Venus Orbiter (PVO) observations (Nagy et al. 1981) should be of few hundreds Rayleighs at 1000 km. It could be a new first opportunity, since PVO, to detect the hot oxygen on Venus, not detected with SPICAV/Venus Express. Unfortunately, the geometry of this second flyby should prevent PHEBUS to detect other lines (like CO, and C) in the EUV range (Gérard et al. 2011b), but EUV observations faraway from Venus could be done to observe them. Periodicities in the EUV dayglow (O II 83.4 nm, O I 130.4 and 135.6 nm) in the Venusian upper atmosphere, possibly associated to periodicities in oxygen or photoelectrons fluxes in the thermosphere due to planetary waves or gravity waves, have been detected by EXCEED aboard Hisaki (Masunaga et al. 2015) and will be investigated using consecutive observations by PHEBUS. Stellar occultations could be performed before or after the closest approach as done during the first flyby.

### 5.1.2 Induced Magnetosphere and Plasma Environment

Unlike the Earth and Mercury, Venus has no intrinsic dipole magnetic field and thus is a perfect prototype of the solar wind interactions with unmagnetized planets (e.g. Futaana et al. 2017). Without the shielding of an intrinsic dipole magnetic field, the impinging solar wind plasma can easily approach the planet and interact directly with the upper atmosphere. Figure 19 illustrates the major features of the solar wind interaction with Venus' upper atmosphere as deduced from the solar cycle maximum measurements obtained by Pioneer Venus Orbiter (PVO). The interaction leads to various regions and boundaries, such as the bow shock, magnetosheath region, induced magnetosphere and ionosphere on the upstream side, and a complex magnetotail in the downstream region.

**Ionosphere** On the dayside of Venus, the upper atmosphere is ionized by solar extreme ultraviolet radiation and solar wind particle impacts, leading to the formation of an ionosphere (Schunk and Nagy 2009), observed for the first time in 1967 by Mariner 5 (Kliore et al. 1967; Fjeldbo and Eshleman 1969). CO₂⁺ is the primary ion created by photoionization, and it quickly transforms into O₂⁺ and O⁺ due to charge exchange reactions. As a result, O₂⁺ is the major ion species at low altitudes while O⁺ becomes most abundant at high altitudes. The peak density of O₂⁺ at low solar zenith angles is observed at 130–150 km
<figure><img src="image_24.png" /><figcaption>**Fig. 18** PVO measured altitude variations of magnetic field strength (thin line) and electron density (circles) in the Venusian ionosphere for low to moderate solar wind dynamical pressure on orbit 186 (left), moderate to high solar wind dynamical pressure on orbit 177 (middle), and very high solar wind dynamical pressure on orbit 176 (right) (from Elphic et al. 1981)</figcaption></figure>

altitude, and above 200 km for $O^+$. Mariner 5 observations showed, surprisingly, a highly variable but significant ionosphere on the nightside of Venus, which was later quantified by observations from the PVO in 1978–1992 (Miller et al. 1980). More interestingly, PVO also observed ionospheric density holes on the Venusian nightside, possibly caused by the pressure balance between observed radial magnetic fields and ionospheric thermal pressure (Brace et al. 1982). Though Venus has no significant intrinsic magnetic field, at times of high solar wind dynamic pressure a significant (~100 nT) horizontal magnetic field is induced in the ionosphere (examples of both situations are shown in Fig. 18, Elphic et al. 1981). The interplanetary magnetic field (IMF) penetrates into the ionosphere, especially during solar cycle minimum, and it can significantly modify the ionospheric density profile. When the thermal pressure of the ionosphere is higher than the solar wind pressure (especially during the solar cycle maximum), the magnetic field drops sharply to nearly zero across the ionopause (Fig. 18, left panel). As a consequence, also the position of the ionopause may vary significantly from about 250 to 450 km altitude: the Venus Express data showed that the ionosphere is highly magnetized (about 95% of the time), but only 15% during solar cycle maximum (Zhang et al. 2008).

**The Flow of Ions Across the Terminator** Since ions are expected to recombine with the ambient electrons during the night of Venus (lasting 58 Earth days), a steady supply of ions from the dayside must exist to maintain the nightside ionosphere. It was already suggested by McElroy and Strobel (1969) that the transterminator ion flow can deliver this supply. Such a significant quantity of ions is indeed transported from the dayside to the nightside across the terminator of Venus and was first observed by the Orbiter Retarding Potential Analyzer (ORPA) onboard PVO using only a few orbits of observations (Knudsen et al. 1980). These first observations showed a chaotic but generally tailward motion of $O^+$ ions with velocities of 1–8 km/s. Above 200 km altitude $O^+$ is the dominant ion species at the terminator, while $O^+$ and $H^+$ ions contribute less than 10% to the total ion density (Miller et al. 1980). The analysis of these data from 3.5 years of PVO observations showed that the
<figure><img src="image_25.png" /><figcaption>**Fig. 19** Schematic diagram of the solar wind interaction with Venus induced magnetosphere and upper atmosphere (from Russell et al. 2007)</figcaption></figure>

transterminator ion flow at Venus has a mean velocity of about 1–2 km/s at 300 km altitude and increases to 4 km/s at 600 km altitude (Knudsen et al. 1982; Knudsen and Miller 1992). Transient velocities up to 8 km/s have been observed too. The mean transterminator ion flux observed ($2 \cdot 10^9 \text{ cm}^{-2} \text{ s}^{-1}$) is sufficient to maintain the nightside ionosphere of Venus. The observations of the PVO ORPA experiment summarized in Miller and Whitten (1991) showed that the velocity of ions on the dayside increases with solar zenith angle and that there is a dawn-dusk asymmetry in the flow-pattern for which at the time no explanation was found. More recently, Fox and Kasprzak (2007) reported a dawn-dusk asymmetry in the thermosphere that may be related to the asymmetry in the ionospheric flow. Also, Fox (2008) showed that new models of the dayside ion production at Venus are in agreement with the transterminator flux at Venus but cannot explain the loss rates of more than $10^{25}$ ions/s as observed in the tail of Venus by the Venus Express ASPERA-4 instrument (Barabash et al. 2007a). It was suggested by Cravens et al. (1982) and Whitten et al. (1982) and shown by Miller and Knudsen (1987) that the transterminator flux at Venus is controlled by the terminator ionopause altitude and thus by the solar wind pressure. Comprehensive reviews of the interaction of the Venus atmosphere with the solar wind can be found in Dubinin et al. (2011) and Futaana et al. (2017).

**The Ionotail** As mentioned above, observations by the PVO ORPA instrument (Miller et al. 1980) showed that above 200 km of altitude $O^+$ is the dominant ion species at the terminator, while $O^+$ and $H^+$ ions contribute less than 10% to the total ion density. Different
<figure><img src="image_26.png" /><figcaption>**Fig. 20** Composition and energy distribution of the escaping plasma as observed by ASPERA-4 on VEX over 33 orbits during the period from 18 May to 30 December 2006: (a) position of the spatial regions (b), (c) and (d); (b) plasma sheet matrix for the +Y lobe; (c) a plasma sheet matrix for the -Y lobe and (d) a matrix for the IMB layer. Vertical axes are energy/charge (E/q), the horizontal axes are the position of the ion impact on the detector (a sensor mass ring, $R_m$), and the color codes the counts that were accumulated over all directions and averaged over occurrences. The red lines represent constant mass lines (from Barabash et al. 2007a)</figcaption></figure>

estimates based on observations by PVO at distances up to $10 R_V$ put the ion loss rate from Venus at about $10^{26}$ H$^+$ ions/s and $1.4 \cdot 10^{25}$ O$^+$ ions/s (Knuksen and Miller 1992). But ionospheric models (Fox 2008) suggest that depending on the strength of solar EUV also C$^+$, N$^+$, He$^+$, CO$^+$, CO$_2^+$ may be observed at altitudes above 400 km. In a rare crossing of the Venus ionotail at the L1 point ($750 R_V$ downstream) by the SOHO spacecraft in June 1996, Grünwaldt et al. (1997) identified C$^+$, N$^+$, and He$^+$, while the SOHO/CELIAS instrument was not sensitive enough to detect heavier ions.

The MIMI CHEMS instrument onboard Cassini spacecraft, when passed Venus at a minimum altitude of 603 km on 24 June 1999, observed pick-up of O$^+$ and C$^+$ ions at pericenter (Krimigis et al. 1999) but no quantitative observations have been published. The energetic particle spectrometer EPS onboard MESSENGER spacecraft, during its flyby at Venus on 5 June 2007, did not discover significant fluxes. The Fast Imaging Plasma Spectrometer (FIPS) discovered pick-up ions of possible Venusian origin, but unfortunately no further analysis of these observations has been published. Hence, the most recent observations of the loss of ions from Venus come from the ASPERA-4 instrument (Barabash et al. 2007b) on Venus Express (VEX, 2006–2014). Barabash et al. (2007a) report observations of O$^+$, H$^+$ and He$^+$ at flux ratios of $Q(\text{He}^+)/Q(\text{O}^+) = 0.07$ and $Q(\text{H}^+)/Q(\text{O}^+) = 1.9$ (see Fig. 20
<figure><img src="image_27.png" /><figcaption>**Fig. 21** Heavy ion flux intensity and direction around Venus in the VSO frame as observed by ASPERA-4 on VEX (years 2006–2012) in XZ plane (left) and XY plane (right). Color scale is the logarithm of the total average omnidirectional flux in each bin in [m⁻² s⁻¹]. White arrows indicate the average direction of the flow, and the arrow length corresponds to net flux (from Nordström et al. 2013)</figcaption></figure>

left side). The Q(H⁺)/Q(O⁺) was recently re-assessed by Persson et al. (2018), where they found a significant solar cycle dependence: in the solar minimum the ratio is 2.6, while the solar maximum it becomes 1.1. Nevertheless, the ratio of Q(H⁺)/Q(O⁺) of 2 on average is in agreement with the predictions of Lammer et al. (2006) that with any lost oxygen atom, two hydrogen atoms must be lost because they derive from water molecules. The total loss rates were later better quantified by Fedorov et al. (2011) to be Q(H⁺) = 7.1·10²⁴ /s, Q(O⁺) = 2.7·10²⁴ /s and Q(He⁺) = 7.9·10²² /s and updated by Nordström et al. (2013) to the values to Q(H⁺) = 1.4·10²⁵ /s and Q(O⁺) = 5.2·10²⁴ /s. Dubinin et al. (2013) derived the density and velocity of the O⁺ in Venus's tail, as shown in Fig. 21 (note that the density of the O⁺ remains well below 1 cm⁻³). More details of the Venus Express observations are reviewed in Futaana et al. (2017). Unfortunately the ASPERA-4 IMA instrument was not designed to separate minor species like C⁺ or N⁺ from O⁺. Also the separation between O⁺, O₂⁺, CO₂⁺ and heavier species is ambiguous, and therefore loss of ions heavier than O⁺ has not been reported so far. Recently ASPERA-4 measurements in the Venusian tail revealed an unexpected reverse flow, namely, a bulk plasma flow towards Venus. The amount of the reverse flow is significant, which influences the net outflow flux from the Venusian ionosphere to space (Kollmann et al. 2016; Persson et al. 2018). The reverse flow is frequently observed in the solar maximum period, and because of that the net outflow flux is reduced in the solar maximum to 2.0·10²⁴ cm⁻² s⁻¹ (compared to 2.9·10²⁴ cm⁻² s⁻¹ in the solar minimum). Interestingly, there are frequently cases where the proton and the oxygen have different bulk flow directions. The mechanism that produces the reverse flow is not yet identified, while several mechanisms can be considered, for example, the magnetic reconnection (Zhang et al. 2012) or manifestation of single particle motions.

**High-Energy Particle Shadowing by an Induced Magnetosphere** One question that can potentially be studied with BepiColombo at Venus is magnetic shadowing of high-
energy particles caused by the induced magnetosphere. This effect was studied at Mars by McKenna-Lawlor et al. (2012), who found shadowing to be effective up to MeV energies by comparing global hybrid simulations to data recorded by the SLED instrument aboard Phobos-2 (McKenna-Lawlor et al. 1990). Different parts of the induced magnetosphere are expected to contribute differently to particle shadowing. For energetic particles with speeds much larger than the bulk speeds of the flow, the most important agents providing shadowing are the enhanced magnetic field and the planet itself. We, therefore, expect the greatest effects from the magnetosheath, the magnetic pileup boundary and the planet's shadow, depending on the energy, particle species and conditions in the solar wind. Besides the energy range and the orbit, also the directional resolution and coverage of the instrument observing the flux will be important. BepiColombo/SIXS observes the electron (60–3000 keV) and proton (1–30 MeV) flux in five orthogonal directions, of which three remain unobstructed in cruise configuration. With a single spacecraft flying by the planet, the study has to be backed up by a global simulation of the planetary plasma environment to provide context and a prediction to be verified by observations. In addition to an instrument capable of measuring magnetic shadowing effects, high fluxes of energetic particles will be necessary for successful studies of magnetic shadowing. Solar activity during the flybys is not expected to be at maximum, but the second flyby will most likely occur already well into the rising phase of Solar Cycle 25, increasing the chances of observing high fluxes of solar energetic particles during the flyby.

**Solar Wind Charge Exchange Induced X-Rays** Interpretation of in situ space plasma measurements is complicated by mixed spatial and temporal variations and orbital coverage. Potential remote sensing techniques of planetary plasma environments addressing these issues include the Solar Wind Charge eXchange (SWCX) induced X-rays. The planetary SWCX X-rays are generated when multiple charged solar wind heavy ions are highly excited in charge transfer reactions with exospheric neutral species and transition to lower energy states emitting soft X-rays occurs (Dennerl 2010). The main source of X-rays at Venus is the fluorescent scattering of the solar X-rays in the atmosphere, but the SWCX induced X-rays can be distinguished from the total flux during favorable low solar X-ray conditions typical near solar minimum. The Venus SWCX X-rays are very confined to the limb making it challenging to distinguish them from the fluorescent scattering, compared to the more extended exosphere of Mars and around comets (Gunell et al. 2007). The X-Ray Spectrometer (XRS) observed a decrease in the X-ray count rate as Venus blocked the background sky during MESSENGER's Venus flyby on 5 June 2007, but no emission related to the planetary environment was seen, which was explained by the instrument's low-energy cutoff at ~1 keV (McNutt et al. 2008). SIXS-X may experience the same challenge due to its spectral range of 1–20 keV, but the detection probability is also influenced by the solar X-ray irradiance as well as the flyby geometry. Other ion sensors, as SERENA MIPA and PICAM, and plasma sensors from Mio may significantly contribute.

**First Venus Flyby** During the first Venus flyby, BepiColombo was embedded in the solar wind magnetoplasma upstream of the bow shock. Therefore, various aspects of solar wind turbulence can be studied depending on the available data. They include the scaling properties of the magnetic energy power spectra and their comparison to known features at 1 AU (see e.g. Bruno et al. 2009; Narita et al. 2011; Sahraoui et al. 2013), the nature of high frequency turbulence and possible wave-particle interactions processes (Sahraoui et al. 2009, 2010), and the nature of the waves (e.g., cyclotron wave) in the upstream pickup ions (or foreshock) region (Delva et al. 2015). From Venus Express data it is expected that when the
spacecraft moves towards the bow shock the magnetic field strength will increase as well as the wave activity, as shown by Guicking et al. (2010). Downstream of the bow shock, two different situations can be met depending on the orientation of the interplanetary magnetic field with respect to the shock normal: a quasi-perpendicular or a quasi-parallel bow shock. In the case of the quasi-perpendicular bow shock it is expected that downstream ion temperature anisotropy will develop, which can trigger various ion-scale instabilities such as mirror modes (Volwerk et al. 2008a,b, 2016; Bader et al. 2019). These waves have been observed beyond Venus's terminator and they can either be generated locally or are transported by the plasma flow in the magnetosheath (Sahraoui et al. 2003). In case of a quasi-parallel bow shock various ultra-low frequency waves can be generated in the foreshock region and propagate or be transported downstream (Fränz et al. 2017). The spacecraft crossed the bow shock around the time of closest approach and will enter the ionotail about 40 min later at a distance of $5 R_V$. This offers the opportunity to study the magnetotail dynamics. The magnetotail of Venus is prone to similar dynamics as the Earth's magnetotail, where magnetic reconnection plays an important role in magnetic flux and plasma particle transport (see e.g. Volwerk et al. 2009, 2010; Zhang et al. 2012; Huang et al. 2018). Therefore, it can be expected to see plasmons travelling tailward during this flyby as also shown by the MHD model with a heavy ion flux of about $10^8$ [amu/(cm$^2$s)]. These relatively high heavy ion fluxes should have allowed PICAM to measure the magnetospheric ions during the first flyby. Also, the MSA ion spectrometer on Mio may measure low count rates of heavy ions from its restricted field of view within the MOSIF shield. This may have allowed for the first time to measure the escape of minor ions like $C^+$ or $N^+$. The magnetic reconnection studies can be completed by identifying the nature of the high frequency waves (e.g., whistlers) that are responsible for the processes of wave-particle interaction near the reconnection site (ion and electron diffusion regions), as generally done in similar studies in the Earth's magnetosphere (Huang et al. 2016). Such identification can be achieved by measuring high frequency magnetic fluctuations, which is possible thanks to the Search-Coil (SC) magnetometers (Low-Frequency Search Coils LF-SC and Dual-Band Search Coil DB-SC) of the Mio/PWI (Plasma Wave Investigation) instrument (Kasaba et al. 2020). In particular, the DB-SC instrument can measure the Z-component of the magnetic field fluctuations up to 640 kHz, which should allow covering all the local characteristic frequencies of the plasma (e.g., plasma frequency, electron cyclotron frequency). This type of studies can actually be achieved across all the regions (shock, magnetosheath and ionotail) explored by BepiColombo. So far, no such high frequency wave measurements exist in the Venusian environment. MPO-MAG data are fundamental to such investigations.

The outcome of such studies will be further strengthened by any additional plasma data (ions and electrons) that would come from the Mio/MPPE instrument. In addition to observations of magnetic reconnection and heavy ion fluxes during the first Venus flyby, it is also interesting to further study the Venusian magnetotail current sheet flapping by the magnetometer, which has been previously observed by Venus Express (Rong et al. 2015). Moreover, the Venusian magnetotail topology (i.e., open, closed or draped field lines) can be investigated by using MPPE electron data and MAG data in conjunction with MHD model results as done at Mars (Luhmann et al. 2015).

## Second Venus Flyby

The CA closer distance and geometry of the second flyby will bring BepiColombo much closer to the planet and it will cross the bow shock far behind Venus's terminator, where the question about its nature, a shock- or a wave-like structure, can be addressed.

The situation for a quasi-parallel bow shock is quite different. Luhmann (1995), using PVO data studied selected orbits which probed the subsolar point of the magnetosheath and
found that the waves were mainly linearly polarized and transverse, with no evidence for regularly occurring slow mode like waves structures standing in the flow of the magnetosheath. The second flyby configuration is also very favorable to investigate the nature of plasma turbulence (incompressible vs. compressible) and the spatial evolution of its properties (e.g., scaling laws, intermittency, energy dissipation rate), from the shock to the terminator, as done in the Saturnian (e.g. Hadid et al. 2015), Jovian (Tao et al. 2015), Earth (Huang et al. 2017), Martian (Ruhunusiri et al. 2017) and Hermean systems (Huang et al. 2020). A recent similar study was carried out using Venus Express data in the Venusian magnetosheath (Xiao et al. 2018). However, the potential availability of data from the magnetometers and Mio/PWI (LF-SC and DB-SC), and plasma data from the Mio/MPPE instrument, should allow us to address this problem in a more complete way and over a broader range of MHD and kinetic scales/frequencies. The spacecraft will probably enter the ionosphere during closest approach at an altitude of ~550 km where BepiColombo is inside the magnetic pileup boundary (MPB). If either PICAM or MSA will have the spacecraft ram vector in their field of view, the instruments would encounter high fluxes of planetary heavy ions. MAG will also observe the horizontal magnetic field lines as IMF penetrates into relatively low altitudes. As shown in Fig. 21 stronger heavy ion fluxes can sometimes not only be observed in the tail region of the planet but also upstream, picked-up and accelerated by the induced electric field. This creates the so-called “ion plume”, discussed in, e.g., Dubinin et al. (2011), frequently observed also at Mars by NASA’s MAVEN mission (e.g., Brain et al. 2015; Curry et al. 2014).

The second BepiColombo Venus flyby which will scan the flanks of the Venus-Solar wind interaction region. As can be seen in the same figure, total heavy ion fluxes can reach up to $10^{11} \text{ cm}^{-2} \text{ s}^{-1} = 10^{7} \text{ cm}^{-2} \text{ s}^{-1}$ at distances below 2 Venus radii. In the magnetosheath there is a high chance to observe planetary ions picked-up by the solar wind flow. Thus, even if the spacecraft will not cross the ionotail during this flyby the chances to observe the planetary ion composition are high. From the point of view of the foreseen level of solar activity, the magnetic shadowing of energetic particles are more likely to be successful during the second flyby. The rapid crossing of the inner region of the induced magnetosphere improves the chances to relate the observed temporal variations with the magnetospheric structures rather than with variations in the ambient magnetic conditions, which can be assessed within minutes of traversing the inner parts of the system. It is likely that solar activity and X-ray irradiance is on average higher during the second flyby, favoring the detection of the SWCX induced X-rays in a geometry looking away from the Sun towards the limb of Venus.

In conclusion, the most significant plasma observations at Venus would be the detection of minor ion species since this was not possible with the limited mass resolution of previous ion spectrometers crossing the Venus plasma environment. Considering that the single-anode geometry factor of Mio/MSA instrument is about $10^{-5} \text{ cm}^2 \text{ sr}$, we expect up to 100 heavy ions per second to be observed by MSA on the second Venus flyby at pericenter distance of about $1.2 R_V$. On the first flyby pericenter distance of BepiColombo was ~16770 km (less than $2 R_V$). At this distance expected heavy ion fluxes are only about $10^{9.5} \text{ m}^{-2} \text{ s}^{-1} = 3.1 \cdot 10^5 \text{ cm}^{-2} \text{ s}^{-1}$ such that the MSA heavy ion count rate may drop to a few counts per second.

However, during the long crossing in the ionotail mass fluxes are higher and count rates could increase when the spacecraft-Z-axis would point more towards Venus. MSA field of view will be restricted to about 5 anodes looking outward from the MOSIF shell. Ideal observing direction will be parallel to the electric field induced by the solar wind, that is perpendicular and towards the Sun-Venus axis. This look direction will also not suffer from saturation by solar wind ions. The expected maximum velocities of pickup ions are twice the
solar wind velocity. For heavy ions this means an energy of more than 10 keV, but the dis-
tribution may vary over a broad energy range well separated from solar wind. The PICAM
instrument on MPO is optimized to observe ions with energies less than the solar wind en-
ergy but with higher fluxes; hence, PICAM is well suited to analyze the ion composition both
in the ionotail and during the ionosphere passage in the second flyby. With the high mass
resolution of MSA and PICAM, the instruments may deliver the first quantified observation
of minor ion escape from Venus as C+, N+, CO+ or CO2+. Another interesting observation
would be the plasma turbulence (at MHD and kinetic scales) and the identification of differ-
ent plasma wave modes in the different regions of the Venus plasma environment. This can
be done primarily using the data from the MPO magnetometer (MAG). Moreover, thanks to
the Mio/PWI AC magnetometers (LF-SC and DB-SC), it will be possible for the first time
to measure high frequency turbulence waves in the Venusian system (if their amplitudes are
large enough to exceed the noise level of the interferences under the limited configuration
during the cruise phase). These studies will be further consolidated by the plasma data from
MPPE and PICAM.

## 5.2 Coordinated Observations

The Venus atmosphere observation opportunities of BepiColombo (Sect. 5.1.1) will become
more comprehensive thanks to the coordinated observation plan with the JAXA's Venus or-
biter Akatsuki, and with Earth-based telescopes. These three independent platforms will pro-
vide different viewing geometries, so that day and night sides of Venus will be captured si-
multaneously. Previously reported phenomena, such as thermal tides (Kouyama et al. 2019)
and local time dependent SO₂ gas distributions (Encrenaz et al. 2019; Marcq et al. 2020),
will be better understood in terms of spatial and temporal variabilities.

BepiColombo's Venus atmosphere observations will be performed by the UV spectrom-
eter PHEBUS, and the mid-IR radiometer and spectrometer MERTIS. The coordinated ob-
servations at the overlapped UV wavelength will allow cross validation of measured UV
data (Lee et al. 2019), between PHEBUS and a UV Imager onboard Akatsuki (UVI, Ya-
mazaki et al. 2018). That at the overlapped mid-IR wavelength will allow a comparison
between close-up spectral features of the cloud top level atmosphere obtained by MERTIS
and global thermal structure acquired by the Longwave Infrared Camera onboard Akatsuki
(LIR, Taguchi et al. 2007). A similar wavelength will be used also from Earth-based co-
ordinated observations, allowing detection of specific targets in the Venus atmosphere. For
example, UV to visible range wind vectors derived from Akatsuki and ground-based obser-
vations can be compared with thermal winds derived from the temperature map of MERTIS
(Peralta et al. 2017; Machado et al. 2017); SO₂ abundance measurements can be supported
with UV and IR observations from Akatsuki and the ground-based observations too (Jes-
sup et al. 2015; Encrenaz et al. 2019). When short time scale features of Venus will be
captured during approach and flybys of BepiColombo, Akatsuki will continuously monitor
Venus, and ground-based telescope will expand the spatial and spectral coverages of the
observations. These planned comparisons will contribute to a fundamental understanding
of Venusian clouds, chemistry, dynamics, radiative energy balance (Sánchez-Lavega et al.
2017; Limaye et al. 2018; Titov et al. 2018; Marcq et al. 2020), and climate and decadal
scale albedo trends (Lee et al. 2019).

**First Venus Flyby** Akatsuki has been operated since December 2015 (Nakamura et al. 2016), and its first extended operation phase is currently scheduled until March 2021. The Akatsuki UVI and LIR cameras are planned to perform Venus observations in coordination with PHEBUS and MERTIS, respectively.
<figure><img src="image_28.png" /><figcaption>**Fig. 22** Relative locations of the Sun, Venus, BepiColombo (blue), and Akatsuki (black) on 28 Aug–2 Sep 2020. The +X is assigned to the Sun. The grey curve is the trajectory of BepiColombo from August to November 2020. Red dots are the locations of BepiColombo and Akatsuki on 28 August 2020. The green arrows indicate the direction towards the Earth</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>**Fig. 23** Venus as seen from BepiColombo (left), Akatsuki (center), and Earth (right) on 28 August 2020. The size of Venus is not scaled, but the apparent sizes in degree or arcsec are written in the bottom. Solar phase angle (α, degree) and the distance to the center of Venus are also indicated</figcaption></figure>

UVI has two filters centered at 283 and 365 nm (Yamazaki et al. 2018): it monitored the Venus dayside, when PHEBUS performed dayside faraway observation on 28 Aug–2 Sep (see Figs. 22 and 23, and Sect. 5.1.1). The dayside was captured at three different solar phase
<figure><img src="image_30.png" /><figcaption>**Fig. 24** Relative locations of the Sun, Venus, BepiColombo (blue), and Akatsuki (black) on 10–20 October 2020. The +X is assigned to the Sun. Red dots are the locations of BepiColombo and Akatsuki at the time of the 1st Venus flyby CA (03:58 UT 15 October 2020). Green arrow indicates the direction of the Earth</figcaption></figure>

angles, including the one from the Earth (Fig. 23) on 20 Aug–3 Sep: this fact allowed a wide range of solar phase angles to be observed at the same time (for the first time ever). The coordinated faraway disk-integrated albedo observations performed are important to study the long-term variations of Venus's UV albedo (Lee et al. 2019) and to validate the observed radiance. These coordinated observations will be also used to retrieve the absorption spectrum of the “unknown absorber” in the Venusian clouds in the EUV-FUV range, by mean of the simultaneous 365 nm imaging of UVI, which reveals the most pronounced absorption of the “unknown absorber” (Pérez-Hoyos et al. 2018; Marcq et al. 2020). In addition, these faraway broad wavelength observations will extend the spectral coverage of Venus analysis as an exoplanet (Lee et al. 2020), which is a valuable reference to characterize Venus-like cloudy exoplanets.

LIR's filter is centered at 10 µm (8–12 µm bandwidth) (Taguchi et al. 2007). Recent discovery of planetary scale waves on the cloud level induced by the surface mountains were discovered (Fukuhara et al. 2017). Clear local time dependence of the waves had been reported (Kouyama et al. 2017), and the interaction between the surface and atmosphere implies a possibility of slowing down of solid body rotation rate (Navarro et al. 2018). At the time of BepiColombo's CA, LIR was monitoring the entire globe near the apocenter of Akatsuki's orbit (Figs. 24 and 25). Therefore, global scale thermal structures at the cloud top across the evening terminator could be obtained by LIR. Contemporary MERTIS's close-up observations over the same location can provide insights of the planet's atmosphere from its spectral analysis.

During the first flyby closest approach, Earth-based facilities have imaged the Venus' morning side, which was partially overlapped with MERTIS observations (late morning to noon area, see Fig. 25), but not by Akatsuki. Therefore, Earth-based observations can extend local time coverage to understand instantaneous global thermal tide, in addition to its long-term average (Kouyama et al. 2019). Also, SO₂ gas abundances can be compared between the global measurements from the Earth and the close-up measurement from MERTIS (Encrenaz et al. 2019). While a cloud-tracking method using UVI images can provide wind vectors over the afternoon side, Earth-based imaging and Doppler shift measurements can provide those over the morning side (Sánchez-Lavega et al. 2016; Machado et al. 2017).
<figure><img src="image_31.png" /><figcaption>**Fig. 25** Venus as seen from BepiColombo (left), Akatsuki (center), and Earth (right) at the time of the 1st Venus flyby CA (03:58 UT 15 October 2020). The size of Venus is not scaled, but the apparent sizes in degree or arcsec are written at the bottom. Solar phase angle (α, degree) and the distance to the center of Venus are also indicated. The background Venus image is LIR's observed mean thermal tide, taken from Kouyama et al. (2019)</figcaption></figure>

Many professional astronomers, and even amateur Venus observers, contributed to this effort (Garate-Lopez et al. 2020).

**Second Venus Flyby** If the Akatsuki mission extension will be determined, the coordinated observations with BepiColombo will be performed again during the 2nd flyby. In August 2021 flyby, BepiColombo and Earth-based observations, on the other hand, will observe the same area of Venus. This is a condition favourable to compare different observational target over the same area simultaneously, allowing comparison for different altitudes within the atmosphere using different wavelengths.

In addition to coordinated observations with Akatsuki, the ESA Solar Orbiter spacecraft will cross Venus during the same days, and have its own flyby at Venus on August 9th, 1 day before BepiColombo.

Solar Orbiter is an ESA spacecraft dedicated to the study of the Sun and of the inner heliosphere by means of remote and in situ instrumentation. Solar Orbiter in situ instruments are planned to be operative during Venus flybys, and will contribute to coordinated observations by investigating the variability of magnetic field and ionized environment, by means of measurements with the following instrumentation: Energetic Particle Detector (EPD), Radio and Plasma Waves (RPW), Solar Wind plasma Analyser (SWA) and the magnetometer (MAG). These instruments will coordinate with the BepiColombo sensors measuring electromagnetic field, plasma waves, and ions and electron populations at different energies (see Sect. 5.1.2 for details). Hence, during the second flyby, not only the atmosphere may be studied with more details thanks to the much lower altitude of BepiColombo at CA, but also the Venus' ionized environment and the planetary interactions with the solar wind can be better investigated, thanks to the very close (in time) flyby of a spacecraft devoted to the study of the inner heliosphere.

# 6 Mercury Flybys

Before entering into orbit around Mercury, BepiColombo will perform six flybys at Mercury that will allow progressive modifications to the interplanetary trajectory. The flybys will take place in 2021, 2022, 2023, 2024 (twice), and at the beginning of 2025.
<table><caption>**Table 3** Mercury flybys dates and their main characteristics. No eclipse will occur during Flybys 4 and 5; no occultation during Flybys 1 and 5 (then, classified as N/A in the table)</caption><thead><tr><th></th><th colspan="6">Mercury flybys</th></tr><tr><th></th><th>#1</th><th>#2</th><th>#3</th><th>#4</th><th>#5</th><th>#6</th></tr></thead><tbody><tr><td>Date</td><td>1 Oct 2021</td><td>23 Jun 2022</td><td>20 Jun 2023</td><td>5 Sep 2024</td><td>2 Dec 2024</td><td>9 Jan 2025</td></tr><tr><td>Closest Approach [km]</td><td>200</td><td>200</td><td>200</td><td>200</td><td>40000</td><td>393</td></tr><tr><td>Electric propulsion end</td><td>CA-38d</td><td>CA-39d</td><td>CA-44d</td><td>CA-85d</td><td>CA-38d</td><td>CA-75d</td></tr><tr><td>Electric propulsion start</td><td>CA+7d</td><td>CA+7d</td><td>CA+51d</td><td>CA+7d</td><td>CA+45d</td><td>CA+8d</td></tr><tr><td>Sun distance [AU]</td><td>0.38</td><td>0.38</td><td>0.33</td><td>0.32</td><td>0.32</td><td>0.45</td></tr><tr><td>True Anomaly Angle [deg]</td><td>260</td><td>263</td><td>312</td><td>331</td><td>331</td><td>150</td></tr><tr><td>Eclipse duration [mn]</td><td>13</td><td>18</td><td>15</td><td>N/A</td><td>N/A</td><td>24</td></tr><tr><td>Occultation duration [mn]</td><td>N/A</td><td>12</td><td>15</td><td>20</td><td>N/A</td><td>9</td></tr><tr><td>Sun-Earth-MCS angle [deg]</td><td>15</td><td>22</td><td>13</td><td>18</td><td>9</td><td>18</td></tr><tr><td>Earth distance [AU]</td><td>0.69</td><td>0.96</td><td>1.22</td><td>0.94</td><td>0.70</td><td>1.27</td></tr></tbody></table>

Table 3 provides a summary of the Mercury flybys main characteristics and Fig. 26 the spacecraft angular distances from the Sun, Mercury, and the Earth in the 48 hours around CA. Regarding attitude constraints, at a Sun distance < 0.7 AU, it is possible to offset the Sun direction in the spacecraft composite +YZ plan, in the range of +14° to -9.1° (see Table 9 in the Appendix). Eclipses will occur during Flyby 1, 2, 4 and 6, and the closest approaches will take place in the middle of eclipse and last 13, 18, 15 and 24 minutes, respectively. In eclipse, the link between MTM and MPO is cut; this is the reason why each module is provided on its own battery. The MPO battery will be charged to 100% for the Mercury flybys with eclipses. Earth occultations will occur in the Flyby 2, 3, 4, and 6 and will last 12, 15, 20, and 9 minutes, respectively.

Figure 27 (upper panel) shows the flybys geometries in the Mercury Solar Magnetic (MSM) frame, $X_{MSM}-\rho_X$ plane where $\rho_X$ is:

$$ \rho_X = \sqrt{Y_{MSM}^2 + Z_{MSM}^2} \qquad (1) $$

All close-in flybys (except for Flyby 4) have pericenters at the nightside of the planet. Flybys 1, 2 and 3 are roughly equatorial. Flybys 4 and 6 are polar. Figure 27 (bottom right) shows the flyby positions as seen from the Sun (YZ plane, MSM frame). Figure 27 (bottom left) shows BepiColombo Mercury flyby trajectories in the Mercury Solar Orbital (MSO) frame, Z-ρZ plane, compared to the MESSENGER flybys (dashed lines) and during its orbital mission phase (grey area).

All Mercury flybys (perhaps with the exception of Flyby 5) will be close enough to the planet to allow to take a wide range of novel and relevant measurements of the Hermean magnetosphere, exosphere, surface and interior.
<figure><img src="image_32.png" /><figcaption>**Fig. 26** BepiColombo linear and angular distances during the six Mercury flybys: to the planet (green), to Sun and Mercury (blue), to Sun and Earth (orange)</figcaption></figure>

**6.1 Scientific Objectives**

Full scientific objectives of Mercury investigations will be addressed in the nominal phase of BepiColombo, starting in 2026, when the two spacecraft will be orbiting around the planet in its final configuration. Despite the limitations in the operative instruments and performances, the 6 flybys will surely help to provide first BepiColombo data of the Hermean environment already five years in advance. In the following, some highlights of achievable scientific objectives are described.

### 6.1.1 Magnetosphere and Plasma Environment

Mercury is a natural laboratory for the comparative study of magnetospheric dynamics (Slavin et al. 2019). Mercury's internal magnetic field (190 nT $R_M^3$) is much weaker than
<figure><img src="image_33.png" /><figcaption>**Fig. 27** Top: the flyby trajectories of BepiColombo at Mercury projected onto the MSM frame, $X-\rho(x)$ plane. Black line is both for the planet profile and the magnetopause average position (Korth et al. 2017). Bottom right: BepiColombo Mercury flybys as seen from the Sun ($Y-Z$ plane, MSM frame). Bottom left: BepiColombo Mercury flyby trajectories in the MSO frame, $Z-\rho(z)$ plane, compared to the MESSENGER flybys (dashed lines). The light grey area indicates the satellite positions of MESSENGER in its orbital mission phase. The two dashes on the planetary profile mark the area studied by Winslow et al. (2014) using proton reflectometry data</figcaption></figure>

the Earth's (30400 nT $R_M^3$) (Anderson et al. 2012). Despite the weakness of the internal magnetic field of Mercury, the interaction with the solar wind leads to the formation of the bow shock and magnetopause boundaries as at the Earth. At Mercury heliocentric distance the solar wind density is 10 times higher than at the Earth, and the IMF intensity is about 10 times higher, while its velocity remains the same. As a consequence, the solar wind dynamic pressure at Mercury is on the average ~10 times larger. The high solar wind dynamic pressure, IMF and radiation environments at Mercury may mimic what exoplanets experience in close-in habitable zones of M dwarfs (e.g., Dong et al. 2018a; Airapetian et al. 2020), thus the investigation of solar wind interaction with Mercury may offer fresh insights concerning the exoplanetary habitability too. Meanwhile, the extreme solar wind conditions result in substantial differences between the Hermean and the Terrestrial magnetospheres interactions with the solar wind. The Hermean shock is much weaker, as the solar wind features lower Mach number and plasma $\beta$ (the ratio between plasma pressure and magnetic pressure) at Mercury than at Earth. Consequently, the stream of ions back-
<figure><img src="image_34.png" /><figcaption>**Fig. 28** Sketch of Mercury magnetosphere in active conditions, as seen from MESSENGER during its second flyby on 6 October 2008, showing the possible magnetospheric phenomena that could be observed during BepiColombo flybys investigations (adapted from Slavin et al. 2009)</figcaption></figure>

streaming into the incoming solar wind is weaker and the foreshock volume is smaller and structures have less time to evolve and steepen before they are carried to the shock (Le et al. 2013). The foreshock is more often located upstream of the subsolar shock, though the Parker spiral IMF angle is typically small (about 20°) (e.g. Slavin and Holzer 1981). Low β conditions are enhanced due to the extended Plasma Depletion Layer (PDL) that is commonly present upstream of the Hermean magnetopause (Gershman et al. 2013) for almost all IMF orientations. PDL contributes more significantly to the high reconnection rate at Mercury in comparison to Earth, leading to rapid successions of flux transfer events (FTEs) at the magnetopause (Slavin et al. 2012b; Imber et al. 2014; Imber and Slavin 2017; Zhong et al. 2020b), magnetotail flux ropes and plasmons (Slavin et al. 2009, 2012a; DiBraccio et al. 2015; Smith et al. 2017; Zhong et al. 2019), and dipolarizations or reconnection fronts (Sundberg et al. 2012; Zhong et al. 2018), as observed during the second Mercury flyby by MESSENGER (see Fig. 28).

**Mercury Flybys Investigations** Despite the recent MESSENGER observations, a full understanding of the reconnection and particles acceleration mechanisms at Mercury is still not achieved. In this sense, BepiColombo double spacecraft mission will greatly improve our understanding of the dynamics of the small magnetosphere and of plasma environment of Mercury, providing – for the first time – simultaneous multi-spacecraft observations during the nominal mission phase (Milillo et al. 2020, this journal). Nevertheless, each one of the six Mercury flybys, though with BepiColombo spacecraft in stacked configuration, will gather cross-sectional measurements of the environment, by crossing the foreshock, the bow shock, the magnetosheath, magnetopause boundaries and different magnetospheric regions.
Most of these flybys will cross the low-latitudes magnetospheric flanks that are less covered by the polar orbits of the two BepiColombo spacecraft (see Fig. 27).

Several scientific issues have been identified during BepiColombo Mercury flybys. The instruments that can be operated for the investigations of these regions are:

- MPO-MAG (fully operative) providing magnetic field at high time resolution,

- MGF and PWI providing magnetic field (at less time resolution since the booms are not deployed),

- the sensors of SERENA and of MPPE that will provide ion and electron full energy but partial angular distributions (due to their position in the stacked configuration),

- SIXS-P and BERM that will detect the high energy protons and electrons,

- and SIXS-X providing X-rays (see Table 1).

As shown in Fig. 27, during each of the Mercury flybys (except for Flyby 5), the spacecraft will cross the bowshock providing a good opportunity to measure the ion waves at the foreshock (with PWI experiment and the magnetometers MAG and MGF). Also the close-in flybys will be an opportunity to investigate the magnetopause reconnection and the dynamics and evolution of FTEs (with the measurements of magnetic field (MAG and MGF), electrons and ions (MPPE and SERENA). During the magnetospheric crossing SIXS-X and SIXS-P together with the BERM radiation monitor will provide information of possible solar events, like Solar Energetic Particles or Flares. These measurements, coupled to magnetic field and particle data, will help our understanding of Mercury's dynamic magnetosphere in relation to the Sun activity.

As already noticed, all close-in flybys (except for Flyby 4) have pericenters at the night-side of the planet, and Flybys 1, 2, and 3 are nearly equatorial. As a consequence, Flybys 1, 2, and 3 can be used to study dawn-to-dusk asymmetries in the tail, down-tail evolution of plasma sheet, lobes, and the nightside magnetotail radius from direct encounters (e.g. Slavin et al. 2010). These flybys are also uniquely designed for the study of the low-latitude boundaries, including the equatorial diamagnetic depressions expected on the dawn side (Müller et al. 2012), nightside reconnection (including x-line location (e.g. Slavin et al. 2009), re-connection dynamics (Zhong et al. 2018, 2020a), potentially substorm loading-unloading phases (e.g. Imber and Slavin 2017), and Kelvin-Helmholtz vortices at the magnetopause. Such measurements had been challenging from high-latitude MESSENGER orbits and will be challenging also for BepiColombo orbits during the nominal mission.

Flybys 4 and 6 are nearly polar. During Flyby 4, when approaching the Hermean geographical northern pole, quasi-steady "Region 1" field-aligned currents is expected (Anderson et al. 2014): these currents flow along magnetic field lines from the magnetopause region at dawn, down to the planet and emerge again on the dusk side flowing again towards the magnetopause. A big open issue is whether these currents are closed in between, either by the planet itself and/or by the magnetosphere. This depends on the available electrical conductivity, which has not been yet fully constrained.

Flyby 4 could also cross the northern cusp (e.g. He et al. 2017; Zhong et al. 2015) which provides an opening window for magnetosheath plasma to the magnetosphere.

Flyby 6 has the closest approach in the nightside high latitudes, and it could be used to study a magnetic Z-profile of the Hermean magnetotail to be compared with the average profile presently known from MESSENGER data (e.g. Rong et al. 2018). On the contrary, "Region 2" field-aligned current systems associated with substorm current wedges have not been observed so far with MESSENGER data (e.g. Glassmeier 2000; Poh et al. 2017).

Finally, the BepiColombo Flybys 4 and 6 could provide low-altitude measurements of magnetic field, electrons, and ions in the nightside, giving a chance to detect these signatures for the first time ever.
### 6.1.2 Intrinsic and Crustal Magnetic Field

The planet's weak internal global magnetic field (Ness et al. 1975), with an intensity of about only 1% of that of the Earth was discovered by the NASA mission Mariner 10 in the 1970s. The MESSENGER spacecraft, during its orbit around the planet in the years 2011–2015, characterized the magnetic field widely (Anderson et al. 2011), but not entirely. Due to MESSENGER's eccentric orbit, only data from the northern hemisphere close to the planet are available. The lack of information in overall the southern hemisphere makes the task of modeling the internal magnetic field very difficult and leads to correlations between the Gauss coefficients describing the magnetic field multipoles (see discussions in Iess et al. 2021 and in Heyner et al. 2021, this journal). The exact determination of the internal magnetic field is further complicated at Mercury, because the reconnection and induction change has very significant impacts on the amount of magnetic flux above the surface of Mercury (Jia et al. 2019; Dong et al. 2019). As a result Mercury's dipolar magnetic moment varies by up to at least 25% during solar wind transients (especially CMEs) (Slavin et al. 2019). Very low altitude magnetic field measurements obtainable, especially during the flybys, will allow a better determination of the dipole moment. The models of the Hermean core magnetic field have been proposed when studying the MESSENGER magnetometer data can be summarized as:

- a global core magnetic field model using the magnetic equator crossing positions as an initial constraint (Anderson et al. 2012);

- a regional model using a discrete regional method (Oliveira et al. 2015);

- a regional model using a continuous regional method (Thébault et al. 2018).

Two common features emerge in all models: the strong axisymmetry of the core field, and the location of the magnetic equator in the northern hemisphere. Only rough constraints could be established for the secular variation between Mariner 10 flybys and the MESSENGER orbital phase (Philpott et al. 2014) leaving the possibility that the field could not change within the 40 years that separated the 2 missions.

For latitudes between 35°N and 75°N, present day models of the Hermean crustal magnetic field were calculated by using the MESSENGER low altitude campaign data (Hood 2016; Hood et al. 2018). It was revealed that magnetic anomalies are distributed heterogeneously over the planet's surface, where the strongest anomaly (~14 nT at 40 km altitude) is related with the Caloris basin.

**Mercury Flybys Investigations** In 2021, six years after the end of MESSENGER's orbital mission, the first BepiColombo flybys to Mercury will allow us to look for hints of secular variation in the dynamo generated field (if they exist) and to confirm the asymmetry in the magnetic field. This would further constrain the type of dynamo operating in the Hermean core. Figure 27 (bottom left) displays the trajectories of the six BepiColombo flybys compared to the three MESSENGER flybys in the $Z-\rho(z)$ plane in the MSO frame.

The MESSENGER positions (Fig. 27, light grey area) from its orbital mission phase are shown for comparison. It can be noticed that, at the pericenters of Flybys 1 and 2, BepiColombo will move close to the planet in a region not previously explored by MESSENGER and will be able to sample the Hermean magnetic field in the southern hemisphere. Hence, the Hermean internal magnetic dynamo models established with MESSENGER magnetic field data (e.g. Wicht and Heyner 2014) can be critically tested. Also, these new measurements can be compared to the study by Winslow et al. (2014), who used the proton reflectometry method to assess the surface magnetic field at the southern hemisphere for the first time.
<figure><img src="image_35.png" /><figcaption>**Fig. 29** BepiColombo Mercury Flybys (here named M1 to M6) footprints on the Hermean surface. (top) BepiColombo sub-spacecraft altitude (coloured tracks), overlapping the topography (grey scale). The coloured map represents the crustal magnetic field strength measured by the MESSENGER magnetometer at 40 km altitude over the northern mid-latitudes (Hood et al. 2018). (bottom) BepiColombo sub-spacecraft illumination angle (coloured tracks) superimposed on the MDIS enhanced colour global mosaic. Illumination angles below 90° correspond to the dayside, and the closest approach of each flyby (below 1000 km) will mostly occur on the nightside. The maximum altitude represented here is 60000 km. The topography and enhanced colour global mosaics is gathered from USGS</figcaption></figure>

More challenging but still possible, will be the measurement of crustal magnetic fields during the flybys, in the case of stronger anomalies than the ones measured in the northern hemisphere by MESSENGER. Figure 29 displays the projected ground track of the BepiColombo flybys on top of a topographic surface map. The crustal magnetic field at 40 km altitude above the ground is depicted as a color overlay in a northern mid-latitude band. In particular, BepiColombo flybys 4 and 6 will pass over the crustal magnetic field regions mapped by MESSENGER; measurements during these flybys could be used to test the crustal magnetization analysis by Oliveira et al. (2019), because these flybys have appropriate pericenters altitudes. Magnetic field data from the remaining flybys could constrain crustal fields at regions not probed so far by MESSENGER, years before the beginning of the orbital science operations in 2026. Hence, detailed sampling of the planetary magnetic field by MPO-MAG, with sampling ≥ 1 Hz, should be planned to give crucial information to constrain both core field and crustal field models, and thereby shed more light on the dynamo processes operating in the deep interior.
### 6.1.3 Other Science Goals

**Exosphere** Also the study of the exosphere will benefit of the possible investigation during the flybys. The flybys geometries will allow detection and precise measurements of the exospheric components even at altitudes well below the final orbit of MPO during nominal phase (200 km to be compared to minimum MPO cruise altitude of ~450 km). Hence, it will be possible to derive vertical density profiles of different exospheric species already detected (Na, Ca, K, Mg, H, He, Fe, Al, Mn) but also to detect, for the first time, new species not yet discovered. The equatorial flybys (1, 2 and 3) can be used to investigate the dawn-dusk asymmetries, while polar ones (4 and 6) to monitor North-South asymmetries. The elongated trajectories of the flybys can allow scanning the elongated exospheric tail at distances up to few tens of Mercury radii and provide information on the acceleration processes, escape rates, and dynamics occurring in the tail. Integrated dayside and nightside exosphere emission and simultaneous imaging of different regions of the exosphere (as high latitude and equatorial, for example), can be additional measurements that will be worth performing during flyby.

All these measurements can be performed by PHEBUS (for details, see Quémerais et al. 2020). The environment sensors onboard Mio and MPO (primarily, MPPE and SERENA) may provide additional information on the ionized populations counterpart (at different energies) along the flybys trajectory paths, and investigate the link with the intrinsic magnetosphere and the interplanetary magnetic field to build a comprehensive view of the Mercury environment.

**Dust** The dust monitor MDM, if operating and with the reduced field of view in favourable configuration, can make measurements of the micrometeoroid population of the region crossed by the MCS during flybys, and investigate the spatial variability of the dust along the Mercury orbit (Kobayashi et al. 2020). Note that true anomaly angles of Mercury during flybys will be ranging from ~150° to ~330° (see Table 3), providing a wide range of heliocentric position to be investigated.

**Surface** Surface measurements will also benefit of the flyby campaign of observations. They will primarily come from MERTIS and MGNS.

MERTIS measurements will provide a first imaging of the Mercury global surface composition in the 7–14 µm spectral range, a region not covered by previous missions. In particular, feldspars can be spectrally identified and characterized, by means of several diagnostic spectral features (i.e. the Christiansen frequency, the Reststrahlen bands, and the transparency feature). In the thermal infrared range at wavelengths longer than 7 µm, spectral signatures in silicates result from characteristic fundamental Si–O vibrations. Therefore FeO- and TiO₂-free silicates (e.g., feldspars, Fe-free pyroxenes and Fe-free olivines), which are almost undetectable in the visible–NIR region, can be identified already during the flybys observation. Finally, MERTIS will implement tests for the proper instrumental set-up and calibration to be used during the nominal mission.

The MGNS measurements will make relevant tests for proper observation of the surface composition and mapping when in nominal phase around the planet. The major objectives for these measurements are for calibration purposes. The flux of the planetary emission is expected to increase with increasing solid angle of the planet (with the maximum value at pericenter), while the background of gamma-rays and neutrons of the spacecraft due to bombardment by Galactic Cosmic Rays will correspondingly decrease. By analyzing opposite variations of gamma-ray nuclear lines of different elements (like Si, Al, Mg, etc.), it will
allow discriminating between lines emitted by the planetary surface and lines emitted by the metallic structure of the spacecraft itself. The MGNS data of flybys will also allow making the direct comparison of the local background data (in terms of gamma rays and neutrons) for two quite different configurations of the spacecraft – the stacked configuration MCS and MPO – while they are at the same positions above the planet. In particular, also the spacecraft crossing of the Hermean magnetosphere boundaries with its high density of ionized particles may contribute to the counts of the MGNS detectors; hence, it should be known in advance before nominal mapping of Mercury surface to remove it (but it will be important information for the magnetospheric studies). Finally, the first proper data of the Hermean surface composition may also be obtained during the flybys (Mitrofanov et al. 2021, this journal).

## 6.2 Coordinated Observations

As for the case of Venus, coordinated Earth-based observations of Mercury during flybys (and then during the nominal mission) will be a valid support for in-situ measurements for exospheric and magnetospheric science. Many discoveries and studies of the exosphere of Mercury during the last three decades were done by means of Earth-based observations. They have also proven that the exosphere is very responsive to the surrounding environment, and in particular to the solar activity, even directly or through the mediation of the magnetosphere.

Coordinated observations of the exosphere are possible from the ground when Mercury has ample solar elongation angles, typically > 20°. Such efforts can be optimally timed utilizing solar telescopes to measure bright sodium emission at 5890 and 5896 Å occurring during the day from the exosphere. Sodium exosphere (discovered from the Earth in 1985, Potter and Morgan 1985) shows distinctive features, still not fully understood. In particular, the sodium exospheric shows typical 2-peaks intensity at high latitudes, and a highly variable dynamics over time scales ranging from minutes (Massetti et al. 2017) up to the planetary orbit time (Cassidy et al. 2015).

Historically, the solar telescope McMath-Pierce performed a long series of Na exospheric observations starting from the '80s, evidencing the two peak intensity pattern for the first time, as well as dawn-dusk asymmetry.

More recently, the solar telescope THEMIS has proven to be a particularly valuable facility for performing daily observations of the exosphere above the disk of Mercury for many hours per day, allowing studies on morphology and variability, as connected to the Interplanetary Magnetic Field, as well as on Solar activity and intrinsic magnetosphere (e.g. Leblanc et al. 2008, 2009; Mangano et al. 2013, 2015; Massetti et al. 2017; Orsini et al. 2018; Milillo et al. 2020). In addition, THEMIS works perfectly also for elongations down to 10° and it has recently installed an adaptive optics system that will greatly improve the spatial resolution of the Na exospheric maps produced. First Mercury observations are scheduled in early 2021. On the contrary, the McMath-Pierce Solar telescope has been now decommissioned (e.g. Potter and Killen 2008; Potter and Morgan 1990). It's role for future Mercury observations will be taken by the new generation solar telescope DKIST at Hawaii (the first call for proposals opened in summer 2020). DKIST, with the ViSP spectrograph, offers a promising follow-on: this next-generation telescope will initially have a 0.5°-25° exclusion zone in solar separation angles until safe operations with off-axis sunlight are characterized.

Observations of the Na tail are also of great interest: they show distinct, time variable north-south asymmetries within a few planetary radii and offer a several hour history of Na atmospheric escape rates at greater distances. These measurements are best performed
<table><caption>**Table 4** Earth-based observation parameters (elongation, illuminated disk fraction, phase angle) at the days of the Mercury flybys, and best suited observatories for daylight and twilight observations within 1 hour of CA. Suitable telescopes are: in Hawaii (Mauna Kea) Keck in evening twilight, CFHT in twilight; in Hawaii (Haleakala): DKIST in daylight and AEOS in twilight; in New Mexico: Dunn Solar Telescope in Daylight, Apache Point 3.5 in Twilight; in Arizona: Perkins in evening twilight, Discovery Channel Telescope in twilight; at the Canaries: THEMIS, GREGOR, and VTT telescopes in Daylight</caption><thead><tr><th>Flyby date</th><th>TAA</th><th>Elongation</th><th>Illuminated disk</th><th>Daylight obs. site</th><th>Twilight obs. site</th></tr></thead><tbody><tr><td>2021-10-01</td><td>260°</td><td>16°</td><td>12%</td><td>Hawaii</td><td>Hawaii</td></tr><tr><td>2022-06-23</td><td>263°</td><td>21°</td><td>52%</td><td>Canaries, New Mexico</td><td>Arizona</td></tr><tr><td>2023-06-20</td><td>312°</td><td>13°</td><td>86%</td><td>Hawaii</td><td>—</td></tr><tr><td>2024-09-05</td><td>331°</td><td>18°</td><td>47%</td><td>Canaries</td><td>—</td></tr><tr><td>2024-12-02</td><td>331°</td><td>9°</td><td>5%</td><td>Hawaii</td><td>Arizona</td></tr><tr><td>2025-01-09</td><td>150°</td><td>18°</td><td>86%</td><td>—</td><td>—</td></tr></tbody></table>

within short ~30 min viewing windows in twilight, which severely limits concurrent timings for coordinated mission support. Na tail observations can be made with traditional spectrographs within a few tens of planetary radii, or coronagraphic imaging at greater distances (e.g. Baumgardner et al. 2008; Schmidt et al. 2010). Observations of Mercury's potassium tail have been informally reported at conference proceedings and exhibit similar spatial structure to the ~60 times brighter Na tail. All other species in the exosphere tail require twilight observations at high spectral resolution. These must be measured off the bright planetary disk which scatters tens of MR/Å in solar continuum. The HIRES instrument on the Keck 1 telescope can observe calcium 4227 Å emissions of a few kiloRayleighs (Bida et al. 2000; Killen et al. 2005). HIRES has also made tentative detections of iron and aluminum emissions of order tens of Rayleighs (Bida and Killen 2017). Comparison with MESSENGER UVVS data suggested that Ca⁺ and Mn measurements may be within reach of HIRES as both at comparable brightness levels in the far blue (Vervack et al. 2016). All the more recently discovered faint species remain poorly characterized, being only studied so far as their initial detection.

Table 4 lists the dates of Mercury flybys together with the main Mercury parameters (as seen from the Earth), and potential optimal sites of Earth-based observations. Unfortunately during flybys, MCS configuration does not allow observations with SERENA neutral sensors (ELENA and STROFIO) perfectly suited for exospheric investigation; on the contrary, PHEBUS spectrometer will be operative. In addition, ion sensors will be able to operate (SERENA/MIPA and PICAM; SIXS, MPPE/MIA, and MPPE/MSA) and may give important information. Also MPPE/ENA may give additional important information, though partially shielded, of the energetic neutral atoms around Mercury, signature of high energy processes occurring in the interactions between the solar wind ions and the exospheric neutrals. Hence, Earth-based observations of the neutral components of Na, for example, will be a fundamental complement to the measurements of in-situ ion component of the exosphere and magnetic field, for a first BepiColombo coordinated detection of the particle environment of Mercury.

According to Doppler shift observed at the Earth due to relative velocities of Mercury and the Earth, optimal observation of the Na exospheric emission could be performed in the Flybys 2 and 4. Flyby 2 (23 June 2022) is of particular interest. In fact, it is the one that will have wider coverage from the Earth, with daylight observations possible from the Canaries and New Mexico, and twilight ones from Arizona, with a 50% of the disk illuminated, and highest elongation (and dawn terminator in view) (see Table 4). In addition, at
<figure><img src="image_36.png" /><figcaption>**Fig. 30** Mercury disk during Flyby 2 on 23 June 2022 as observed from BepiColombo (left) and from Earth (right) at three different times: 15:00 (top), 15:22 (center, time of CA) and 16:00 UTC (bottom). At the bottom of each planetary disk, first number is the apparent size of Mercury (in degrees or arcsecs) and α is the solar phase angle (in degrees)</figcaption></figure>

the Flyby 2 time, MCS is expected to approach Mercury from the night side and cross over the same portion of the surface viewed from the Earth shortly after. This means that, after in-situ measurements of the tail portion of the exosphere, direct comparison of global view images of the exosphere above the disk from the Earth with local in-situ measurements from PHEBUS can be performed. Figure 30 shows the illuminated portions of the disk and the expected longitudes as imaged by BepiColombo (left) and Earth-based telescopes (right) at three different times (before, during and after CA). The occurrence of both daylight and twilight observing site will allow coverage of both the tail and the disk exosphere, along the whole flyby.
# 7 Summary

BepiColombo, the ESA/JAXA cornerstone mission to Mercury launched on 20 October 2018, is an ambitious project for the study of such a mysterious planet, still full of open questions (even after the recent MESSENGER mission). Though BepiColombo will begin its nominal mission in early 2026 with its two spacecraft finally orbiting in two different orbits around Mercury, it will be able to provide interesting and valuable measurements also during its long cruise (7.2 years) into the inner Solar System, even in its packed configuration that obstructs part of the scientific payload. The almost continuous measurements of MAG, ISA, MGNS and BERM onboard MPO from the launch up to now, as well as the overall data obtained during the recent flybys at the Earth and at Venus, clearly proved the scientific value as well as the feasibility of dedicated measurements also during cruise, and in planetary environments significantly different from the one for which the instrumentation has been designed.

Cruise science is divided into three main strands, related to the “deep space” cruise science (far from the planets), to the 2 Venus flybys, and to the 6 flybys of Mercury.

BepiColombo cruise science is very promising, thanks to the fact that a big part of the available instruments can provide good measurements also for interplanetary investigations.

In fact, BepiColombo is provided with two magnetometers, one of which (MPO-MAG) has the boom fully deployed and can work all the time allowed by ESA constraints on thrust arcs or other operations. In addition, electron and ion detectors (SERENA, SIXS, BERM, Mio/MPP), as well as neutron and gamma-ray spectrometers by MGNS (at different energy ranges) can be operative too. Finally, coordination with other spacecrafts and satellites presently operating in other regions of the inner Solar System may give additional vanage points of observations, and may provide multi-point analysis of different interplanetary medium features. In Table 5, potential science cases during the cruise phase are summarized.

Venus flyby investigations are the other key-point of the interplanetary cruise science, for Venus will be approached twice, on 15 October 2020, and on 10 August 2021, with two interesting flyby trajectories. In both cases, the evening terminator region will be approached at CA time, once by coming from the dayside and at high altitude (about 10700 km), and a second time by arriving from the nightside, at a much lower altitude (down to about 550 km). The two flyby configurations can allow valuable investigations of many different regions and phenomena of the interplanetary-planetary induced magnetosphere, plus of the Venus' atmospheric composition and thermal profiles. In the following, Tables 6, 7, and 8 show a summary of the potential science cases to be investigated by the BepiColombo instrumentation during the flybys at Venus (divided by scientific goal: atmosphere, magnetosphere and mixed/of support to the previous two).

Finally, Mercury itself will be the last (and primary) focus of the interplanetary investigations of BepiColombo. The main target of the BepiColombo mission will be approached for the first time already in October 2021, when the first of six flybys will occur to allow progressive modifications to the interplanetary cruise trajectory to reach the final orbit around Mercury in December 2025. Hence, the six flybys will allow the first in-situ measurements of the different instruments of BepiColombo four years in advance before its orbit insertion. In addition, the different geometry of the Mercury flybys trajectories (and the low altitude of CA of many of them down to 200 km) will allow investigations of regions that have never been investigated by previous missions, and first BepiColombo measurements of magnetosphere, exosphere, surface and interior will be addressed.
<table><caption>**Table 5** Schematic diagram of the cruise science cases as related to the instrument devoted for the investigation, and to possible multi-spacecraft measurements</caption><thead><tr><th>CRUISE science</th><th>What should be measured</th><th>Instrument</th><th>When? Necessary constellation?</th></tr></thead><tbody><tr><td>Small-scale solar wind turbulence: how they evolve with radial distance</td><td>B field, high-res (min – sec – or less)</td><td>MPO-MAG</td><td>radial alignment, behind shock waves, within planetary magnetic sheats</td></tr><tr><td>Small-scale solar wind structures: magnetic holes, what is their relation to mirror mode waves</td><td>B field, high-res (min – sec – or less)</td><td>MPO-MAG</td><td>behind shock waves, within planetary magnetic sheats</td></tr><tr><td>Large-scale solar wind structures: transient CMEs, ICMEs properties and evolution, interaction of multiple CMEs, CME+SIR interaction</td><td>B field, 5-min – hourly</td><td>MPO-MAG (particle detector of MMO-MPPE is NOT directly exposed to solar wind), SIXS-P</td><td>radial alignment</td></tr><tr><td>SEPs: acceleration</td><td>fluxes of energetic particles (energy spectrum, pitch angle distribution), magnetic field configuration</td><td>SIXS-P, MPO-MAG, BERM</td><td>radial alignment + widely separated longitudes, close to shock waves</td></tr><tr><td>SEPs: transport (cross-field diffusion?)</td><td>fluxes of energetic particles (energy spectrum, pitch angle distribution), magnetic field configuration</td><td>SIXS-P, MPO-MAG, BERM</td><td>alignment along Parker spiral + widely separated field lines</td></tr><tr><td>Large-scale solar wind structures: background solar wind, CIRs – radial evolution, backmapping</td><td>fluxes of energetic particles (energy spectrum, pitch angle distribution), B-field, 5-min – hourly</td><td>MPO-MAG, SIXS-P</td><td>radial alignment</td></tr><tr><td>CIRs – magnetic sector boundaries</td><td>fluxes of energetic particles (energy spectrum, pitch angle distribution), B-field, 5-min – hourly</td><td>MPO-MAG, SIXS-P</td><td>same Parker field line</td></tr><tr><td>CIRs – 3d evolution</td><td>fluxes of energetic particles (energy spectrum, pitch angle distribution), B-field, 5-min – hourly</td><td>MPO-MAG, SIXS-P</td><td>latitudinal separation</td></tr><tr><td>Comets – water production rates in the hydrogen coma (+OH) – Panstarrs</td><td>EUV (Lyman-alpha) + FUV</td><td>PHEBUS</td><td>IMEX will calculate when each target intersects PHEBUS FOV</td></tr><tr><td>Comets – particles in the cometary trails</td><td>UV, dust impact flux</td><td>PHEBUS, MMO-MDM</td><td>constantly</td></tr></tbody></table>

<figure><img src="image_37.png" /></figure>
<table><caption>**Table 5** (Continued)</caption><thead><tr><th>CRUISE science</th><th>What should be measured</th><th>Instrument</th><th>When? Necessary constellation?</th></tr></thead><tbody><tr><td>Dust – interplanetary dust distribution in the inner solar system</td><td>dust impact flux within the MDM FOV</td><td>MMO-MDM (surrounded by MOSIF, but detection possible through MOSIF aperture)</td><td>constantly</td></tr><tr><td>Solar corona studies – density fluctuations, velocity of SW</td><td>radio occultation measurements, TEC</td><td>MORE</td><td>BC + Akatsuki, occultation</td></tr><tr><td>General relativity tests</td><td>time delay of radio signals</td><td>MORE</td><td>superior conjunction</td></tr><tr><td>Galactic Cosmic Rays</td><td>MCS emission due to bombardment (local radiation background)</td><td>MGNS, SIXS-P</td><td>constantly</td></tr><tr><td>Gamma-ray bursts: localize the sources</td><td>flux density</td><td>MNGS</td><td>constantly</td></tr><tr><td>Gamma-ray of solar flares</td><td>flux density</td><td>MNGS</td><td>constantly</td></tr></tbody></table>

<figure><img src="image_38.png" /></figure>
<table><caption>**Table 6** Schematic diagram of the Venus atmospheric investigations planned by MERTIS and PHEBUS</caption><thead><tr><th>Instrument</th><th>Channel/sensor</th><th>Sensitivity</th><th>Scientific investigation</th><th>Mode</th><th>Measurements</th></tr></thead><tbody><tr><td rowspan="10">MERTIS</td><td>radiometer TIS</td><td>7–14 µm</td><td>Venus as an exoplanet – track phase curve and rotation rate</td><td>inertial mode</td><td>radiance</td></tr><tr><td>spectrometer TIR</td><td>7–40 µm</td><td>Size of Venus similar for BC and JAXA</td><td>inertial mode</td><td>cloud feature, CO<sub>2</sub> temperature profiles</td></tr><tr><td>radiometer TIS</td><td>7–14 µm</td><td> Akatsuki</td><td>inertial mode</td><td>SO<sub>2</sub> abundance, cloud feature, CO<sub>2</sub> temperature profiles</td></tr><tr><td>spectrometer TIR</td><td>7–40 µm</td><td>Tracking of SO<sub>2</sub> abundance</td><td>inertial mode</td><td>temperature profile (day and night)</td></tr><tr><td>spectrometer TIR</td><td>7–40 µm</td><td>Daytime: detect mid-infrared thermal emission from cloud top level; Nighttime: detect mid-infrared thermal emission from below cloud top</td><td>inertial mode</td><td>spatial variation of thermal structure and trace gases</td></tr><tr><td>radiometer TIS</td><td>7–14 µm</td><td>altitude range 55–100 km (middle altitude atmosphere) to detect absorption of CO<sub>2</sub>, SO<sub>2</sub>, H<sub>2</sub>SO<sub>4</sub> and other trace gases</td><td>inertial mode</td><td>temperature profile (day and night)</td></tr><tr><td>spectrometer TIR</td><td>7–40 µm</td><td>Daytime: detect mid-infrared thermal emission from cloud top level; Nighttime: detect mid-infrared thermal emission from below cloud top</td><td>inertial mode</td><td>temperature profile (day and night)</td></tr><tr><td>FUV channel</td><td>145–315 nm</td><td>star occultation at CA toward nightside (both flybys)</td><td>inertial mode</td><td>atmospheric CO<sub>2</sub>, O<sub>3</sub> and SO<sub>2</sub></td></tr><tr><td>EUV channel</td><td>55–155 nm</td><td>not operating</td><td>—</td><td>—</td></tr><tr><td>NUV channel</td><td>404 nm</td><td>star occultation at CA toward nightside (both flybys)</td><td>inertial mode</td><td>other trace gases from the atmosphere</td></tr><tr><td>NUV channel</td><td>422 nm</td><td>star occultation at CA toward nightside (both flybys)</td><td>inertial mode</td><td>other trace gases from the atmosphere</td></tr></tbody></table>

<figure><img src="image_39.png" /></figure>
<table><caption>**Table 7** Schematic diagram of the Venus investigations on magnetic environment planned by MAG, SERENA, SIXS and Mio sensors: MPPE, PWI, MGF</caption><thead><tr><th>Instrument</th><th>Channel/sensor</th><th>Sensitivity</th><th>Scientific investigation</th><th>Mode</th><th>Measurements</th></tr></thead><tbody><tr><td>MAG</td><td>IB/OB</td><td></td><td>-measurement of bow-shock and other plasma boundaries<br/>-measurement of draped dayside magnetic field<br/>-study of low frequency wave activity inside the Ion Composition Boundary<br/>-potential measurement of flu ropes<br/>-confirmation of the tail lobes with opposite field polarity</td><td>high resolution:<br/>128 Hz range:<br/>128 nT resolution:<br/>3.9 pT</td><td>time series of magnetic field vectors</td></tr><tr><td>SERENA</td><td>MIPA</td><td>15 eV-15 keV</td><td>solar wind-Venus interactions and Venus magnetotail: better study boundaries and provide higher resolution than Venus Express. This will be at the expense of mass resolution, focusing on plasma dynamics</td><td>high resolution between CA-1h and CA+3h (mode 5); then mode 8</td><td>ion flux and energy before, during, and after crossing the bow-shock and the magnetopause</td></tr><tr><td>PICAM</td><td></td><td>1 eV-3 keV</td><td>solar wind-Venus interaction and Venus magnetotail</td><td>no attitude;<br/>Mode 1 Mode 4.<br/>Mode 7 (single pulse and Hadamard)</td><td></td></tr><tr><td>SIXS</td><td>SIX-X</td><td>1-20 keV</td><td>Measure X-ray emission from Venus atmosphere, X-ray sky background and particle background in X-ray sensors</td><td>science mode</td><td>Nominal monitoring with change of detector a few times</td></tr><tr><td></td><td>SIX-P</td><td>electrons: 0.05-3 MeV<br/>protons: 1-30 MeV</td><td>Monitor energetic electrons and protons + SEP effects detection on atmospheric escape</td><td>science mode 2 (high res data downlink)</td><td>Nominal continuous monitoring of protons and electrons from all unobstructed directions</td></tr></tbody></table>

<figure><img src="image_40.png" /></figure>
<table><caption>Table 7 (Continued)</caption><thead><tr><th>Instrument</th><th>Channel/sensor</th><th>Sensitivity</th><th>Scientific investigation</th><th>Mode</th><th>Measurements</th></tr></thead><tbody><tr><td rowspan="7">MIO (all)</td><td>MPPE/ENA</td><td>10 eV-3.3 keV</td><td>ENA (Energetic Neutral Atoms) detection on the dayside due to processes of backscattering and/or sputtering of solar wind ions over the exobase</td><td>FOV directed to Venus</td><td>energy input/particle supply/loss mechanisms between the solar wind and atmosphere of non-magnetized planet (tailward flux of ENAs)</td></tr><tr><td>MPPE/HEP-ele</td><td>30-700 keV</td><td>High energy electron distribution of the shock origin around Venus; Shielding effect of the high energy electrons by Venus</td><td>N/A</td><td>Density and energy of high energy electrons</td></tr><tr><td>MPPE/MSA</td><td>1 eV-26 keV</td><td>C<sup>+</sup>/N<sup>+</sup> pick up ions</td><td>N/A</td><td>Density and energy of each ion component</td></tr><tr><td>MPPE/MLA</td><td>15 eV-29 keV</td><td>atmospheric pick up ions</td><td>N/A</td><td>Density and energy of ions</td></tr><tr><td>MPPE/MEA</td><td>3 eV-26 keV</td><td>electron acceleration signature</td><td>N/A</td><td>Density and energy (temperature) of low energy electrons</td></tr><tr><td>PWI/EWO</td><td>0.1 Hz-20 kHz (B-field)</td><td>Lighting EM pulse detection (below 500 or 1000 km altitude)</td><td>N/A</td><td>Perturbation of magnetic fields around Venus; Bow shock (position, structure); Field line draping; Low altitude field enhancement/decay; Tail structure, reconnection</td></tr><tr><td>PWI/SORBET</td><td>2.5-640 kHz (B-field)</td><td>DC at low frequency field: Bow shock (position, structure); Field line draping; Low altitude field enhancement/decay; Tail structure, reconnection.</td><td>N/A</td><td>L mode</td></tr><tr><td>MGF/MGF-1</td><td>DC-64 Hz</td><td>Upstream foreshock waves; Proton cyclotron waves; Magnetosheath turbulence</td><td>M mode (AC field); L mode (Rem data)</td><td></td></tr><tr><td>MGF/MGF-0</td><td>DC-64 Hz</td><td></td><td></td><td></td></tr></tbody></table>

<figure><img src="image_41.png" /></figure>
<table><caption>**Table 8** Schematic diagram of the Venus investigations planned by MGNS, ISA and BIERM to support both atmospheric and magnetospheric science</caption><thead><tr><th>Instrument</th><th>Channel/sensor</th><th>Sensitivity</th><th>Scientific investigation</th><th>Mode</th><th>Measurements</th></tr></thead><tbody><tr><td rowspan="5">MGNS</td><td>neutron spectrometer STN</td><td>10<sup>-3</sup> eV-1 keV</td><td rowspan="5"><strong>Atmosphere:</strong> elemental composition of the Venus atmosphere;<br/><strong>Magnetic Environment:</strong> interaction between the solar wind and Venus atmosphere</td><td rowspan="5">angle between the direction to the centre of Venus and the -X axis &lt; 20 degrees; background science configuration during cruise is OK</td><td rowspan="5">C, O, N, H<sub>2</sub>O</td></tr><tr><td>neutron spectrometer SETN</td><td>0.4 eV-1 keV</td></tr><tr><td>neutron spectrometer MD</td><td>0.4 eV-500 keV</td></tr><tr><td>neutron spectrometer SHEN</td><td>300 keV-10 MeV</td></tr><tr><td>neutron spectrometer GRS</td><td>100 keV-10 MeV</td></tr><tr><td rowspan="3">ISA</td><td rowspan="3"></td><td rowspan="3"></td><td><strong>Interior:</strong> gravity gradient<br/><strong>Magnetic Environment:</strong> eventual acceleration effect due to Bow-shock crossing (entrance)<br/>eventual acceleration effect due to Bow-shock crossing (exit)</td><td rowspan="3">nominal and thermalized</td><td rowspan="3"></td></tr><tr><td><strong>Magnetic Environment:</strong> radiation environment</td></tr><tr><td></td></tr><tr><td rowspan="3">BERM</td><td>Electron Channel</td><td>0.3-10 MeV</td><td rowspan="3"><strong>Magnetic Environment:</strong> radiation environment</td><td rowspan="3">nominal</td><td rowspan="3">energetic electrons, protons, and ions</td></tr><tr><td>Proton Channel</td><td>1-200 MeV</td></tr><tr><td>Heavy Ions Channel</td><td>1-50 MeV/mg/cm<sup>2</sup></td></tr></tbody></table>

<figure><img src="image_42.png" /></figure>
## 8 Conclusions

In conclusion, in the present paper we collected a wide range of possible science cases that the BepiColombo mission will approach in the present years during its cruise phase towards Mercury. Science cases include mainly four different environments: the Earth, the planet Venus and planet Mercury during the flybys, and the interplanetary environment during cruise (far from planets). The science goals will be addressed by mean of strong synergy among the onboard instrumentation, but also in coordination with other spacecraft presently operating in the inner Solar System, and with Earth-based telescopes observations.

**Acknowledgements** The authors would like to thank the two anonymous referees who kindly reviewed the first version of this manuscript and provided valuable suggestions and comments to improve its quality.

A special thanks to Prof. Dario Del Moro (University of Rome Tor Vergata) for fruitful suggestions and support on space weather modeling section, and on editing.

In addition:

A. Madár and M. Dósa were funded by the National Research, Development and Innovation Office – NKFIH, grant FK128548.

Y.J. Lee has received funding from EU Horizon 2020 MSCA-IF No. 841432.

J.S. Oliveira is funded by the ESA Research Fellowship program in Space Science.

C. Dong was supported by NASA under Proposal Grants 80NSSC19K0161 and 80NSSC20K1416.

C. Schmidt acknowledges support from NASA programs 80NSSC19K0790, 80NSSC18K0857 and 80NSSC21K0051.

The MGNS instrument team in Russia is funded by Ministry of Science and Higher Education of the Russian Federation, grant AAAA-A18-118012290370-6.

The IMEX Dust Streams in Space model was developed under ESA funding (contract 4000106316/12/NL/AF-IMEX).

**Funding Note** Open Access funding provided by Istituto Nazionale di Astrofisica within the CRUI-CARE Agreement.

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## Appendix

<table><caption>**Table 9** Sun offset allowed range in Y direction during cruise. SAA is the Sun
Avoidance Angle; SEP here is for Solar Electric Propulsion system; (*) when not operational</caption><thead><tr><td>SAA details</td><td>SAA</td><td>Constraints</td></tr></thead><tbody><tr><td>minimum value</td><td>-9.1°</td><td>MOSIF cut-off plane at -18°</td></tr><tr><td>maximum at 0.8 AU</td><td>+47°</td><td>SEP illumination up to 50° (*)</td></tr><tr><td>maximum between 0.7 and 0.8 AU</td><td>+27°</td><td>SEP illumination up to 30° (*)</td></tr><tr><td>maximum at &lt; 0.7 AU</td><td>+14°</td><td>SEP illumination up to 17° (*)</td></tr></tbody></table>