Search

Exoplanets with short orbit period reside very close to their host stars. They transition very rapidly between different sectors of the circumstellar space environment along their orbit, leading to large variations of the magnetic field in the vicinity of the planet on short timescales. This rapid change of the magnetic flux through the conducting and resistive layer of the planetary upper atmosphere may drive currents that dissipate in the form of Joule Heating. Here, we estimate the amount of Joule Heating dissipation in the upper atmosphere of Trappist-1e, and two hypothetical planets orbiting the Sun in close-in orbits. We find that the rapid orbital motion could drive a significant amount of atmospheric heating and could significantly affect the planetary atmosphere escape rate. Thus, the process should be accounted for when studying the long-term evolution of exoplanetary atmospheres., Comment: 8 pages, 4 figures, accepted to ApJ

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

The habitability of the surface of any planet is determined by a complex evolution of its interior, surface, and atmosphere. The electromagnetic and particle radiation of stars drive thermal, chemical and physical alteration of planetary atmospheres, including escape. Many known extrasolar planets experience vastly different stellar environments than those in our Solar system: it is crucial to understand the broad range of processes that lead to atmospheric escape and evolution under a wide range of conditions if we are to assess the habitability of worlds around other stars. One problem encountered between the planetary and the astrophysics communities is a lack of common language for describing escape processes. Each community has customary approximations that may be questioned by the other, such as the hypothesis of H-dominated thermosphere for astrophysicists, or the Sun-like nature of the stars for planetary scientists. Since exoplanets are becoming one of the main targets for the detection of life, a common set of definitions and hypotheses are required. We review the different escape mechanisms proposed for the evolution of planetary and exoplanetary atmospheres. We propose a common definition for the different escape mechanisms, and we show the important parameters to take into account when evaluating the escape at a planet in time. We show that the paradigm of the magnetic field as an atmospheric shield should be changed and that recent work on the history of Xenon in Earth's atmosphere gives an elegant explanation to its enrichment in heavier isotopes: the so-called Xenon paradox.

The search for exoplanets in the radio bands has been focused on detecting radio emissions produced by the interaction between magnetized planets and the stellar wind (auroral emission). Here we introduce a new tool, which is part of our MHD stellar corona model, to predict the ambient coronal radio emission and its modulations induced by a close planet. For simplicity, the present work assumes that the exoplanet is stationary in the frame rotating with the stellar rotation. We explore the radio flux modulations using a limited parameter space of idealized cases by changing the magnitude of the planetary field, its polarity, the planetary orbital separation, and the strength of the stellar field. We find that the modulations induced by the planet could be significant and observable in the case of hot Jupiter planets --- above 100% modulation with respect to the ambient flux in the $10-100~MHz$ range in some cases, and 2-10% in the frequency bands above $250~MHz$ for some cases. Thus, our work indicates that radio signature of exoplanets might not be limited to low-frequency radio range. We find that the intensity modulations are sensitive to the planetary magnetic field polarity for short-orbit planets, and to the stellar magnetic field strength for all cases. The new radio tool, when applied to real systems, could provide predictions for the frequency range at which the modulations can be observed by current facilities., Comment: Accepted to AJ

We compare a global high resolution resistive magnetohydrodynamics (MHD) simulation of Earth's magnetosphere with observations from the Magnetospheric Multiscale (MMS) constellation for a southward IMF magnetopause crossing during October 16, 2015 that was previously identified as an electron diffusion region (EDR) event. The simulation predicts a complex time-dependent magnetic topology consisting of multiple separators and flux ropes. Despite the topological complexity, the predicted distance between MMS and the primary separator is less than 0.5 Earth radii. These results suggest that global magnetic topology, rather than local magnetic geometry alone, determines the location of the electron diffusion region at the dayside magnetopause., Comment: 13 pages, 4 figures; Submitted to Journal of Geophysical Research

The Community Coordinated Modeling Center (CCMC) at NASA Goddard Space Flight Center is a multi-agency partnership to enable, support and perform research and development for next-generation space science and space weather models. CCMC currently hosts nearly 100 numerical models and a cornerstone of this activity is the Runs on Request (RoR) system which allows anyone to request a model run and analyze/visualize the results via a web browser. CCMC is also active in the education community by organizing student research contests, heliophysics summer schools, and space weather forecaster training for students, government and industry representatives. We present a generic magnetohydrodynamic (MHD) model - PAMHD - that has been added to the CCMC RoR system which allows the study of a variety of fluid and plasma phenomena in one, two and three dimensions using a dynamic point-and-click web interface. Flexible initial and boundary conditions allow experimentation with a variety of plasma physics problems such as shocks, instabilities, planetary magnetospheres and astrophysical systems. Experimentation with numerical effects, e.g. resolution, solution method and boundary conditions, is also possible and can provide valuable context for space weather forecasters when interpreting observations or modeling results. We present an overview of the C++ implementation and show example results obtained through the CCMC RoR system, including the first to our knowledge MHD simulation of the interaction of the magnetospheres of Jupiter and Saturn in two dimensions., Comment: 28 pages, 11 figures, submitted to Computer Physics Communications

