With upcoming missions such as the James Webb Space Telescope (JWST), the European Extremely Large Telescope (ELT), and the Atmospheric Remote-sensing Infrared Exoplanet Large-survey (ARIEL), we soon will be on the verge of detecting and characterizing Earth-like exoplanetary atmospheres for the first time. These planets are most likely to be found around smaller and cooler K- and M-type stars. However, recent observations showed that their radiation environment might be much harsher than that of the Sun. Thus, the exoplanets are most likely exposed to an enhanced stellar radiation environment, which could affect their habitability, for example, in the form of a hazardous flux of energetic particles. Knowing the stellar radiation field, and being able to model the radiation exposure on the surface of a planet is crucial to assess its habitability. In this study, we present 3D magnetohydrodynamic (MHD)-based model efforts investigating M-stars, focusing on V374 Peg, Proxima Centauri, and LHS 1140, chosen because of their diverse astrospheric quantities. We show that V374 Peg has a much larger astrosphere (ASP) than our Sun, while Proxima Centauri and LHS 1140 most likely have ASPs comparable or even much smaller than the heliosphere, respectively. Based on a 1D transport model, for the first time, we provide numerical estimates of the modulation of Galactic cosmic rays (GCRs) within the three ASPs. We show that the impact of GCRs on the Earth-like exoplanets Proxima Centauri b and LHS 1140 b cannot be neglected in the context of exoplanetary habitability.

Jovian electrons serve as an important test-particle distribution in the inner heliosphere and have been used extensively in the past to study the (diffusive) transport of cosmic rays in the inner heliosphere. With new limits on the Jovian source function (i.e. the particle intensity just outside the Jovian magnetosphere), and a new set of in-situ observations at 1 AU for both cases of good and poor magnetic connection between the source and observer, we revisit some of these earlier simulations. We aim to find the optimal numerical set-up that can be used to simulate the propagation of 6 MeV Jovian electrons in the inner heliosphere. Using such a set-up, we further aim to study the residence (propagation) times of these particles for different levels of magnetic connection between Jupiter and an observer at Earth (1 AU). Using an advanced Jovian electron propagation model based on the stochastic differential equation (SDE) approach, we calculate the Jovian electron intensity for different model parameters. A comparison with observations leads to an optimal numerical set-up, which is then used to calculate the so-called residence (propagation) times of these particles. Comparing to in-situ observations, we are able to derive transport parameters that are appropriate to study the propagation of 6 MeV Jovian electrons in the inner heliosphere. Moreover, using these values, we show that the method of calculating the residence time applied in former literature is not suited to being interpreted as the propagation time of physical particles. This is due to an incorrect weighting of the probability distribution. We propose and apply a new method, where the results from each pseudo-particle are weighted by its resulting phase-space density (i.e. the number of physical particles that it represents). Thereby we obtain more reliable estimates for the propagation time.

The atmospheres of small, potentially rocky exoplanets are expected to cover a diverse range in composition and mass. Studying such objects therefore requires flexible and wide-ranging modeling capabilities. We present in this work the essential development steps that lead to our flexible radiative transfer module, REDFOX, and validate REDFOX for the Solar system planets Earth, Venus and Mars, as well as for steam atmospheres. REDFOX is a k-distribution model using the correlated-k approach with random overlap method for the calculation of opacities used in the $\delta$-two-stream approximation for radiative transfer. Opacity contributions from Rayleigh scattering, UV / visible cross sections and continua can be added selectively. With the improved capabilities of our new model, we calculate various atmospheric scenarios for K2-18b, a super-Earth / sub-Neptune with $\sim$8 M$_\oplus$ orbiting in the temperate zone around an M-star, with recently observed H$_2$O spectral features in the infrared. We model Earth-like, Venus-like, as well as H$_2$-He primary atmospheres of different Solar metallicity and show resulting climates and spectral characteristics, compared to observed data. Our results suggest that K2-18b has an H$_2$-He atmosphere with limited amounts of H$_2$O and CH$_4$. Results do not support the possibility of K2-18b having a water reservoir directly exposed to the atmosphere, which would reduce atmospheric scale heights, hence too the amplitudes of spectral features inconsistent with the observations. We also performed tests for H$_2$-He atmospheres up to 50 times Solar metallicity, all compatible with the observations.

