Search

The interaction between a Coronal Mass Ejection (CME) and a comet has been observed several times by in-situ observations from the Rosetta Plasma Consortium (RPC), which is designed to investigate the cometary magnetosphere of comet 67P/Churyumov-Gerasimenko (CG). Goetz et al. (2019) reported a magnetic field of up to 300 nT measured in the inner coma, which is among the largest interplanetary magnetic fields observed in the solar system. They suggested the large magnetic field observations in the inner coma come from magnetic field pile-up regions, which are generated by the interaction between a CME and/or corotating interaction region and the cometary magnetosphere. However, the detailed interaction between a CME and the cometary magnetosphere of comet CG in the inner coma has not been investigated by numerical simulations yet. In this manuscript, we will use a numerical model to simulate the interaction between comet CG and a Halloween class CME and investigate its magnetospheric response to the CME. We find that the plasma structures change significantly during the CME event, and the maximum value of the magnetic field strength is more than 500nT close to the nucleus. Virtual satellites at similar distances as Rosetta show that the magnetic field strength can be as large as 250nT, which is slightly less than what Goetz et al. (2019) reported., Comment: 19 pages, 5 figures, submitted to ApJ

Deriving properties of cometary nuclei from coma data is of significant importance for our understanding of cometary activity and has implications beyond. Ground-based data represent the bulk of measurements available for comets. Yet, to date these observations only access a comet's gas and dust coma at rather large distances from the surface and do not directly observe its surface. In contrast, spacecraft missions are one of the only tools that gain direct access to surface measurements. However, such missions are limited to roughly one per decade. We can overcome these challenges by recognising that the coma contains information about the nucleus's properties. In particular, the near-surface gas environment is most representative of the nucleus. It can inform us about the composition, regionality of activity, and sources of coma features and how they link to the topography, morphology, or other surface properties. The inner coma data is a particularly good proxy because it has not yet, or only marginally, been contaminated by coma chemistry or secondary gas sources. Additionally, the simultaneous observation of the innermost coma with the surface provides the potential to make a direct link between coma measurements and the nucleus. If we hope to link outer coma measurements obtained by Earth-based telescopes to the surface, we must first understand how the inner coma measurements are linked to the surface. Numerical models that describe the flow from the surface into the immediate surroundings are needed to make this connection. This chapter focuses on the advances made to understand the flow of the neutral gas coma from the surface to distances up to a few tens of nuclei radii. We describe both simple/heuristic models and state-of-the-art physically consistent models. Model limitations and what they each are best suited for are discussed., Comment: 29 pages, 9 figures, Chapter in press for the book Comets III, edited by K. Meech and M. Combi, University of Arizona Press

The need for quantitative characterization of the solar energetic particle (SEP) dynamics goes beyond being an academic discipline only. It has numerous practical implications related to human activity in space. The terrestrial magnetic field shields the International Space Station (ISS) and most uncrewed missions from exposure to SEP radiation. However, extreme SEP events with hard energy spectra are particularly rich in hundreds of MeV to several GeV protons that can reach the altitudes of the Low Earth Orbit (LEO). These protons have a high penetrating capability, thus producing significant radiation hazards for human spaceflight. SEPs also have a significant effect on the atmosphere. Sudden ionization of the upper atmosphere at high latitudes that occurs during polar cap absorption (PCA) events can block high frequency (HF) communication for hours, affecting communication with aircraft on intercontinental high-altitude flights. Another effect of SEPs in the atmosphere is creating NOx molecules in the upper atmosphere that can deplete the atmospheric ozone population. The paper also presents an analysis of (1) how various pitch angle diffusion coefficient approximations affect the properties of the simulated SEPs population and (2) discusses how pitch angle scattering when SEPs are beyond 1 AU affects a SEP event decay phase at the Earth's orbit.

Exoplanets orbiting M-dwarfs within habitable zones are exposed to stellar environments more extreme than that terrestrial planets experience in our Solar System, which can significantly impact the atmospheres of the exoplanets and affect their habitability and sustainability. This study provides the first prediction of hot oxygen corona structure and the associated photochemical loss from a 1 bar CO2-dominated atmosphere of a Venus-like rocky exoplanet, where dissociative recombination of O2+ ions is assumed to be the major source reaction for the escape of neutral O atoms and formation of the hot O corona (or exospheres) as on Mars and Venus. We employ a 3D Monte Carlo code to simulate the exosphere of Proxima Centauri b (PCb) based on the ionosphere simulated by a 3D magnetohydrodynamic model. Our simulation results show that variability of the stellar wind dynamic pressure over one orbital period of PCb does not affect the overall spatial structure of the hot O corona but contributes to the change in the global hot O escape rate that varies by an order of magnitude. The escape increases dramatically when the planet possesses its intrinsic magnetic fields as the ionosphere becomes more extended with the presence of a global magnetic field. The extended hot O corona may lead to a more extended H exosphere through collisions between thermal H and hot O, which exemplifies the importance of considering nonthermal populations in exospheres to interpret future observations., Comment: Accepted for publication in The Astrophysical Journal