We present high resolution Hall MHD simulations of Ganymede's magnetosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection layers have significant global effects not captured in resistive MHD simulations. Consistent with local kinetic simulations of magnetic reconnection, our global simulations show the development of intense field-aligned currents along the magnetic separatrices. These currents extend all the way down to the moon's surface, where they may contribute to Ganymede's aurora. Within the magnetopause and magnetotail current sheets, Hall currents in the reconnection plane accelerate ions to the local Alfv\'en speed in the out-of-plane direction, producing a global system of ion drift belts that circulates Jovian magnetospheric plasma throughout Ganymede's magnetosphere. We discuss some observable consequences of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian "double magnetopause" structure, an Alfv\'enic ion jet extending across the upstream magnetopause and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves., Comment: 14 pages, 12 figures; presented at Geospace Environment Modeling (GEM) workshop (June, 2014) and Fall American Geophysical Union (AGU) meeting (December, 2014); submitted to Journal of Geophysical Research, December 2014

Solar eruptions, known as Coronal Mass Ejections (CMEs), are frequently observed on our Sun. Recent Kepler observations of superflares on G-type stars have implied that so called super-CMEs, possessing kinetic energies 10 times of the most powerful CME event ever observed on the Sun, could be produced with a frequency of 1 event per 800-2000 yr on solar-like slowly rotating stars. We have performed a 3D time-dependent global magnetohydrodynamic simulation of the magnetic interaction of such a CME cloud with the Earth's magnetosphere. We calculated the global structure of the perturbed magnetosphere and derive the latitude of the open-closed magnetic field boundary. We also estimated energy fluxes penetrating the Earth's ionosphere and discuss the consequences of energetic particle fluxes on biological systems on early Earth., Comment: To be published in Proceedings of 18th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun Proceedings of Lowell Observatory (9-13 June 2014) Edited by G. van Belle & H. Harris

Energetic processes often produce transversely-heated angular distributions of the magnetized core (lowest energy) plasma. This characteristic is found in solar wind ion pickup, resulting from cometary or interstellar gas ionization, in Earths' ionosphere, and with hot ions formed around the Space Transportation System during gas releases. We investigate the thermalization of O+ ion pickup using the 2.5D hybrid simulation method (with fluid electrons and kinetic ions) of the ion pickup (ring) distributions, formed in the auroral ionosphere, with a range of ring velocities and thermal to magnetic pressure ratios. We find that in the unstable collisonless regime the anisotropy of the non-thermal distribution produces the ion-cyclotron instability, and the nonlinear relaxation is accompanied by wave-particle scattering that results in an emitted power of EMIC waves. We conclude that ionospheric pickup thermalization is slow due to the small ring speed compared to the thermal and Alfven speeds, while in the solar wind and other space plasmas regions with larger ion-ring velocity the collisionless relaxation and thermalization is rapid in terms of O+ ion gyro-period.

Many space plasmas (especially electrons generated in planetary ionospheres) exhibit fine-detailed structures that are challenging to fully resolve with the energy resolution of typical space plasma analyzers (10% → 20%). While analyzers with higher resolution have flown, generally this comes at the expense of sensitivity and temporal resolution. We present a new technique for measuring plasmas with extremely high energy resolution through the combination of a top-hat Electrostatic Analyzer (ESA) followed by an internally mounted Retarding Potential Analyzer (RPA). When high resolutions are not required, the RPA is grounded, and the instrument may operate as a typical general-purpose plasma analyzer using its ESA alone. We also describe how such an instrument may use its RPA to remotely vary the geometric factor (sensitivity) of a top hat analyzer, as was performed on the New Horizons Solar Wind at Pluto and MAVEN SupraThermal and Thermal Ion Composition instruments. Finally, we present results from laboratory testing of our prototype, showing that this technique may be used to construct an instrument with 1.6% energy resolution, constant over all energies and angles.