Due to its Earth-like minimum mass of 1.27 M$_{\text{E}}$ and its close proximity to our Solar system, Proxima Centauri b is one of the most interesting exoplanets for habitability studies. Its host star, Proxima Centauri, is however a strongly flaring star, which is expected to provide a very hostile environment for potentially habitable planets. We perform a habitability study of Proxima Centauri b assuming an Earth-like atmosphere under high stellar particle bombardment, with a focus on spectral transmission features. We employ our extensive model suite calculating energy spectra of stellar particles, their journey through the planetary magnetosphere, ionosphere, and atmosphere, ultimately providing planetary climate and spectral characteristics, as outlined in Herbst et al. (2019). Our results suggest that together with the incident stellar energy flux, high particle influxes can lead to efficient heating of the planet well into temperate climates, by limiting CH$_4$ amounts, which would otherwise run into anti-greenhouse for such planets around M-stars. We identify some key spectral features relevant for future spectral observations: First, NO$_2$ becomes the major absorber in the visible, which greatly impacts the Rayleigh slope. Second, H$_2$O features can be masked by CH$_4$ (near infra-red) and CO$_2$ (mid to far infra-red), making them non-detectable in transmission. Third, O$_3$ is destroyed and instead HNO$_3$ features become clearly visible in the mid to far infra-red. Lastly, assuming a few percent of CO$_2$ in the atmosphere, CO$_2$ absorption at 5.3 $\mu$m becomes significant (for flare and non-flare cases), strongly overlapping with a flare related NO feature in Earth's atmosphere.

The interpretation of recent observations of bow shocks around O-stars and the creation of corresponding models require a detailed understanding of the associated (magneto-)hydrodynamic structures. We base our study on three-dimensional numerical magneto-)hydrodynamical models, which are analyzed using the dynamically relevant parameters, in particular, the (magneto)sonic Mach numbers. The analytic Rankine-Hugoniot relation for HD and MHD are compared with those obtained by the numerical model. In that context we also show that the only distance which can be approximately determined is that of the termination shock, if it is a hydrodynamical shock. For MHD shocks the stagnation point does not, in general, lie on the inflow line, which is the line parallel to the inflow vector and passing through the star. Thus an estimate via the Bernoulli equation as in the HD case is, in general, not possible. We also show that in O-star astrospheres, distinct regions exist in which the fast, slow, Alfv\'enic, and sonic Mach numbers become lower than one, implying sub-slow magnetosonic as well as sub-fast and sub-sonic flows. Nevertheless, the analytic MHD Rankine Hugoniot relations can be used for further studies of turbulence and cosmic ray modulation.

The Atmospheric Radiation Interaction Simulator (AtRIS) was used to model the altitude-dependent Venusian absorbed dose and the Venusian dose equivalent. For the first time, we modeled the dose rates for different shape-, size-, and composition-mimicking detectors (phantoms): a CO$_2$-based phantom, a water-based microbial cell, and a phantom mimicking human tissue. Based on a new model approach, we give a reliable estimate of the altitude-dependent Venusian radiation dose in water-based microorganisms here for the first time. These microorganisms are representative of known terrestrial life. We also present a detailed analysis of the influence of the strongest ground-level enhancements measured at the Earth's surface, and of the impact of two historic extreme solar events on the Venusian radiation dose. Our study shows that because a phantom based on Venusian air was used, and because furthermore, the quality factors of different radiation types were not taken into account, previous model efforts have underestimated the radiation hazard for any putative Venusian cloud-based life by up to a factor of five. However, because we furthermore show that even the strongest events would not have had a hazardous effect on putative microorganisms within the potentially habitable zone (51 km - 62 km), these differences may play only a minor role.