Space plasma simulations have seen an increase in the use of magnetohydrodynamic (MHD) with embedded Particle-in-Cell (PIC) models. This combined MHD-EPIC algorithm simulates some regions of interest using the kinetic PIC method while employing the MHD description in the rest of the domain. The MHD models are highly efficient and their fluid descriptions are valid for most part of the computational domain, thus making large-scale global simulations feasible. However, in practical applications, the regions where the kinetic effects are critical can be changing, appearing, disappearing and moving in the computational domain. If a static PIC region is used, this requires a much larger PIC domain than actually needed, which can increase the computational cost dramatically. To address the problem, we have developed a new method that is able to dynamically change the region of the computational domain where a PIC model is applied. We have implemented this new MHD with Adaptively Embedded PIC (MHD-AEPIC) algorithm using the BATS-R-US Hall MHD and the Adaptive Mesh Particle Simulator (AMPS) as the semi-implicit PIC models. We describe the algorithm and present a test case of two merging flux ropes to demonstrate its accuracy. The implementation uses dynamic allocation/deallocation of memory and load balancing for efficient parallel execution. We evaluate the performance of MHD-AEPIC compared to MHD-EPIC and the scaling properties of the model to large number of computational cores.

MHD-based global space weather models have mostly been developed and maintained at academic institutions. While the "free spirit" approach of academia enables the rapid emergence and testing of new ideas and methods, the lack of long-term stability and support makes this arrangement very challenging. This paper describes a successful example of a university-based group, the Center of Space Environment Modeling (CSEM) at the University of Michigan, that developed and maintained the Space Weather Modeling Framework (SWMF) and its core element, the BATS-R-US extended MHD code. It took a quarter of a century to develop this capability and reach its present level of maturity that makes it suitable for research use by the space physics community through the Community Coordinated Modeling Center (CCMC) as well as operational use by the NOAA Space Weather Prediction Center (SWPC)., Comment: 105 pages, 36 figures, in press

Neutral hydrogen has been shown to greatly impact the plasma flow in the heliopshere and the location of the heliospheric boundaries. We present the results of the Solar-wind with Hydrogen Ion Exchange and Large-scale Dynamics (SHIELD) model, a new, self-consistent, kinetic-MHD model of the outer heliosphere within the Space Weather Modeling Framework. The charge-exchange mean free path is on order of the size of the heliosphere; therefore, the neutral atoms cannot be described as a fluid. The SHIELD model couples the MHD solution for a single plasma fluid to the kinetic solution from for neutral hydrogen atoms streaming through the system. The kinetic code is based on the Adaptive Mesh Particle Simulator (AMPS), a Monte Carlo method for solving the Boltzmann equation. The SHIELD model accurately predicts the increased filtration of interstellar neutrals into the heliosphere. In order to verify the correct implementation within the model, we compare the results of the SHIELD model to other, well-established kinetic-MHD models. The SHIELD model matches the neutral hydrogen solution of these studies as well as the shift in all heliospheric boundaries closer to the Sun in comparison the the multi-fluid treatment of the neutral hydrogen atoms. Overall the SHIELD model shows excellent agreement to these models and is a significant improvement to the fluid treatment of interstellar hydrogen., Comment: 21 pages, 6 figures. ApJS submitted

Simulation of plasmas in electromagnetic fields requires numerical solution of a kinetic equation that describes the time evolution of the particle distribution function. In this paper we propose a finite volume scheme based on integral relation for Poisson brackets to solve the Liouville equation, the most fundamental kinetic equation. The proposed scheme conserves the number of particles, maintains the total-variation-diminishing (TVD) property, and provides high-quality numerical results. Other types of kinetic equations may be also formulated in terms of Poisson brackets and solved with the proposed method including the transport equations describing the acceleration and propagation of Solar Energetic Particles (SEPs), which is of practical importance, since the high energy SEPs produce radiation hazards. The proposed scheme is demonstrated to be accurate and efficient, which makes it applicable to global simulation systems analyzing space weather., Comment: 50 pages, 9 figures

The Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) suite of instruments operated throughout the over two years of the Rosetta mission operations in the vicinity of comet 67P/Churyumov-Gerasimenko. It measured gas densities and composition throughout the comet's atmosphere, or coma. Here we present two-years' worth of measurements of the relative densities of the four major volatile species in the coma of the comet, H2O. CO2, CO and O2, by one of the ROSINA sub-systems called the Double Focusing Mass Spectrometer (DFMS). The absolute total gas densities were provided by the Comet Pressure Sensor (COPS), another ROSINA sub-system. DFMS is a very high mass resolution and high sensitivity mass spectrometer able to resolve at a tiny fraction of an atomic mass unit. We have analyzed the combined DFMS and COPS measurements using an inversion scheme based on spherical harmonics that solves for the distribution of potential surface activity of each species as the comet rotates, changing solar illumination, over short intervals and as the comet changes distance from the sun and orientation of its spin axis over long time intervals. We also use the surface boundary conditions derived from the inversion scheme to simulate the whole coma with our fully kinetic Direct Simulation Monte Carlo model and calculate the production rates of the four major species throughout the mission. We compare the derived production rates with revised remote sensing observations by the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) as well as with published observations from the Microwave Instrument for the Rosetta Orbiter (MIRO). Finally we use the variation of the surface production of the major species to calculate the total mass loss over the mission and, for different estimates of the dust/gas ratio, calculate the variation of surface loss over the nucleus., Comment: 144 pages, 9 figures, 8 tables, 1 supplemental table

The nucleus of the Jupiter-family comet 67P/Churyumov-Gerasimenko was discovered to be bi-lobate in shape when the European Space Agency spacecraft Rosetta first approached it in July 2014. The bi-lobate structure of the cometary nucleus has led to much discussion regarding the possible manner of its formation and on how the composition of each lobe might compare with that of the other. During its two-year-long mission from 2014 to 2016, Rosetta remained in close proximity to 67P/Churyumov-Gerasimenko, studying its coma and nucleus in situ. Based on lobe-specific measurements of HDO and H2O performed with the ROSINA DFMS mass spectrometer on board Rosetta, the Deuterium-to-Hydrogen ratios in water from the two lobes could be compared. No appreciable difference was observed, suggesting that both lobes formed in the same region and are homogeneous in their Deuterium-to-Hydrogen ratios., Comment: 10 pages, 6 figures. Submitted to MNRAS

In this Letter, we make use of sophisticated 3D numerical simulations to assess the extent of atmospheric ion and photochemical losses from Mars over time. We demonstrate that the atmospheric ion escape rates were significantly higher (by more than two orders of magnitude) in the past at $\sim 4$ Ga compared to the present-day value owing to the stronger solar wind and higher ultraviolet fluxes from the young Sun. We found that the photochemical loss of atomic hot oxygen dominates over the total ion loss at the current epoch whilst the atmospheric ion loss is likely much more important at ancient times. We briefly discuss the ensuing implications of high atmospheric ion escape rates in the context of ancient Mars, and exoplanets with similar atmospheric compositions around young solar-type stars and M-dwarfs., Comment: 6 pages, 4 figures, 1 table, accepted for publication in ApJ Letters

We study roles of the thermosphere and exosphere on the Martian ionospheric structure and ion escape rates in the process of the solar wind-Mars interaction. We employ a four-species multifluid MHD (MF-MHD) model to simulate the Martian ionosphere and magnetosphere. The $cold$ thermosphere background is taken from the Mars Global Ionosphere Thermosphere Model (M-GITM) and the $hot$ oxygen exosphere is adopted from the Mars exosphere Monte Carlo model - Adaptive Mesh Particle Simulator (AMPS). A total of four cases with the combination of 1D (globally averaged) and 3D thermospheres and exospheres are studied. The ion escape rates calculated by adopting 1D and 3D atmospheres are similar; however, the latter are required to adequately reproduce MAVEN ionospheric observations. In addition, our simulations show that the 3D hot oxygen corona plays an important role in preventing planetary molecular ions (O$_2^+$ and CO$_2^+$) escaping from Mars, mainly resulting from the mass loading of the high-altitude exospheric O$^+$ ions. The $cold$ thermospheric oxygen atom, however, is demonstrated to be the primary neutral source for O$^+$ ion escape during the relatively weak solar cycle 24., Comment: 21 pages, 10 figures, 4 tables, accepted for publication in Journal of Geophysical Research-Space Physics

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.