Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by di fferent physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and e fficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solarwind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a blockadaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional blockadaptive grid. We describe the algorithms of BATL and demonstrate its e fficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that di fferent time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.https://ntrs.nasa.gov/search.jsp?R=20110005631 2018-05-12T08:50:45+00:00Z

Nitrogen is a critical ingredient of complex biological molecules 1. Molecular nitrogen, however, which was outgassed into the Earth's early atmosphere 2 , is relatively chemically inert and nitrogen fixation into more chemically reactive compounds requires high temperatures. Possible mechanisms of nitrogen fixation include lightning, atmospheric shock heating by meteorites, and solar ultraviolet radiation 3,4. Here we show that nitrogen fixation in the early terrestrial atmosphere can be explained by frequent and powerful coronal mass ejection events from the young Sun-so-called superflares. Using magnetohydrodynamic simulations constrained by Kepler Space Telescope observations, we find that successive superflare ejections produce shocks that accelerate energetic particles, which would have compressed the early Earth's magnetosphere. The resulting extended polar cap openings provide pathways for energetic particles to penetrate into the atmosphere and, according to our atmospheric chemistry simulations, initiate reactions converting molecular nitrogen, carbon dioxide and methane to the potent greenhouse gas nitrous oxide as well as hydrogen cyanide, an essential compound for life. Furthermore, the destruction of N 2 , CO 2 and CH 4 suggests that these greenhouse gases cannot explain the stability of liquid water on the early Earth. Instead, we propose that the e cient formation of nitrous oxide could explain a warm early Earth. Here we develop a new concept for the rise of prebiotic chemistry on early Earth that suggests abiotic nitrogen fixation mediated by the energy flux from palaeo-solar eruptive events. The flare statistics of Kepler data suggest that the frequency of occurrence of superflares with energies >5 × 10 34 erg observed on G-type dwarfs follows a power-law distribution with a spectral index between α = −2.0 and −2.3, which is comparable to those observed on active M-type red-dwarf stars and the Sun 5,6. If the occurrence rate of superflares on young solar-like stars is ∼0.1 events per day 6 , then the frequency of super Carrington-type flare events with energy

Atmospheres of exoplanets in the habitable zones around active young G-K-M stars are subject to extreme X-ray and EUV (XUV) 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 of Myr, 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.

[1] Ionospheric outflow has been shown to be a significant contributor to the plasma population of the magnetosphere during active geomagnetic conditions. We present the results of new efforts to model the source and effects of out-flowing plasma in the Space Weather Modeling Framework (SWMF). In particular, we develop and use the Polar Wind Outflow Model (PWOM), a field-aligned, multifluid, multifield line polar wind code to simulate the ionospheric outflow. The PWOM is coupled to the ionosphere electrodynamics and global magnetosphere components of the SWMF, so we can calculate the outflow and its resulting impact on magnetospheric composition and dynamics. By including the outflow as part of a coupled system, we study the consequences of outflow on the larger space environment system. We present our methodology for the magnetosphere-ionosphere coupling, as well as the effect of outflow on the magnetosphere during two geomagnetic storms. Moreover, we explore the use of multispecies MHD to track the resulting plasma composition in the magnetosphere. We find that, by including ionospheric outflow during geomagnetic storms, we can reduce the RMS error in the simulated magnetic field as compared with various GOES satellites by as much as 50%. Additionally, we find that the outflow causes a strong decrease in Dst and in the cross-polar cap potential.

