The solar wind (SW) and the extreme ultraviolet (EUV) radiation modulate fluxes of interstellar and heliospheric particles inside the heliosphere both in time and in space. Understanding this modulation is necessary to correctly interpret measurements of particles of interstellar origin inside the heliosphere. We present a revision of heliospheric ionization rates and provide the Sun-Heliosphere Observation-based Ionization Rates (SHOIR) model based on the currently available data. We calculate the total ionization rates using revised SW and solar EUV data. We study the in-ecliptic variation of the SW parameters, the latitudinal structure of the SW speed and density, and the reconstruction of the photoionization rates. The revision affects the most the SW out of the ecliptic plane during solar maximum and the estimation of the photoionization rates. The latter, due to a change of the reference data. The revised polar SW is slower and denser during the solar maximum of SC 24. The currently estimated total ionization rates are higher than the previous ones for H, O, and Ne, and lower for He. The changes for the in-ecliptic total ionization rates are less than 10% for H and He, up to 20% for O, and up to 35% for Ne. Additionally, the changes are not constant in time and vary as a function of time and latitude.

We present a modification of a model of solar cycle evolution of the solar Lyman-alpha line profile, along with a sensitivity study of interstellar neutral H hydrogen to uncertainties in radiation pressure level. The line profile model, originally developed by Kowalska-Leszczynska et al. 2018a, is parametrized by the composite solar Lyman-alpha flux, which recently was revised Machol et al. 2019. We present modified parameters of the previously-developed model of solar radiation pressure for neutral hydrogen and deuterium atoms in the heliosphere. The mathematical function used in the model, as well as the fitting procedure, remain unchanged. We show selected effects of the model modification on ISN H properties in the heliosphere and we discuss the sensitivity of these quantities to uncertainties in the calibration of the composite Lyman-alpha series.

The interstellar neutral (ISN) gas enters the heliosphere and is detected at a few au from the Sun, as demonstrated by Ulysses and the nterstellar Boundary Explorer (IBEX) missions. Ulysses observed ISN gas from different vantage points in a polar orbit from 1994 to 2007, while IBEX has been observing in an Earth orbit in a fixed direction relative to the Sun from 2009. McComas et al. 2018 reported about an IMAP-Lo detector on board the Interstellar Mapping and Acceleration Probe (IMAP), with an ability to track the ISN flux in the sky. We present observation geometries for ISN gas for a detector with the capability to adjust the boresight direction along the Earth orbit over a year within a multichoice ISN observation scheme. We study science opportunities from the observations as a function of time during a year and the phase of solar activity. We identify observation geometries and determine the observation seasons separately for various ISN species and populations. We find that using an adjustable viewing direction allows for ISN gas observations in the upwind hemisphere, where the signal is not distorted by gravitational focusing, in addition to the viewing of ISN species throughout the entire year. Moreover, we demonstrate that with appropriately adjusted observation geometries, primary and secondary populations can be fully separated. Additionally, we show that atoms of ISN gas on indirect trajectories are accessible for detection, and we present their impact on the study of the ionization rates for ISN species.

Energetic neutral atoms (ENA) are an important tool for investigating the structure of the heliosphere. Recently, it was observed that fluxes of ENAs (with energy $\le$ 55 keV) coming from the upwind and downwind regions of the heliosphere are similar in strength. This led the authors of these observations to hypothesize that the heliosphere is bubble-like rather than comet-like, meaning that it has no extended tail. We investigate the directional distribution of the ENA flux for a wide energy range (3--88 keV) including the observations from IBEX (Interstellar Boundary Explorer), INCA (Ion and Neutral Camera, on board Cassini), and HSTOF (High energy Suprathermal Time Of Flight sensor, on board SOHO, Solar and Heliospheric Observatory). An essential element is the model of pickup ion acceleration at the termination shock (TS) proposed by Zank. We use state of the art models of the global heliosphere, interstellar neutral gas density, and pickup ion distributions. The results, based on the "comet-like" model of the heliosphere, are close in flux magnitude to ENA observations by IBEX, HSTOF and partly by INCA (except for the 5.2-13.5 keV energy channel). We find that the ENA flux from the tail dominates at high energy (in agreement with HSTOF, but not INCA). At low energy, our comet-like model produces the similar strengths of the ENA fluxes from the upwind and downwind directions, which, therefore, removes this as a compelling argument for a bubble-like heliosphere.

