# Dr Nicola Baresi

## About

### Biography

Nicola Baresi started his Astrodynamics career with a MSc thesis on “Optimal Control of Formation Flying Satellites”. After graduating full marks in Physics from the University of Padova in 2011, he moved to Israel where he worked as a postgraduate researcher for the Distributed Space Systems Laboratory of the Technion, Israeli Institute of Technology. Starting from 2013, Nicola moved to the United States of America as a US-Italy Fulbright scholarship awardee and pursued his PhD studies on spacecraft formation flying and dynamical systems theory. He was later awarded with a PhD in Astrodynamics and Satellite Navigation Systems from the University of Colorado Boulder, as well as with an outstanding graduate research award from the department of aerospace engineering sciences of the same university. Following graduation, Nicola was employed at the Japanese Aerospace eXploration Agency (JAXA), working on small and large scale satellite missions to the Moon and Mars. He eventually joined the University of Surrey in 2019, first as a Surrey Research Fellow and now as a lecturer in Orbital Mechanics at Surrey Space Centre.

## Research

### Research interests

- Astrodynamics
- Dynamical Systems Theory
- Spacecraft Formation Flying
- Autonomous Guidance, Navigation & Control
- Trajectory Optimization
- Orbit Determination
- Asteroids & Comets

## Teaching

- EEE3039: SPACE DYNAMICS AND MISSIONS
- EEEM009: ADVANCED GUIDANCE, NAVIGATION AND CONTROL

## Publications

Quasi-periodic invariant tori are emerging as a powerful tool to populate the phase space and enable a better understanding of key astrodynamics problems. This chapter introduces modern numerical continuation techniques for generating two-dimensional invariant tori in the elliptical restricted three-body problem of the Earth–Moon system. As a test case and baseline trajectories for the upcoming Lunar Orbital Platform Gateway, the (2:9) and (1:4) synodic resonant near rectilinear halo orbits are hereby considered. Their dynamical substitutes are first calculated and later analyzed using two-dimensional torus maps that facilitate the creation of accurate initial guesses for real-ephemeris orbits. The transition process is demonstrated with a trajectory optimization procedure that successfully generates continuous ballistic arcs around the Moon for more than one year.

Successive convex programming is a promising technique for onboard applications thanks to its speed and guaranteed convergence. Hence it can be an enabler for future missions where spacecraft autonomy plays a key role. The definition of a good value of the trust region plays a vital role in the successful convergence of SCVX algorithms. This work presents an improved trust region algorithm based on a differential algebra technique that relies on the information given by the nonlinearities of the constraints and does not depend on the user for the initialization of the trust region.

Periodic orbits in the Restricted Three-Body Problem are widely adopted as nominal trajectories by di↵erent missions. To maintain periodic orbits in a three-body regime, a stationkeeping strategy based on a high-order Target Point Approach (TPA) is proposed, where fuel-optimal and error-robust TPA parameters are acquired from stochastic global optimization. Accurate TPA maneuvers are calculated in a high-order fashion enabled by Di↵erential Algebra techniques. Orbit determination epoch is selected using a sensitivity analysis based on the convergence radius of a stroboscopic map. Stochasticity is handled by incorporating Monte Carlo simulations in the process of optimization and the evaluation of high-order ODE expansions is employed to supplant the time-consuming numerical integration. Two specific types of periodic orbits, Near Rectilinear Halo Orbits and Quasi-Satellite Orbits, are investigated to demonstrate the validity and eciency of the strategy.

For deep-space missions that remain in the vicinity of a target body, solar eclipses might become a source of major concern due to thermal and power constraints. This paper presents a method for minimizing the eclipse duration of periodic trajectories in three-body systems using synodic resonant periodic orbits. The proposed methodology successfully captures the global eclipse structure of synodic resonant periodic orbits in terms of driving parameters such as elongation and phase angles. Two-dimensional Eclipse maps are introduced to identify optimal orbit insertion conditions that avoid or minimize eclipse intervals. The validity and applicability of the proposed method is tested in the full-ephemeris model (DE430) of the Earth–Moon and Mars–Phobos systems, respectively. As a result of these investigations, we propose new science minimum-eclipse trajectories that are under consideration for the upcoming JAXA missions EQUULEUS and MMX. •Eclipses have the potential to endanger the safety and operations of space missions.•We derive relation between eclipse of periodic orbits and mission design parameters.•Minimum eclipse orbits can be achieved through the proposed graphical tool.•The validity of developments has been investigated and applied for two JAXA missions.