On 27 June 1996, the NASA Galileo spacecraft made humanity’s first flyby of Jupiter’s largest moon, Ganymede, discovering that it is the only moon known to possess an internally generated magnetic field. Resurrecting the original Galileo Plasma Subsystem (PLS) data analysis software, we processed the raw PLS data from G01 and for the first time present the properties of plasmas encountered. Entry into the magnetosphere of Ganymede occurred near the confluence of the magnetopause and plasma sheet. Reconnection-driven plasma flows were observed (consistent with an Earth-like Dungey cycle), which may be a result of reconnection in the plasma sheet, magnetopause, or might be Ganymede’s equivalent of a Low-Latitude Boundary Layer. Dropouts in plasma density combined with velocity perturbations afterward suggest that Galileo briefly crossed the cusps into closed magnetic field lines. Galileo then crossed the cusps, where field-aligned precipitating ions were observed flowing down into the surface, at a location consistent with observations by the Hubble Space Telescope. The density of plasma outflowing from Ganymede jumped an order of magnitude around closest approach over the north polar cap. The abrupt increase may be a result of crossing the cusp or may represent an altitude-dependent boundary such as an ionopause. More diffuse, warmer field-aligned outflows were observed in the lobes. Fluxes of particles near the moon on the nightside were significantly lower than on the dayside, possibly resulting from a diurnal cycle of the ionosphere and/or neutral atmosphere.

We examine the mechanisms responsible for the dropout of energetic electron flux during 31 May to 1 June 2013 using Van Allen Probe (Radiation Belt Storm Probes (RBSP)) electron flux data and simulations with the Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. During the storm main phase, L-shells at RBSP locations are greater than ~8, which are connected to open drift shells. Consequently, diminished electron fluxes were observed over a wide range of energies. The combination of drift shell splitting, magnetopause shadowing, and drift loss all results in butterfly electron pitch angle distributions (PADs) at the nightside. During storm sudden commencement, RBSP observations display electron butterfly PADs over a wide range of energies. However, it is difficult to determine whether there are butterfly PADs during the storm main phase since the maximum observable equatorial pitch angle from RBSP is not larger than ~40° during this period. To investigate the causes of the dropout, the CIMI model is used as a global 4-D kinetic inner magnetosphere model. The CIMI model reproduces the dropout with very similar timing and flux levels and PADs along the RBSP trajectory for 593 keV. Furthermore, the CIMI simulation shows butterfly PADs for 593 keV during the storm main phase. Based on comparison of observations and simulations, we suggest that the dropout during this event mainly results from magnetopause shadowing.

During the 14 November 2012 geomagnetic storm, the Van Allen Probes spacecraft observed a number of sharp decreases ('dropouts') in particle fluxes for ions and electrons of different energies. In this paper, we investigate the global magnetosphere dynamics and magnetosphere- ionosphere (M-I) coupling during the dropout events using multipoint measurements by Van Allen Probes, TWINS, and AMPERE together with the output of the two-way coupled global BATS-R-US-CRCM model. We find different behavior for two pairs of dropouts. For one pair, the same pattern was repeated: (1) weak nightside Region 1 and 2 Birkeland currents before and during the dropout; (2) intensification of Region 2 currents after the dropout; and (3) a particle injection detected by TWINS after the dropout. The model predicted similar behavior of Birkeland currents. TWINS low-altitude emissions demonstrated high variability during these intervals, indicating high geomagnetic activity in the near-Earth tail region. For the second pair of dropouts, the structure of both Birkeland currents and ENA emissions was relatively stable. The model also showed quasi-stationary behavior of Birkeland currents and simulated ENA emissions with gradual ring current buildup. We confirm that the first pair of dropouts was caused by large-scale motions of the OCB (open-closed boundary) during substorm activity. We show the new result that this OCB motion was associated with global changes in Birkeland (M-I coupling) currents and strong modulation of low-altitude ion precipitation. The second pair of dropouts is the result of smaller OCB disturbances not related to magnetospheric substorms. The local observations of the first pair of dropouts result from a global magnetospheric reconfiguration, which is manifested by ion injections and enhanced ion precipitation detected by TWINS and changes in the structure of Birkeland currents detected by AMPERE. This study demonstrates that multipoint measurements along with the global model results enable the reconstruction of a more complete system-level picture of the dropout events and provides insight into M-I coupling aspects that have not previously been investigated.

In this response, we address the three main comments by Tsurutani et al. (2018, http://doi.org/10.1002/2017JA024779) namely, unusually high plasma density, interplanetary magnetic field intensity, and fast storm recovery phase. The authors agree that there is room to improve the modeling by taking into account these comments and other aspects that were not fully explored during our initial work. We are already in the process of undertaking a more comprehensive modeling project.

Abstract Parallel electric fields and their associated electric potential structures play a crucial role inionospheric-magnetospheric interactions at any planet. Although there is abundant evidence that parallel electric fields play key roles in Martian ionospheric outflow and auroral electron acceleration, the fields themselves are challenging to directly measure due to their relatively weak nature. Using measurements by the Solar Wind Electron Analyzer instrument aboard the NASA Mars Atmosphere and Volatile EvolutioN(MAVEN) Mars Scout, we present the discovery and measurement of a substantial (Phi) Mars 7.7 +/-0.6 V) parallel electric potential drop on closed magnetic field lines spanning the terminator from day to night above the great impact basin of Utopia Planitia, a region largely free of crustal magnetic fields. A survey of the previous 26 orbits passing over a range of longitudes revealed similar signatures on seven orbits, with a mean potential drop (Phi) Mars of 10.9 +/- 0.8 V, suggestive that although trans-terminator electric fields of comparable strength are not ubiquitous, they may be common, at least at these northerly latitudes.