We study the distribution of the interstellar neutral (ISN) gas density and the pick-up ion (PUI) density of hydrogen, helium, neon, and oxygen in the heliosphere for heliocentric distances from inside 1 au up to the solar wind termination shock (TS), both in and out of the ecliptic plane. We discuss similarities and differences in the large-scale structures of the ISN gas and PUIs formed in the heliosphere between various species. We discuss the distribution of ISN gas and PUI densities for two extreme phases of the solar activity cycle, it is the solar minimum and the solar maximum. We identify the location of the ISN gas density cavity of various species. We study the relative abundance ratios of Ne/O, H/He, Ne/He, and O/He for ISN gas and PUIs densities and their variation with location in the heliosphere. We also discuss the modulation of relative abundance ratios of ISN gas and PUIs along the TS. We conclude that the preferable locations for detection of He$^+$ and Ne$^+$ PUIs are in the downwind hemisphere within 1~au, whereas for H$^+$ and O$^+$ PUIs the preferable locations for detection are for distances from Jupiter to Pluto orbits.

Following the derivation of a more accurate model of the evolution of the solar Lyman-$\alpha$ line with the changing solar activity by Kowalska-leszczynska et al. 2018 (IKL18) than the formerly used model by Tarnopolski et al. 2009 (ST09), we investigate potential consequences that adoption of the resulting refined model of radiation pressure has for the model distribution of interstellar neutral (ISN) H in the inner heliosphere and on the interpretation of selected observations. We simulated the ISN H densities using the two alternative radiation pressure models and identical models of all other factors affecting the ISN H distribution. We found that during most of the solar cycle, the IKL18 model predicts larger densities of ISN H and PUIs than ST09 in the inner heliosphere, especially in the downwind hemisphere. However, the density of ISN H at the termination shock estimated by Bzowski et al. 2008 obtained using ST09 does not need revision, and the detection of ISN D by IBEX is supported. However, we point out the existence of a considerable absorption of a portion of the solar Lyman-$\alpha$ spectral flux inside the heliosphere. Therefore, the model of radiation pressure for ISN H is still likely to need revision, and hence the available models of ISN H are not self-consistent.

Recent studies of interstellar neutral (ISN) hydrogen observed by the Interstellar Boundary Explorer (IBEX) suggested that the present understanding of the radiation pressure acting on hydrogen atoms in the heliosphere should be revised. There is a significant discrepancy between theoretical predictions of the ISN H signal using the currently used model of the solar Lyman-alpha profile by Tarnopolski et al. 2009 (TB09) and the signal due to ISN H observed by IBEX-Lo. We developed a new model of evolution of the solar Lyman-alpha profile that takes into account all available observations of the full-disk solar Lyman-alpha profiles from SUMER/SOHO, provided by Lemaire et al. 2015 (L15), covering practically the entire 23rd solar cycle. The model has three components that reproduce different features of the profile. The main shape of the emission line that is produced in the chromosphere is modeled by the kappa function; the central reversal due to absorption in the transition region is modeled by the Gauss function; the spectral background is represented by the linear function. The coefficients of all those components are linear functions of the line-integrated full-disk Lyman-alpha irradiance, which is the only free parameter of the model. The new model features potentially important differences in comparison with the model by TB09, which was based on a limited set of observations. This change in the understanding of radiation pressure, especially during low solar activity, may significantly affect the interstellar H and D distributions in the inner heliosphere and their derivative populations.