We present a small-scale mission concept to characterize the permanently shadowed regions of the lunar south pole. MARAUDERS aims to measure in situ for the first time the presence, distribution, and state of volatiles in one permanently shaded crater at a greater resolution than existing orbital measurements using up to 12 deployed impactors. A total of 15 permanently shadowed regions have been characterized as potential landing sites candidates for the probes. The science principle is based on penetrometry, that has proven in the past to be an efficient technique to estimate regolith properties from acceleration profiles. We demonstrate this concept by numerically simulating the surface interaction between our probes and the lunar regolith, thereby demonstrating how deceleration profiles can elucidate information on key regolith properties and help discriminate between two ice-regolith end-members. The preliminary payload design indicates that a good baseline for the impactors would be a spherical shell of 30–40 mm in size and ~90 g in mass per impactor, including electronics and the communication system. This would sum up to an overall payload of ~1 kg contained in a volume of ~15.10−4 m3, which is in agreement with a small-scale payload. Preliminary landing trajectory design enabled the computing of a nominal deployment scenario (with constraint on altitude, ejection velocity and spin rate) that would provide dispersions of the probes from ~250 m down to ~20 m if deployed from orbit, and down to ~10 m if deployed from a carrier lander/rover. Both scenarios will be able to comply with the MARAUDERS’ objectives to assess: (1) the presence (2) the distribution and (3) the surface strength heterogeneity (that can be traced back to the state of volatiles through lab experiments) of water-ice volatiles in permanently shadowed regions at a resolution

Numerical methods have proven to be reliable tools in generating long-term bounded relative trajectories in nontrivial dynamical environments. Yet, none of the existing procedures have a systematic approach to come up with all of the bounded relative trajectories associated with a variety of satellites in Earth orbit. In this paper, such a systematic procedure is developed based on the assumption that any trajectory, except at critical inclinations, is either periodic or quasi periodic in the four-dimensional Routh reduced system describing the motion of a mass particle about an axisymmetric body. This allows key design parameters to be identified, such as the nodal period and the node drift per nodal period, which are later used to constrain a differential corrector scheme looking for families of quasi-periodic invariant tori that yield bounded relative motion in the Earth-centered inertial frame. Numerical simulations show that the computed solutions are good candidates for cluster flight missions at a large range of altitude and inclination values.

Quasi-periodic invariant tori are of great interest in astrodynamics because of their capability to further expand the design space of satellite missions. However, there is no general consent on what is the best methodology for computing these dynamical structures. This paper compares the performance of four different approaches available in the literature. The first two methods compute invariant tori of flows by solving a system of partial differential equations via either central differences or Fourier techniques. In contrast, the other two strategies calculate invariant curves of maps via shooting algorithms: one using surfaces of section, and one using a stroboscopic map. All of the numerical procedures are tested in the co-rotating frame of the Earth as well as in the planar circular restricted three-body problem. The results of our numerical simulations show which of the reviewed procedures should be preferred for future studies of astrodynamics systems.

Modern space missions are frequently targeted towards new and unexplored regions of space, such as the region between the Earth and the Moon, which is denoted as the Cislunar space (NASA, 2020; International Space Exploration Coordination Group – ISECG, 2018); binary asteroid systems (Rivkin et al., 2021); comets and other irregularly shaped celestial objects; satellites of other Solar system planets. In all of these mission scenarios, the spacecraft dynamics is governed by an intriguing, yet complex and chaotic dynamical environment that is driven by the presence of multiple and/or non-spherical massive bodies. The gravitational influence of these objects shall be addressed with methods and techniques that are different from the standard Keplerian tools available in the classic Two-Body Problem. In recent years, the space community has shown a renovated scientific and technological interest in mastering the multi-body non-Keplerian astrodynamics for practical applications. Immediately, new technological and engineering challenges emerged in order to cope with this uninvestigated portion of outer space. In particular, the Guidance, Navigation and Control (GNC) and the Propulsive subsystems developments have been strongly supportive of this endeavor.

Upcoming missions towards remote planetary moons will fly in chaotic dynamical environments that are significantly perturbed by the oblateness of the host planet. Such a dominant perturbation is often neglected when designing spacecraft trajectories in planetary moon systems. This paper introduces a new time-periodic set of equations of motion that is based on the analytical solution of the zonal equatorial problem and better describes the dynamical evolution of a spacecraft subject to the gravitational attraction of a moon and its oblate host planet. Such a system, hereby referred to as the Zonal Hill Problem, remains populated by resonant periodic orbits and families of two-dimensional quasi-periodic invariant tori that are calculated by means of numerical continuation procedures. The resulting periodic and quasi-periodic trajectories are investigated for the trajectory design of future planetary moons explorers.