Simulation studies of the Earth's radiation belts and ring current are very useful in understanding the acceleration, transport, and loss of energetic particles. Recently, the Comprehensive Ring Current Model (CRCM) and the Radiation Belt Environment (RBE) model were merged to form a Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. CIMI solves for many essential quantities in the inner magnetosphere, including ion and electron distributions in the ring current and radiation belts, plasmaspheric density, Region 2 currents, convection potential, and precipitation in the ionosphere. It incorporates whistler mode chorus and hiss wave diffusion of energetic electrons in energy, pitch angle, and cross terms. CIMI thus represents a comprehensive model that considers the effects of the ring current and plasmasphere on the radiation belts. We have performed a CIMI simulation for the storm on 5-9 April 2010 and then compared our results with data from the Two Wide-angle Imaging Neutral-atom Spectrometers and Akebono satellites. We identify the dominant energization and loss processes for the ring current and radiation belts. We find that the interactions with the whistler mode chorus waves are the main cause of the flux increase of MeV electrons during the recovery phase of this particular storm. When a self-consistent electric field from the CRCM is used, the enhancement of MeV electrons is higher than when an empirical convection model is applied. We also demonstrate how CIMI can be a powerful tool for analyzing and interpreting data from the new Van Allen Probes mission.

The precipitation of high‐energy magnetospheric electrons (E ∼ 600 eV–10 KeV) in the diffuse aurora contributes significant energy flux into the Earth's ionosphere. To fully understand the formation of this flux at the upper ionospheric boundary, ∼700–800 km, it is important to consider the coupled ionosphere‐magnetosphere system. In the diffuse aurora, precipitating electrons initially injected from the plasma sheet via wave‐particle interaction processes degrade in the atmosphere toward lower energies and produce secondary electrons via impact ionization of the neutral atmosphere. These precipitating electrons can be additionally reflected upward from the two conjugate ionospheres, leading to a series of multiple reflections through the magnetosphere. These reflections greatly influence the initially precipitating flux at the upper ionospheric boundary (700–800 km) and the resultant population of secondary electrons and electrons cascading toward lower energies. In this paper, we present the solution of the Boltzman‐Landau kinetic equation that uniformly describes the entire electron distribution function in the diffuse aurora, including the affiliated production of secondary electrons (E < 600 eV) and its energy interplay in the magnetosphere and two conjugated ionospheres. This solution takes into account, for the first time, the formation of the electron distribution function in the diffuse auroral region, beginning with the primary injection of plasma sheet electrons via both electrostatic electron cyclotron harmonic waves and whistler mode chorus waves to the loss cone, and including their subsequent multiple atmospheric reflections in the two magnetically conjugated ionospheres. It is demonstrated that magnetosphere‐ionosphere coupling is key in forming the electron distribution function in the diffuse auroral region.

[1] The magnetosphere contains a significant amount of ionospheric O + , particularly during geomagnetically active times. The presence of ionospheric plasma in the magnetosphere has a notable impact on magnetospheric composition and processes. We present a new multifluid MHD version of the Block-Adaptive-Tree Solar wind Roe-type Upwind Scheme model of the magnetosphere to track the fate and consequences of ionospheric outflow. The multifluid MHD equations are presented as are the novel techniques for overcoming the formidable challenges associated with solving them. Our new model is then applied to the May 4, 1998 and March 31, 2001 geomagnetic storms. The results are juxtaposed with traditional single-fluid MHD and multispecies MHD simulations from a previous study, thereby allowing us to assess the benefits of using a more complex model with additional physics. We find that our multifluid MHD model (with outflow) gives comparable results to the multispecies MHD model (with outflow), including a more strongly negative Dst, reduced CPCP, and a drastically improved magnetic field at geosynchronous orbit, as compared to single-fluid MHD with no outflow. Significant differences in composition and magnetic field are found between the multispecies and multifluid approach further away from the Earth. We further demonstrate the ability to explore pressure and bulk velocity differences between H + and O + , which is not possible when utilizing the other techniques considered.

The search for life in the Universe is a fundamental problem of astrobiology and modern science. The current progress in the detection of terrestrial-type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favourable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of global (astrospheric), and local (atmospheric and surface) environments of exoplanets in the habitable zones (HZs) around G-K-M dwarf stars including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favourable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro)physical, chemical and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the HZ to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field in light of presentations and discussions during the NASA Nexus for Exoplanetary System Science funded workshop ‘Exoplanetary Space Weather, Climate and Habitability’ and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.