Observations of energetic neutral atoms (ENAs) allow for remote studies of the plasma condition in the heliosphere and the neighboring local interstellar medium. The first IBEX results revealed an arc-like enhancement of the ENA intensity in the sky, known as the ribbon. The ribbon was not expected from the heliospheric models prior to the IBEX launch. One of the proposed explanations of the ribbon is the mechanism of the secondary ENA emission. The ribbon reveals energy-dependent structure in the relative intensity along its circumference and in the position. Namely, the ribbon geometric center varies systematically by about 10$^\circ$ in the energy range 0.7-4.3 keV. Here, we show by analytic modeling that this effect is a consequence of the helio-latitudinal structure of the solar wind reflected in the secondary ENAs. Along with a recently measured distance to the ribbon source just beyond the heliopause, our findings support the connection of the ribbon with the local interstellar magnetic field by the mechanism of the secondary ENA emission. However, the magnitude of the center shift in the highest IBEX energy channel is much larger in the observations than expected from the modeling. This may be due to another, not currently recognized, process of ENA generation.

Direct-sampling observations of interstellar neutral (ISN) He by Interstellar Boundary Explorer (IBEX) provide valuable insight into the physical state of and processes operating in the interstellar medium ahead of the heliosphere. The ISN He atom signals are observed at the four lowest ESA steps of the IBEX-Lo sensor. The observed signal is a mixture of the primary and secondary components of ISN He and H. Previously, only data from one of the ESA steps have been used. Here, we extended the analysis to data collected in the three lowest ESA steps with the strongest ISN He signal, for the observation seasons 2009-2015. The instrument sensitivity is modeled as a linear function of the atom impact speed onto the sensor's conversion surface separately for each ESA step of the instrument. We found that the sensitivity increases from lower to higher ESA steps, but within each of the ESA steps it is a decreasing function of the atom impact speed. This result may be influenced by the hydrogen contribution, which was not included in the adopted model, but seems to exist in the signal. We conclude that the currently accepted temperature of ISN He and velocity of the Sun through the interstellar medium do not need a revision, and we sketch a plan of further data analysis aiming at investigating ISN H and a better understanding of the population of ISN He originating in the outer heliosheath.

The Interstellar Boundary Explorer (IBEX) has now operated in space for 7 years and returned nearly continuous observations that have led to scientific discoveries and reshaped our entire understanding of the outer heliosphere and its interaction with the local interstellar medium. Here we extend prior work, adding the 2014-2015 data for the first time, and examine, validate, initially analyze, and provide a complete 7-year set of Energetic Neutral Atom (ENA) observations from ~0.1 to 6 keV. The data, maps, and documentation provided here represent the 10th major release of IBEX data and include improvements to various prior corrections to provide the citable reference for the current version of IBEX data. We are now able to study time variations in the outer heliosphere and interstellar interaction over more than half a solar cycle. We find that the Ribbon has evolved differently than the globally distributed flux (GDF), with a leveling off and partial recovery of ENAs from the GDF, owing to solar wind output flattening and recovery. The Ribbon has now also lost its latitudinal ordering, which reflects the breakdown of solar minimum solar wind conditions and exhibits a greater time delay than for the surrounding GDF. Together, the IBEX observations strongly support a secondary ENA source for the Ribbon, and we suggest that this be adopted as the nominal explanation of the Ribbon going forward.

We simulated the modulation of the interstellar neutral (ISN) He, Ne, and O density and pick-up ion (PUI) production rate and count rate along the Earth's orbit over the solar cycle from 2002 to 2013 to verify if solar cycle-related effects may modify the inferred ecliptic longitude of the ISN inflow direction. We adopted the classical PUI model with isotropic distribution function and adiabatic cooling, modified by time- and heliolatitude-dependent ionization rates and non-zero injection speed of PUIs. We found that the ionization losses have a noticeable effect on the derivation of the ISN inflow longitude based on the Gaussian fit to the crescent and cone peak locations. We conclude that the non-zero radial velocity of the ISN flow and the energy range of the PUI distribution function that is accumulated are of importance for a precise reproduction of the PUI count rate along the Earth orbit. However, the temporal and latitudinal variations of the ionization in the heliosphere, and particularly their variation on the solar cycle time-scale, may significantly modify the shape of PUI cone and crescent and also their peak positions from year to year and thus bias by a few degrees the derived longitude of the ISN gas inflow direction.