Atmospheres of exoplanets in the habitable zones around active young G-K-M stars are subject to extreme X-ray and EUV (XUV) (Extreme Ultraviolet) fluxes from their host stars that can initiate atmospheric erosion. Atmospheric loss affects exoplanetary habitability in terms of surface water inventory, atmospheric pressure, the efficiency of greenhouse warming, and the dosage of the UV surface irradiation. Thermal escape models suggest that exoplanetary atmospheres around active K-M stars should undergo massive hydrogen escape, while heavier species including oxygen will accumulate forming an oxidizing atmosphere. Here, we show that non-thermal oxygen ion escape could be as important as thermal, hydrodynamic H escape in removing the constituents of water from exoplanetary atmospheres under supersolar XUV irradiation. Our models suggest that the atmospheres of a significant fraction of Earth-like exoplanets around M dwarfs and active K stars exposed to high XUV fluxes will incur a significant atmospheric loss rate of oxygen and nitrogen, which will make them uninhabitable within a few tens to hundreds million years, given a low replenishment rate from volcanism or cometary bombardment. Our non-thermal escape models have important implications for the habitability of the Proxima Centauri's terrestrial planet.

Understanding what processes govern atmospheric escape and the loss of planetary water is of paramount importance for understanding how life in the universe can exist. One mechanism thought to be important at all planets is an ambipolar electric field that helps ions overcome gravity. We report the discovery and first quantitative extraterrestrial measurements of such a field at the planet Venus. Unexpectedly, despite comparable gravity, we show the field to be five times stronger than in Earths similar ionosphere. Contrary to our understanding, Venus would still lose heavy ions (including oxygen and all water-group species) to space, even if there were no stripping by the solar wind. We therefore find that it is possible for planets to lose heavy ions to space entirely through electric forces in their ionospheres and such an electric wind must be considered when studying the evolution and potential habitability of any planet in any star system.

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.

In this work, we demonstrate how to track magnetic separators in three-dimensional simulated magnetic fields with or without magnetic nulls, apply these techniques to enhance our understanding of reconnection at the magnetopause. We present three methods for locating magnetic separators and apply them to 3-D resistive MHD simulations of the Earth's magnetosphere using the Block-Adaptive-Tree Solar-wind Roe-type Upwind Scheme code. The techniques for finding separators and determining the reconnection rate are insensitive to interplanetary magnetic field (IMF) clock angle and can in principle be applied to any magnetospheric model. Moreover, the techniques have a number of advantages over prior separator finding techniques applied to the magnetosphere. The present work examines cases of high and low resistivity for two clock angles. We go beyond previous work examine the separator during Flux Transfer Events (FTEs). Our analysis of reconnection on the magnetopause yields a number of interesting conclusions: Reconnection occurs all along the separator even during predominately northward IMF cases. Multiple separators form in low-resistivity conditions, and in the region of an FTE the separator splits into distinct branches. Moreover, the local contribution to the reconnection rate, as determined by the local parallel electric field, drops in the vicinity of the FTE with respect to the value when there are none.

The precipitation of high-energy magnetospheric electrons (E greater than 500-600 electronvolts) in the diffuse aurora contributes significant energy flux into Earth's ionosphere. In the diffuse aurora, precipitating electrons initially injected from the plasmasheet via wave-particle interaction processes degrade in the atmosphere toward lower energies and produce secondary electrons via impact ionization of the neutral atmosphere. These initially precipitating electrons of magnetospheric origin can be additionally reflected back into the magnetosphere by the two magnetically conjugated atmospheres, leading to a series of multiple reflections that can greatly influence the initially precipitating flux at the upper ionospheric boundary (700-800 kilometers) and the resultant population of secondary electrons and electrons cascading toward lower energies. We present the solution of the Boltzmann.Landau kinetic equation that uniformly describes the entire electron distribution function in the diffuse aurora, including the affiliated production of secondary electrons (E is less than or equal to 600 electronvolts) and their energy interplay in the magnetosphere and two conjugated ionospheres. This solution takes into account the role of multiple atmospheric reflections of the precipitated electrons that were initially moved into the loss cone via wave.particle interaction processes in Earth's plasmasheet.