The first opportunity to detect indications for life outside the Solar System may be provided already within the next decade with upcoming missions such as the James Webb Space Telescope (JWST), the European Extremely Large Telescope (E-ELT) and/or the Atmospheric Remote-sensing Infrared Exoplanet Large-survey (ARIEL) mission, searching for atmospheric biosignatures on planets in the habitable zone of cool K- and M-stars. Nevertheless, their harsh stellar radiation and particle environment could lead to photochemical loss of atmospheric biosignatures. We aim to study the influence of cosmic rays on exoplanetary atmospheric biosignatures and the radiation environment considering feedbacks between energetic particle precipitation, climate, atmospheric ionization, neutral and ion chemistry, and secondary particle generation. We describe newly-combined state-of-the-art modeling tools to study the impact of the radiation and particle environment on atmospheric particle interaction, the influence on the atmospheric chemistry, and the climate-chemistry coupling in a self-consistent model suite. To this end, models like the Atmospheric Radiation Interaction Simulator (AtRIS), the Exoplanetary Terrestrial Ion Chemistry model (ExoTIC), and the updated coupled climate-chemistry model are combined. Amongst others, we model the atmospheric response during quiescent solar periods and during a strong solar energetic particle event as well as the scenario-dependent terrestrial transit spectra, as seen by the NIR-Spec infrared spectrometer onboard the JWST. We find that the comparatively weak solar event drastically increases the spectral signal of HNO$_3$, while significantly suppressing the spectral feature of ozone. Because of the slow recovery after such events, the latter indicates that ozone might not be a good biomarker for planets orbiting stars with high flaring rates.

On July 25 2017 a multi-step Forbush decrease (FD) with the remarkable total amplitude of more than 15\% was observed by MSL/RAD at Mars. We find that these particle signatures are related to very pronounced plasma and magnetic field signatures detected in situ by STEREO-A on July 24 2017, with a higher than average total magnetic field strength reaching more than 60 nT. In the observed time period STEREO-A was at a relatively small longitudinal separation (46 degrees) to Mars and both were located at the back side of the Sun as viewed from Earth. We analyse a number of multi-spacecraft and multi-instrument (both in situ and remote-sensing) observations, and employ modelling to understand these signatures. We find that the solar sources are two CMEs which erupted on July 23 2017 from the same source region on the back side of the Sun as viewed from Earth. Moreover, we find that the two CMEs interact non-uniformly, inhibiting the expansion of one of the CMEs in STEREO-A direction, whereas allowing it to expand more freely in the Mars direction. The interaction of the two CMEs with the ambient solar wind adds up to the complexity of the event, resulting in a long, sub-structured interplanetary disturbance at Mars, where different sub-structures correspond to different steps of the FD, adding-up to a globally large-amplitude FD.

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 favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs 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 favorable 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 habitable zone 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 and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

On their way through the heliosphere, Galactic Cosmic Rays (GCRs) are modulated by various effects before they can be detected at Earth. This process can be described by the Parker equation, which calculates the phase space distribution of GCRs depending on the main modulation processes: convection, drifts, diffusion and adiabatic energy changes. A first order approximation of this equation is the force field approach, reducing it to a one-parameter dependency, the solar modulation potential $\phi$. Utilizing this approach, it is possible to reconstruct $\phi$ from ground based and spacecraft measurements. However, it has been shown previously that $\phi$ depends not only on the Local Interstellar Spectrum (LIS) but also on the energy range of interest. We have investigated this energy dependence further, using published proton intensity spectra obtained by PAMELA as well as heavier nuclei measurements from IMP-8 and ACE/CRIS. Our results show severe limitations at lower energies including a strong dependence on the solar magnetic epoch. Based on these findings, we will outline a new tool to describe GCR proton spectra in the energy range from a few hundred MeV to tens of GeV over the last solar cycles. In order to show the importance of our modification, we calculate the global production rates of the cosmogenic radionuclide $^{10}$Be which is a proxy for the solar activity ranging back thousands of years.