The modeling of the heliosphere requires continuous three-dimensional solar wind data. The in-situ out-of-ecliptic measurements are very rare, so that other methods of solar wind detection are needed. We use the remote-sensing data of the solar wind speed from observations of interplanetary scintillation (IPS) to reconstruct spatial and temporal structures of the solar wind proton speed from 1985 to 2013. We developed a method of filling the data gaps in the IPS observations to obtain continuous and homogeneous solar wind speed records. We also present a method to retrieve the solar wind density from the solar wind speed, utilizing the invariance of the solar wind dynamic pressure and energy flux with latitude. To construct the synoptic maps of the solar wind speed we use the decomposition into spherical harmonics of each of the Carrington rotation map. To fill the gaps in time we apply the singular spectrum analysis to the time series of the coefficients of spherical harmonics. We obtained helio-latitudinal profiles of the solar wind proton speed and density over almost three recent solar cycles. The accuracy in the reconstruction is, due to computational limitations, about 20%. The proposed methods allow us to improve the spatial and temporal resolution of the model of the solar wind parameters presented in our previous paper (Sok\'o{\l} et al. 2013) and give a better insight into the time variations of the solar wind structure. Additionally, the solar wind density is reconstructed more accurately and it fits better to the in-situ measurements from Ulysses.

This paper is one of three companion papers presenting the results of our in-depth analysis of the interstellar neutral helium (ISN He) observations carried out using the IBEX-Lo during the first six Interstellar Boundary Explorer (IBEX) observation seasons. We derive corrections for losses due to the limited throughput of the interface buffer and determine the IBEX spin-axis pointing. We develop an uncertainty system for the data, taking into account the resulting correlations between the data points. This system includes uncertainties due to Poisson statistics, background, spin-axis determination, systematic deviation of the boresight from the prescribed position, correction for the interface buffer losses, and the expected Warm Breeze (WB) signal. Subsequently, we analyze the data from 2009 to examine the role of various components of the uncertainty system. We show that the ISN He flow parameters are in good agreement with the values obtained by the original analysis. We identify the WB as the principal contributor to the global $\chi^2$ values in previous analyses. Other uncertainties have a much milder role and their contributions are comparable to each other. The application of this uncertainty system reduced the minimum $\chi^2$ value 4-fold. The obtained $\chi^2$ value, still exceeding the expected value, suggests that either the uncertainty system may still be incomplete or the adopted physical model lacks a potentially important element, which is likely an imperfect determination of the WB parameters. The derived corrections and uncertainty system are used in the accompanying paper by Bzowski et al. in an analysis of the data from six seasons.

Neutral interstellar helium has been observed by the Interstellar Boundary Explorer (IBEX) since 2009 with a signal-to-noise ratio well above 1000. Because of the geometry of the observations, the signal observed from January to March each year is the easiest to identify. However, as we show via simulations, the portion of the signal in the range of intensities from 10^{-3} to 10^{-2} of the peak value, previously mostly left out from the analysis, may bring important information about the details of the distribution function of interstellar He gas in front of the heliosphere. In particular, these observations may inform us about possible departures of the parent interstellar He population from equilibrium. We compare the expected distribution of the signal for the canonical assumption of a single Maxwell-Boltzmann population with the distributions for a superposition of the Maxwell-Boltzmann primary population and the recently discovered Warm Breeze, and for a single primary population given by a kappa function. We identify the regions on the sky where the differences between those cases are expected to be the most visible against the background. We discuss the diagnostic potential of the fall peak of the interstellar signal, reduced by a factor of 50 due to the Compton-Getting effect but still above the detection limit of IBEX. We point out the strong energy dependence of the fall signal and suggest that searching for this signal in the data could bring an independent assessment of the low-energy measurement threshold of the IBEX-Lo sensor.

We have developed a refined and optimized version of the Warsaw Test Particle Model of interstellar neutral gas in the heliosphere, specially tailored for analysis of IBEX-Lo observations. The former version of the model was used in the analysis of neutral He observed by IBEX that resulted in an unexpected conclusion that the interstellar neutral He flow vector was different than previously thought and that a new population of neutral He, dubbed the Warm Breeze, exists in the heliosphere. It was also used in the reanalysis of Ulysses observations that confirmed the original findings on the flow vector, but suggested a significantly higher temperature. The present version model has two strains targeted for different applications, based on an identical paradigm, but differing in the implementation and in the treatment of ionization losses. We present the model in detail and discuss numerous effects related to the measurement process that potentially modify the resulting flux of ISN He observed by IBEX, and identify those of them that should not be omitted in the simulations to avoid biasing the results. This paper is part of a coordinated series of papers presenting the current state of analysis of IBEX-Lo observations of ISN He. Details of the analysis method are presented by Swaczyna et al. 2015, and results of the analysis are presented by Bzowski et al. 2015.

The NIS He gas has been observed at a few AU to the Sun almost from the beginning of the space age. To model its flow an estimate of the loss rates due to ionization by solar extreme-ultraviolet (EUV) flux is needed. The EUV irradiance has been measured directly from mid 1990-ties, but with high temporal and spectral resolution only from 2002. Beforehand only EUV proxies are available. A new method of reconstruction of the Carrington rotation averaged photoionization rates for neutral interstellar helium (NIS He) in the ecliptic at 1 AU to the Sun before 2002 is presented. We investigate the relation between the solar rotation averaged time series of the ionization rates for NIS He at 1 AU derived from TIMED measurements of EUV irradiance and the solar 10.7 cm flux (F10.7) only. We perform a weighted iterative fit of a nonlinear model to data split into sectors. The obtained formula allows to reconstruct the solar rotation averages of photoionization rates for He between ~1947 and 2002 with an uncertainty ranging from less than 10% during solar minimum up to 20% for solar maximum.

We investigate the signals from neutral He atoms observed from Earth orbit in 2010 by IBEX. The full He signal observed during the 2010 observation season can be explained as a superposition of pristine neutral interstellar He gas and an additional population of neutral He that we call the Warm Breeze. The Warm Breeze is approximately two-fold slower and 2.5 times warmer than the primary interstellar He population, and its density in front of the heliosphere is ~7% that of the neutral interstellar helium. The inflow direction of the Warm Breeze differs by ~19deg from the inflow direction of interstellar gas. The Warm Breeze seems a long-term feature of the heliospheric environment. It has not been detected earlier because it is strongly ionized inside the heliosphere, which brings it below the threshold of detection via pickup ion and heliospheric backscatter glow observations, as well as by the direct sampling of GAS/Ulysses. Possible sources for the Warm Breeze include (1) the secondary population of interstellar helium, created via charge exchange and perhaps elastic scattering of neutral interstellar He atoms on interstellar He+ ions in the outer heliosheath, or (2) a gust of interstellar He originating from a hypothetic wave train in the Local Interstellar Cloud. A secondary population is expected from models, but the characteristics of the Warm Breeze do not fully conform to modeling results. If, nevertheless, this is the explanation, IBEX-Lo observations of the Warm Breeze provide key insights into the physical state of plasma in the outer heliosheath. If the second hypothesis is true, the source is likely to be located within a few thousand of AU from the Sun, which is the propagation range of possible gusts of interstellar neutral helium with the Warm Breeze characteristics against dissipation via elastic scattering in the Local Cloud.

Direct sampling of neutral interstellar (NIS) atoms by the Interstellar Boundary Explorer (IBEX) can potentially provide a complementary method for studying element abundances in the Local Interstellar Cloud and processes in the heliosphere interface.}{We set the stage for abundance-aimed in-depth analysis of measurements of NIS He, Ne, and O by IBEX and determine systematic differences between abundances derived from various calculation methods and their uncertainties.}{Using a model of ionization rates of the NIS species in the heliosphere, based on independent measurements of the solar wind and solar EUV radiation, we develop a time-dependent method of calculating the survival probabilities of NIS atoms from the termination shock (TS) of the solar wind to IBEX. With them, we calculate densities of these species along the Earth's orbit and simulate the fluxes of NIS species as observed by IBEX. We study pairwise ratios of survival probabilities, densities and fluxes of NIS species at IBEX to calculate correction factors for inferring the abundances at TS.}{The analytic method to calculate the survival probabilities gives acceptable results only for He and Ne during low solar activity. For the remaining portions of the solar cycle, and at all times for O, a fully time dependent model should be used. Electron impact ionization is surprisingly important for NIS O. Interpreting the IBEX observations using the time dependent model yields the LIC Ne/O abundance of $0.16\pm40%$. The uncertainty is mostly due to uncertainties in the ionization rates and in the NIS gas flow vector.}{The Ne/He, O/He and Ne/O ratios for survival probabilities, local densities, and fluxes scaled to TS systematically differ and thus an analysis based only on survival probabilities or densities is not recommended, except the Ne/O abundance for observations at low solar activity.

The abundance of deuterium in the interstellar gas in front of the Sun gives insight into the processes of filtration of neutral interstellar species through the heliospheric interface and potentially into the chemical evolution of the Galactic gas. We investigate the possibility of detection of neutral interstellar deuterium at 1 AU from the Sun by direct sampling by the Interstellar Boundary Explorer (IBEX). We simulate the flux of neutral interstellar D at IBEX for the actual measurement conditions. We assess the number of interstellar D atom counts expected during the first three years of IBEX operation. We also simulate observations expected during an epoch of high solar activity. In addition, we calculate the expected counts of D atoms from the thin terrestrial water layer, sputtered from the IBEX-Lo conversion surface by neutral interstellar He atoms. Most D counts registered by IBEX-Lo are expected to originate from the water layer, exceeding the interstellar signal by 2 orders of magnitude. However, the sputtering should stop once the Earth leaves the portion of orbit traversed by interstellar He atoms. We identify seasons during the year when mostly the genuine interstellar D atoms are expected in the signal. During the first 3 years of IBEX operations about 2 detectable interstellar D atoms are expected. This number is comparable with the expected number of sputtered D atoms registered during the same time intervals. The most favorable conditions for the detection occur during low solar activity, in an interval including March and April each year. The detection chances could be improved by extending the instrument duty cycle, e.g., by making observations in the special deuterium mode of IBEX-Lo.

NASA's Interstellar Boundary Explorer (IBEX) mission has recently opened a new window on the interstellar medium (ISM) by imaging neutral atoms. One "bright" feature in the sky is the interstellar wind flowing into the solar system. Composed of remnants of stellar explosions as well as primordial gas and plasma, the ISM is by no means uniform. The interaction of the local ISM with the solar wind shapes our heliospheric environment with hydrogen being the dominant component of the very local ISM. In this paper, we report on direct sampling of the neutral hydrogen of the local ISM over four years of IBEX observations. The hydrogen wind observed at 1 AU has decreased and nearly disappeared as the solar activity has increased over the last four years; the signal at 1 AU has dropped off in 2012 by a factor of ~8 to near background levels. The longitudinal offset has also increased with time presumably due to greater radiation pressure deflecting the interstellar wind. We present longitudinal and latitudinal arrival direction measurements of the bulk flow as measured over four years beginning at near solar minimum conditions. The H distribution we observe at 1 AU is expected to be different from that outside the heliopause due to ionization, photon pressure, gravity, and filtration by interactions with heliospheric plasma populations. These observations provide an important benchmark for modeling of the global heliospheric interaction. Based on these observations we suggest a further course of scientific action to observe neutral hydrogen over a full solar cycle with IBEX.