This paper presents the trajectory design for EQUilibriUm Lunar-Earth point 6U Spacecraft (EQUULEUS), which aims to demonstrate orbit control capability of CubeSats in the cislunar space. The mission plans to observe the far side of the Moon from an Earth-Moon L2 (EML2) libration point orbit. The EQUULEUS trajectory design needs to react to uncertainties of mission design parameters such as the launch conditions, errors, and thrust levels. The main challenge is to quickly design science orbits at EML2 and low-energy transfers from the post-deployment trajectory to the science orbits within the CubeSat’s limited propulsion capabilities. To overcome this challenge, we develop a systematic trajectory design approach that 1) designs over 13,000 EML2 quasi-halo orbits in a full-ephemeris model with a statistical stationkeeping cost evaluation, and 2) identifies families of low-energy transfers to the science orbits using lunar flybys and solar perturbations. The approach is successfully applied for the trajectory design of EQUULEUS.

Quasi-satellite orbits (QSOs) are stable retrograde orbits in the restricted three-body problem that have gained attention as a viable candidate for future deep-space missions towards remote planetary satellites. JAXA's robotic sample return mission MMX will utilize QSOs to perform scientific observations of the Martian moon Phobos before landing on its surface and attempt sample retrieval. The complex dynamical environment around Phobos makes the proximity operations of MMX quite challenging and requires novel and sophisticated techniques for maintaining and transferring between different quasi-satellite orbits. The present paper explores the application of invariant manifolds of unstable retrograde orbits to design transfer trajectories around Phobos. Starting from the equations of the Circular Hill Problem with ellipsoidal Phobos, we first compute families of three-dimensional QSOs using out-of-plane bifurca-tions near planar orbits. The feasibility of using unstable family members as staging orbits between high-altitude and low-altitude QSOs is later assessed. The final candidates are ranked based on MMX scientific requirements, transfer analyses, and station-keeping costs. It is found that intermediate 3D-QSOs can be maintained with as little as 1 m/s per month. Furthermore, it is discovered that transfer from high-altitude QSOs to low-altitude QSOs can be executed with a total ∆V of less than 40 m/s and total time of flight of less than 5 days.

Landing on Phobos and bringing samples from its surface would settle the debate on the origin of the Martian moons and support future manned exploration to Mars. To fulfill these scientific objectives, JAXA is planning to send a sample return probe to Phobos by the first half of the next decade, named the Martian Moons eXploration (MMX) mission. In order to explore scientifically interesting regions of Phobos, as well as to support the sampling operations of MMX, a number of Deployable CAMera 5 payloads are proposed to be deployed from quasi-satellite orbits (QSOs) around the Martian moon. This paper explores the feasibility of ballistic deployments from QSOs under realistic dynamical environment and surface constraints in order to guarantee surface settlement within the lifespan of DCAM5. First, we analyze the dynamical environment and escape speeds from Phobos by means of the Circular Hill Problem. Then, the surface coefficient of restitution is estimated by generic impacts onto Phobos regolith via discrete element method simulations. By combining these two analyses, maximum allowable impact velocities for surface settling are calculated and applied to downselect the number of feasible ballistic landings from QSOs. It is found that access to Phobos surface is possible especially along the leading and trailing sides of the Martian moon and in agreement with the engineering requirements of DCAM5.

Understanding the solar corona and its structure, evolution and composition can provide new insights regarding the processes that control the transport of energy throughout the solar atmosphere and out into the heliosphere. However, the visible emission coming from the corona is more than a million times weaker than the emission from the photosphere, implying that direct corona observations are only possible when the disk of the Sun is fully obscured. In this paper we perform a feasibility study of a Sun occultation mission using the Earth as a natural occulter. The challenge is that the occultation zone created by the Earth does not follow a Keplerian trajectory, causing satellites placed in this region to quickly drift away from eclipse conditions. To increase the number of revisits while optimizing the propellant budget, we propose optimal trajectories in the Sun–Earth-Spacecraft circular restricted three body problem that account for scientific and engineering constraints such as limited power budget and mission duration. Chemical propulsion, electric propulsion and solar sailing configurations are compared in terms of performance and mission feasibility, revealing how 24 h of corona observations would be possible every 39 days with as little as 199 m/s of í µí»¥í µí±. The feasibility of the solar sail approach is hereby demonstrated, making it a challenging engineering alternative to currently available technologies.

The relative motion about 4179 Toutatis is studied in order to investigate the feasibility of formation flying as an alternative concept for future asteroid exploration missions. In particular, the existence of quasi-frozen orbits about slowly rotating bodies allows us to compute families of periodic orbits in the body-fixed frame of the asteroid. Since these periodic orbits are of the center×center type, quasi-periodic invariant tori are calculated via fully numerical procedures and used to initialize spacecraft formations about the central body. Numerical simulations show that the resulting in-plane and out-of-plane relative trajectories remain bounded over long time spans; i.e., more than 30 days. •We consider two spacecraft flying in a formation about the asteroid 4179 Toutatis.•We derive first order conditions that minimize the drift between the satellites.•We find a full family of stable periodic orbits about the target asteroid.•We compute quasi-periodic invariant tori and use them to initialize the formation.•Numerical simulations prove that the satellites stay bounded for long time spans.

The Martian Moons eXploration (MMX) mission will study the Martian moons Phobos and Deimos, Mars, and their environments. The mission scenario includes both landing on the surface of Phobos to collect samples and deploying a small rover for in situ observations. Engineering safeties and scientific planning for these operations require appropriate evaluations of the surface environment of Phobos. Thus, the mission team organized the Landing Operation Working Team (LOWT) and Surface Science and Geology Sub-Science Team (SSG-SST), whose view of the Phobos environment is summarized in this paper. While orbital and large-scale characteristics of Phobos are relatively well known, characteristics of the surface regolith, including the particle size-distributions, the packing density, and the mechanical properties, are difficult to constrain. Therefore, we developed several types of simulated soil materials (simulant), such as UTPS-TB (University of Tokyo Phobos Simulant, Tagish Lake based), UTPS-IB (Impact-hypothesis based), and UTPS-S (Simpler version) for engineering and scientific evaluation experiments.

Quasi-satellite or distant retrograde orbits are stable orbits in the restricted three-body problem will be used in remote planetary satellite missions like JAXA’s Martian Moons eXploration (MMX). Due to Phobos’ irregular gravity field, the proximity operations of the MMX mission require sophisticated techniques for transfers and maintenance between different quasi-satellite orbits. This paper focuses on the design of transfer trajectories between Quasi-Satellite Orbits (QSO) around Phobos. We utilize horizontal bifurcated Multi-Revolution Periodic QSOs (MP-QSOs) that are found in the vicinity of an ellipsoidal Phobos under the assumptions of the circular Hill Problem. Transfer strategy via MP-QSOs is explored through transfer maps that illustrate transfer design space between different altitudes QSOs. It is found that transfers via MP-QSOs provide insights on minimum ΔV transfers and the parameters determining the transfer cost. The proposed transfer method is explicitly applied to MMX baseline QSOs. This research has identified that transfers between high(100 × 200 km)-to-mid(50 × 100 km), mid-to-low(30 × 50 km), and lower altitude(

Motivated by the near-future re-exploration of 1 the cislunar space, this paper investigates dynamical substitutes of the Earth-Moon’s resonant Near-Rectilinear Halo Orbits (NRHOs) under the Elliptic-Circular Restricted Four-Body Problem formulation of the Earth-Moon-Sun system. This model considers that the Earth and Moon move in elliptical orbits about each other and that a third body, the Sun, moves in a circular orbit about the Earth-Moon barycenter. By making use of this higher-fidelity dynamical model, we are able to incorporate the Sun’s influence and the Moon’s eccentricity, two of the most significant perturbations of the cislunar environment. As a result of these perturbations, resonant periodic NRHOs of the Earth-Moon Circular Restricted Three-Body Problem (CR3BP) are hereby replaced by two-dimensional quasi-periodic tori that better represent the dynamical evolution of satellites near the vicinity of the Moon. We present the steps and algorithms needed to compute these dynamical structures in the Elliptic-Circular model and subsequently assess their utility for spacecraft missions. We focus on the planned orbit for the NASA-led Lunar Gateway mission, a 9:2 synodic resonant L2 southern NRHO, as well as on the 4:1 synodic and 4:1 sidereal resonances, due to the proximity to the nominal orbit and their advantageous dynamical properties. We verify that the dynamical equivalents of these orbits preserve key dynamical attributes such as eclipse avoidance and near-linear stability. Furthermore, we find that the higher dimensionality of quasi-periodic solutions offers interesting alternatives to mission designers in terms of phasing maneuvers and low-altitude scientific observations.

The Martian Moons eXploration mission, currently under development by the Japan Aerospace Exploration Agency (JAXA), will be launched in 2024 with the goal of retrieving pristine samples from the surface of Phobos. Soon after arrival, the spacecraft will inject into retrograde relative trajectories known as quasi-satellite orbits and study the geophysical environment of the Martian moon for more than three years. This paper presents the orbit design and maintenance strategy of the Martian Moons eXploration mission in the framework of the elliptic Hill problem with ellipsoidal secondary. This paper first introduces a numerical continuation procedure on the eccentricity of Phobos to replace purely periodic solutions with families of quasi-periodic invariant tori. Two-dimensional torus maps can be then constructed and used to represent physical quantities of interest, as well as to generate reference trajectories at arbitrary epochs. Sensitivity and stability analyses are carried out to investigate the dynamic properties of retrograde relative trajectories in the elliptic case. Finally, a linear quadratic regulator is implemented in order to assess the robustness of the computed trajectories under injection, navigation, and execution errors. Monte Carlo simulations demonstrate that the baseline quasi-satellite orbits of the Martian Moons eXploration mission can be maintained with as low as 6.265 m/s per month.

Closed-loop feedback-driven control laws can be used to solve low-thrust many-revolution trajectory design and guidance problems with minimal computational cost. Lyapunov-based control laws offer the benefits of increased stability whilst their optimality can be increased by tuning their parameters. In this paper, a reinforcement learning framework is used to make the parameters of the Lyapunov-based Q-law state-dependent, increasing its optimality. The Jacobian of these state-dependent parameters is available analytically and, unlike in other optimisation approaches, can be used to enforce stability throughout the transfer. The results focus on GTO–GEO and LEO–GEO transfers in Keplerian dynamics, including the effects of eclipses. The impact of the network architecture on the behaviour is investigated for both time- and mass-optimal transfers. Robustness to navigation errors and thruster misalignment is demonstrated using Monte Carlo analyses. The resulting approach offers potential for on-board autonomous transfers and orbit reconfiguration.

Although recent numerical studies have demonstrated the possibility to obtain entire families of bounded relative trajectories in the Earth zonal problem, it is still unclear whether such initial conditions can cope with mismodeled dynamics such as atmospheric drag, luni-solar attraction, and solar radiation pressure. In the attempt to deal with the dynamical perturbations caused by the Earth's atmosphere, this paper offers a semi-analytical approach that combines previous work developed by the authors with newly found analytical relationships describing the effects of drag on the nodal and sidereal periods of a satellite. It is shown that both of these quantities change in a secular fashion due to the variations in the semi-major axis of spacecraft. Nevertheless, bounded relative motion can be guaranteed over the course of several days by either initializing the satellites on the same invariant curve of a stroboscopic mapping or by design of appropriate ballistic coefficients. These conclusions are supported by numerical simulations that prove the feasibility of our approach for spacecraft formations in low Earth orbit.

Air-breathing electric propulsion (ABEP) enables long duration missions at very low orbital altitudes through the use of drag compensation. A system-level spacecraft model is developed, using the interaction between thruster, intake and solar arrays, and coupled to a calculation of the drag. A quadratic solution is found for specific impulse and evaluated to identify the thruster performance required for drag-compensation at varying altitudes. An upper altitude limit around 190 km is based on a minimum thruster propellant density, resulting in required thruster performance values of 𝐼𝑠𝑝 > 3000 s and 𝑇 ∕ 𝑃 > 8 mN/kW for a realistic ABEP spacecraft. The orbit of an air-breathing spacecraft is propagated with time, which highlights the prescribed orbit eccentricity due to non-spherical gravity and therefore an increased variability in the atmospheric conditions. A thruster control law is introduced which avoids a divergent altitude behaviour by preventing thruster firings around the orbit periapsis, as well as adding robustness against atmospheric changes due to season and solar activity. Through the use of an initial frozen orbit, thruster control and an augmented 𝑇 ∕ 𝑃 , a stable long-term profile is demonstrated based on the performance data of a gridded-ion thruster tested with atmospheric propellants. An initial mean semi-major axis altitude of 200 km relative to the equatorial Earth radius, a spacecraft mass of 200 kg, 𝐼𝑠𝑝 = 5455 s and 𝑇 ∕ 𝑃 = 23 mN/kW, results in an altitude range of around 10 km at altitudes of 160–183 km during a period of medium to high solar activity.

Missions around small bodies present many challenges from their design to the operations, due to the highly non-linear and uncertain dynamics, the limited ∆v budget and constraints coming from orbit determination and mission design. Within this context, mathematical tools to enhance the understanding of the dynamics behavior can be proven useful to support the mission design process. Chaos indicators are adopted to reveal patterns of time-dependent dynamical systems and to enable the identification of practical stability regions, which are then exploited to design bounded orbits in the proximity of small bodies. The methodology is applied to study the MMX and Hera missions. In the MMX context, the final goal is to obtain bounded orbits useful for the global surface mapping and gravity potential determination of Phobos. On the other hand, concerning the Hera mission, a qualitative analysis of the natural motion about the Didymos binary asteroid system is carried out to compute bounded orbits convenient for the global characterization of the two asteroids and to investigate potential landing trajectories. Sensitivity analyses via Monte Carlo simulations are performed to prove the robustness of the different bounded orbits.

The risk of collisions in Earth’s orbit is growing markedly. In January 2021, SpaceX and OneWeb released an operator-to-operator fact sheet that highlights the critical reliance on conjunction data messages (CDMs) and observations, demonstrating the need for a diverse sensing environment for orbital objects. Recently, the University of Oxford and the University of Surrey developed, in collaboration with Trillium Technologies and the European Space Operations Center, an opensource Python package for modeling the spacecraft collision avoidance process, called Kessler. Such tools can be used for importing/exporting CDMs in their standard format, modeling the current low-Earth orbit (LEO) population and its short-term propagation from a given catalog file, as well as modeling the evolution of conjunction events based on the current population and observation scenarios, hence emulating the CDMs generation process of the Combined Space Operations Center (CSpOC). The model also provides probabilistic programming and ML tools to predict future collision events and to perform Bayesian inference (i.e., optimal use of all available observations). In the framework of a United Kingdom Space Agency-funded project, we analyze and study the impact of megaconstellations and observation models in the collision avoidance process. First, we monitor and report how the estimated collision risk and other quantities at the time of closest approach (e.g. miss distance, uncertainties, etc.) vary, according to different observation models, which emulate different radar observation accuracy. Then, we analyze the impact of future megaconstellations on the number of warnings generated from the increase in the number of conjunctions leading to an increased burden on space operators. FCC licenses were used to identify credible megaconstellation sources to understand how a potential consistent increase in active satellites will impact LEO situational safety. We finally present how our simulations help understand the impact of these future megaconstellations on the current population, and how we can devise better ground observation strategies to quantify future observation needs and reduce the burden on operators.

To maintain the periodic orbits in a three-body regime, a high-order Target Phase Approach (TPhA) is proposed in this work. Two types of polynomial maps, the phase-angle Poincaré map and high-order maneuver map, are established respectively for the determination of stationkeeping epochs and calculation of correction maneuvers. A stochastic optimization framework tailored for the TPhA-based stationkeeping process is leveraged in search of fuel-optimal and error-robust TPhA parameters. Quasi-Satellite Orbits (QSOs) around Phobos are investigated to demonstrate the efficacy of this approach in both low-and high-fidelity models. Monte-Carlo simulations demonstrate that the baseline QSO of JAXA's Martian Moons eXploration (MMX) mission can be maintained with a monthly manuever budget of around 1.13m/s.

Upcoming missions towards remote planetary moons will fly towards chaotic dynamical environments that are significantly perturbed by the oblateness of the host planet. This paper introduces a new time-periodic set of equations of motion that is based on the analytical solution of the zonal equatorial problem. Such a system, hereby referred to as the Zonal Hill Problem, remains populated by resonant periodic orbits and families of two-dimensional quasi-periodic invariant tori that are calculated by means of homotopy continuation procedures. The resulting periodic and quasi-periodic trajectories are investigated for the mission analysis and design of future planetary moons explorers.

Stable Quasi-Satellite Orbits (QSOs) have gained a lot of attention as suitable candidate orbits to explore remote planetary moons. Despite many studies on QSO and its orbital stability are found in literature, the problem of efficient transfer between two QSO is still unsolved. Previous works on transfers between two QSOs includes transfers between planar QSOs using impulsive maneuvers and utilizing the bifurcated multiple-revolutional periodic QSOs (MP-QSOs). Purpose of this paper is to explore the three-dimensional case and design transfer trajectories from lower altitude QSOs to three-dimensional QSOs being considered for upcoming sample return missions like Martian Moons eXploration (MMX) and PHOOTPRINT. Specifically, we utilize the unstable 3D QSO to generate manifolds and connect them to the planar QSO in the framework of Mars-Phobos Circular Hill Problem with Ellipsoidal secondary (HPE). Numerical simulations suggest that transfer between a planar and spatial QSO around Phobos is possible with ΔV as low as 8.277 m/s with TOF of 4.19 days and similarly transfer between a spatial and planar QSO with ΔV of 8.286 m/s with TOF of 6.75 days. As a result, transfer design space between planar and spatial QSOs is established using the invariant manifolds.

Near Rectilinear Halo Orbits (NRHOs) are orbits of great interest for the upcom-ing lunar missions. To maintain NRHOs in a three-body regime, a stationkeeping strategy based on a high-order Target Point Approach (TPA) is proposed, where fuel-optimal and error-robust TPA parameters are acquired from stochastic global optimization. Accurate TPA manevuers are calculated in a high-order fashion enabled by Differential Algebra (DA) techniques. Stochasticity is handled by incorporating Monte Carlo simulations in the process of optimization and the evaluation of high-order ODE expansions is employed to supplant the time-consuming numerical integration. Multiple candidate NRHOs with different stability properties are investigated.

Space exploration has often benefitted from the qualitative analyses of non integrable problems enabled by numerical continuation procedures. Yet, standard approaches based on Newton’s method typically end with discrete representations of family branches that may be subject to misinterpretation and overlook important dynamical features. In this research, we introduce novel continuation procedures based on the differential algebra of Taylor polynomials. Our algorithms aim at generating dense family branches as an atlas of polynomial charts that are locally valid for a range of system and continuation parameters. Examples of particular solutions will be shown within the framework of the Circular Restricted Three-Body Problem, along with fold and period-doubling bifurcations that are efficiently detected using automatic domain splitting and map inversion techniques.

Many Earth-Moon libration point orbits are being evaluated as candidates for future space missions, ranging from lunar gateways and possibly inhabited assets in the vicinity of the Moon to science missions in resonant synodic orbits, as in the cases of EQUULEUS and LUMIO CubeSat missions. The problem of libration point orbit maintenance is investigated here using the Target Point Approach, which leads to a closed LQR formulation of the station-keeping Delta-v once a set of parameters is specified. These parameters strongly depend on the orbit type. In view of the dramatic variety of stability property featured by halo orbits, these parameters must be fine-tuned with ad-hoc Monte-Carlo simulations or via a trial and error procedure. In this paper, a genetic algorithm is used to optimize the parameters of the target point approach for halo orbits with period ranging from 7 to 14 days and to provide a unified framework for orbit maintenance analysis of halo orbits. Numerical simulations in a high-fidelity ephemeris model show reductions in the station-keeping Delta-v budget for libration point orbits compared to previous results found in the literature.

Quasi-satellite orbits (QSO) are stable retrograde parking orbits around Phobos that are currently being considered for JAXA's upcoming robotic sample return mission MMX. During the proximity operations of MMX, the spacecraft inserted in a high altitude QSO will gradually descend to lower altitude QSOs with suitable transfer and station-keeping techniques between different relative QSOs. Preliminary analysis of two-impulsive planar transfers between relative retrograde orbits utilizing the bifurcated QSOs families is studied to estimate the ∆V costs and time of flights of the transfers. In spatial transfer problem, trajectories utilizing the invariant manifolds of unstable 3D-QSOs are weakly to highly unstable and require additional station-keeping strategies to perform MMX scientific observations. These transfer trajectories have a longer flight time and might need minor correction maneuvers along the transfer paths. In this paper, an orbital maintenance strategy that suppresses and eliminates linear dynamical instability of the unstable 3D-QSOs has been considered for shortlisting feasible 3D-QSOs for high-latitude observations. Differently from previous works, we utilize the initial guesses found through the preliminary results that provide two-impulsive transfer ∆V execution points and optimize the transfers between relative QSOs around Phobos. Primer vector theory is applied to investigate the primer vector of the transfer trajectories to evaluate whether intermediate maneuver or initial/final coasting times along the trajectories can minimize the total ∆V cost between the transfers.

Fifty years after the first images of the Martian moons were downlinked from the Mariner 7 and 9, Phobos and Deimos remain mysterious objects of the Solar System. On the one hand, the near-equatorial near-circular orbits of these remote bodies suggest that Phobos and Deimos are likely the product of a giant impact between a protoplanetary object and Mars' ancestor. On the other hand, the spectral and geophysical features observed by several spacecraft missions indicate that both of the Martian moons are extremely dark asteroids captured by the gravity well of the Red planet. The goal of the Martian Moons eXploration (MMX) mission–currently under development by JAXA and international collaborators–is to finally settle the debate on the moons' origin by retrieving pristine samples from the surface of Phobos. The samples will be collected in 2027 after extensive observation campaigns carried out from low-altitude retrograde orbits around the Martian satellite. This paper presents the orbit design and maintenance strategy for the proximity phase of MMX. We will start by reviewing the equations of the elliptical Hill problem in order to account for irregular gravity field of Phobos. The resulting non-autonomous system is populated by families of quasi-periodic orbit invariant tori that are generated through numerical continuation procedures on stroboscopic mappings. Five baseline trajectories are selected with different altitude profiles to meet the scientific instrument requirements of the MMX mission. Each of the candidate trajectories is tested under navigation errors and mismodeled dynamics to assess the operational feasibility of the newly found solutions and identify points of minimum and maximum sensitivity. A two-burn impulse per day strategy is implemented in order to track the reference retrograde orbits. Preliminary results show that the proposed relative trajectories can be maintained around Phobos for more than 30 days with less than 10 m/s.

Future missions to the Moon and beyond are likely to involve low-thrust propulsion technologies due to their propellant efficiency. However, these still present a difficult trajectory design problem. Lyapunov control laws can generate sub-optimal trajectories with minimal computational cost and are suitable for feasibility studies and as initial guesses for optimisation methods. In this work we combine Lyapunov control laws with state-dependent weights trained via reinforcement learning to design low-thrust transfers from GTO towards low-altitude Lunar orbits. The agent is able to explore third-body effects during training and learn to remain stable to perturbations during the different transfer phases. Three different approaches are investigated: backwards propagation, backwards propagation with freed geometry, and forwards propagation including rendezvous capability with the Lunar SOI. The last of these proves to be the most successful, coming within 6.6% of the optimal solution.

Despite the advantages of very-low altitude retrograde orbits around Phobos, questions remain about the efficacy of conventional station-keeping strategies in preventing spacecraft such as the Martian Moons eXploration from escaping or impacting against the surface of the small irregular moon. This paper introduces new high-fidelity simulations in which the output of a sequential Square-Root Information Filter is combined with recently developed orbit maintenance strategies based on differential algebra and convex optimization methods. The position and velocity vector of the spacecraft are first estimated using range, range-rate, and additional onboard data types such as LIDAR and camera images. This information is later processed to assess the necessity of an orbit maintenance maneuver based on the estimated relative altitude of MMX about Phobos. If a maneuver is deemed necessary, the state of the spacecraft is fed to either a successive convex optimization procedure or a high-order target phase approach capable of providing sub-optimal station-keeping maneuvers. The performance of the two orbit maintenance approaches is assessed via Monte Carlo simulations and compared against work in the literature so as to identify points of strength and weaknesses.

The Martian Moons eXploration mission will be the first of its kind to sample and study Mars's moon Phobos for a prolonged period of time. The aim of this work is to show that the adoption of periodic and quasi-periodic retrograde trajectories would be beneficial for the scientific return of MMX. A consider covariance analysis is hereby implemented in order to compare the estimation of high-order gravitational field coefficients from different orbital geometries and processing different sets of observables. It is shown that low-altitude non-planar quasi-satellite orbits would refine the knowledge of the moon's gravity field.

The application of dynamical systems theory in astrodynamics has enabled mission designers to construct innovative fuel-efficient trajectories within chaotic systems. While these solutions often take the form of periodic orbits and two-dimensional invariant manifolds, we can extend these solutions to higher dimensional objects such as quasi-periodic trajectories and their manifolds. Quasi-periodic trajectories allow us to better explore the dynamical environment and expand the design domain of spacecraft missions. By increasing the dimensional space of our solutions, we can obtain trajectories that better fulfill scientific objectives in space missions and that offer significant advantages in orbit maintenance and control. Additionally, these solutions are closer to what we would find in high-fidelity or full-ephemeris models, where perfect periodic solutions do not exist. This is particularly true for space missions flying towards chaotic environments, such as JAXA's Martian Moons eXploration (MMX) mission, which aims to retrieve samples from the largest moon of Mars, Phobos. This paper explores advances in the numerical computation of quasi-periodic tori families and possible applications to the MMX mission around the Mars-Phobos system under the formulation of a restricted three-body model with an ellipsoidal secondary.

This paper investigates the dynamical substitutes of the Moon's synodic and side-real resonant Near-Rectilinear Halo Orbits (NRHOs) under the Circular-Elliptic Restricted Four-Body Problem formulation. This model considers that the Earth and Moon move in elliptical orbits about each other and that a third body, the Sun, moves in a circular orbit about the Earth-Moon barycenter. The resonant periodic NRHOs are replaced by two-dimensional quasi-periodic tori in this model, which better approximate the real dynamics. We present the steps and algorithms needed to compute these dynamical structures as well as their geometry in the Circular-Elliptic model.

Understanding the solar corona and its composition can provide new insights regarding the temperature and the magnetic field of the Sun. The light coming from the corona is more than a million of times weaker than the direct light from the Sun; consequently observing the corona is only possible when the Sun is obscured. From ground, total solar eclipses offer a good opportunity to observe the corona; however, these events only occur every 18 months on average, lasting typically only for a few minutes. The goal of this paper is to perform a feasibility analysis of a Sun occultation mission using Earth as an occulter. However, the occultation zone created by the Earth does not follow a Keplerian trajectory, causing satellites placed in this region to quickly drift apart from the target area. To increase the number of revisits while optimizing the propellant budget, we propose optimal trajectories in the Sun-Earth-Spacecraft circular restricted three body problem that account for scientific and engineering constraints such as limited power budget and mission duration. Chemical Propulsion, Electric Propulsion and Solar Sailing configurations are compared in terms of performance and mission feasibility, revealing how 20 hours of corona observations per cycle would be possible with 0.25 km/s with a revisit of the occultation zone every 35 days. In addition to that the solar sail was proven to be an interesting alternative to chemical and low-thrust propulsion systems.

Closed-loop feedback-driven control laws can be used to solve low-thrust many-revolution trajectory design and guidance problems with minimal computational cost. They treat the problem from a targeting perspective and hence value stability over optimality. The optimality can be increased by making the parameters state-dependent at the cost of reduced stability. In this paper, an actor-critic reinforcement learning framework is used to make the parameters of the Lyapunov-based Q-law state-dependent. A single-layer neural network ensures the Jacobian of these state-dependent parameters can be calculated and used to enforce stability throughout the transfer. The current results focus on GTO-GEO and LEO-GEO transfers in Keplerian dynamics. A trade-off between optimality and stability is observed for the first, but the added stability increases optimality for the later. Robustness to uncertainties in position and velocity are also investigated, along with the effects of eclipses and dynamical perturbations such as J2, Sun and Moon third body attractions.