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.

Coronal mass ejections (CMEs) cause disturbances in the environment of the Earth when they arrive at the Earth. However, the prediction of the arrival of CMEs still remains a challenge. We have developed an interplanetary scintillation (IPS) estimation system based on a global magnetohydrodynamic (MHD) simulation of the inner heliosphere to predict the arrival time of CMEs. In this system, the initial speed of a CME is roughly derived from white light coronagraph observations. Then, the propagation of the CME is calculated by a global MHD simulation. The IPS response is estimated by the three-dimensional density distribution of the inner heliosphere derived from the MHD simulation. The simulated IPS response is compared with the actual IPS observations made by the Institute for Space-Earth Environmental Research, Nagoya University, and shows good agreement with that observed. We demonstrated how the simulation system works using a halo CME event generated by a X9.3 flare observed on September 5, 2017. We find that the CME simulation that best estimates the IPS observation can more accurately predict the time of arrival of the CME at the Earth. These results suggest that the accuracy of the CME arrival time can be improved if our current MHD simulations include IPS data.

We comparatively studied the long-term variation (1992-2017) in polar brightening observed with the Nobeyama Radioheliograph, the polar solar wind velocity with interplanetary scintillation observations at the Institute for Space-Earth Environmental Research, and the coronal hole distribution computed by potential field calculations of the solar corona using synoptic magnetogram data obtained at Kitt Peak National Solar Observatory. First, by comparing the solar wind velocity (V) and the brightness temperature (T_b) in the polar region, we found good correlation coefficients (CCs) between V and T_b in the polar regions, CC = 0.91 (0.83) for the northern (southern) polar region, and we obtained the V-T_b relationship as V =12.6 (T_b-10,667)^{1/2}+432. We also confirmed that the CC of V-T_b is higher than those of V-B and V-B/f, where B and f are the polar magnetic field strength and magnetic flux expansion rate, respectively. These results indicate that T_b is a more direct parameter than B or B/f for expressing solar wind velocity. Next, we analyzed the long-term variation of the polar brightening and its relation to the area of the polar coronal hole (A). As a result, we found that the polar brightening matches the probability distribution of the predicted coronal hole and that the CC between T_b and A is remarkably high, CC = 0.97. This result indicates that the polar brightening is strongly coupled to the size of the polar coronal hole. Therefore, the reasonable correlation of V-T_b is explained by V-A. In addition, by considering the anti-correlation between A and f found in a previous study, we suggest that the V-T_b relationship is another expression of the Wang-Sheeley relationship (V-1/f) in the polar regions.

We analyze the short cosmic ray intensity increase ("cosmic ray burst": CRB) on June 22, 2015 utilizing a global network of muon detectors and derive the global anisotropy of cosmic ray intensity and the density (i.e. the omnidirectional intensity) with 10-minute time resolution. We find that the CRB was caused by a local density maximum and an enhanced anisotropy of cosmic rays both of which appeared in association with Earth's crossing of the heliospheric current sheet (HCS). This enhanced anisotropy was normal to the HCS and consistent with a diamagnetic drift arising from the spatial gradient of cosmic ray density, which indicates that cosmic rays were drifting along the HCS from the north of Earth. We also find a significant anisotropy along the HCS, lasting a few hours after the HCS crossing, indicating that cosmic rays penetrated into the inner heliosphere along the HCS. Based on the latest geomagnetic field model, we quantitatively evaluate the reduction of the geomagnetic cut-off rigidity and the variation of the asymptotic viewing direction of cosmic rays due to a major geomagnetic storm which occurred during the CRB and conclude that the CRB is not caused by the geomagnetic storm, but by a rapid change in the cosmic ray anisotropy and density outside the magnetosphere.

We present the first dedicated observations of Interplanetary Scintillation (IPS) with the Murchison Widefield Array (MWA). We have developed a synthesis imaging technique, tailored to the properties of modern "large-N" low-frequency radio telescopes. This allows us to image the variability on IPS timescales across 900 square degrees simultaneously. We show that for our observations, a sampling rate of just 2Hz is sufficient to resolve the IPS signature of most sources. We develop tests to ensure that IPS variability is separated from ionospheric or instrumental variability. We validate our results by comparison with existing catalogues of IPS sources, and near-contemporaneous observations by other IPS facilities. Using just five minutes of data, we produce catalogues at both 79MHz and 158MHz, each containing over 350 scintillating sources. At the field centre we detect approximately one scintillating source per square degree, with a minimum scintillating flux density at 158MHz of 110mJy, corresponding to a compact flux density of approximately 400mJy. Each of these sources is a known radio source, however only a minority were previously known to contain sub-arcsecond components. We discuss our findings and the prospects they hold for future astrophysical and heliospheric studies.

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 analyze the galactic cosmic ray (GCR) density and its spatial gradient in Forbush Decreases (FDs) observed with the Global Muon Detector Network (GMDN) and neutron monitors (NMs). By superposing the GCR density and density gradient observed in FDs following 45 interplanetary shocks (IP-shocks), each associated with an identified eruption on the sun, we infer the average spatial distribution of GCRs behind IP-shocks. We find two distinct modulations of GCR density in FDs, one in the magnetic sheath and the other in the coronal mass ejection (CME) behind the sheath. The density modulation in the sheath is dominant in the western flank of the shock, while the modulation in the CME ejecta stands out in the eastern flank. This east-west asymmetry is more prominent in GMDN data responding to $\sim$ 60 GV GCRs than in NM data responding to $\sim$ 10 GV GCRs, because of softer rigidity spectrum of the modulation in the CME ejecta than in the sheath. The GSE-y component of the density gradient, $G_y$ shows a negative (positive) enhancement in FDs caused by the eastern (western) eruptions, while $G_z$ shows a negative (positive) enhancement in FDs by the northern (southern) eruptions. This implies the GCR density minimum being located behind the central flank of IP-shock and propagating radially outward from location of the solar eruption. We also confirmed the average $G_z$ changing its sign above and below the heliospheric current sheet, in accord with the prediction of the drift model for the large-scale GCR transport in the heliosphere.

A detailed investigation was carried out to understand a prolonged (~44 hours) weakly southward interplanetary magnetic field (IMF-Bz) condition during May 02--04, 1998. In-situ observations, during the period, showed the passage of an expanding magnetic cloud embedded in an interplanetary coronal mass ejection (ICME), followed up by a shock and an interplanetary discontinuity driven by another ICME. It is the arrival of the ICMEs and the upfront shocks that cause the prolonged southward IMF-Bz condition. The magnetic configuration of the source regions of the IMF associated with the ICME interval were also examined, which showed open magnetic field structures, emanating from a small active region on the north of the heliospheric current sheet (HCS). The structures remained constantly to the north of the HCS, both on April 29 and May 01, suggesting no change in their polarity. The draping of these outward directed radial field lines around the propagating CMEs in the shocked plasma explained the observed polarity changes of the IMF-Bz at 1 AU. In addition, multiple enhancements were also detected in the geomagnetic field variations, which showed a distinct one-to-one correspondence with the density pulses observed at 1 AU, during 0700--1700 UT on May 03. The spectral analysis of both the variations showed the same discrete frequencies of 0.48, 0.65 and 0.75 mHz, demonstrating that the solar wind density enhancements can cause detectable global geomagnetic disturbances. The observations, thus, provide a deeper insight into the possible causes and geomagnetic consequences of a prolonged weakly southward IMF-Bz condition.

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.

We have studied the relationship between the solar-wind speed $[V]$ and the coronal magnetic-field properties (a flux expansion factor [$f$] and photospheric magnetic-field strength [$B_{\mathrm{S}}$]) at all latitudes using data of interplanetary scintillation and solar magnetic field obtained for 24 years from 1986 to 2009. Using a cross-correlation analyses, we verified that $V$ is inversely proportional to $f$ and found that $V$ tends to increase with $B_{\mathrm{S}}$ if $f$ is the same. As a consequence, we find that $V$ has extremely good linear correlation with $B_{\mathrm{S}}/f$. However, this linear relation of $V$ and $B_{\mathrm{S}}/f$ cannot be used for predicting the solar-wind velocity without information on the solar-wind mass flux. We discuss why the inverse relation between $V$ and $f$ has been successfully used for solar-wind velocity prediction, even though it does not explicitly include the mass flux and magnetic-field strength, which are important physical parameters for solar-wind acceleration.

We report observations of a steady 20 year decline of solar photospheric fields at latitudes $\geq$45${^{o}}$ starting from $\sim$1995. This prolonged and continuing decline, combined with the fact that Cycle 24 is already past its peak, implies that magnetic fields are likely to continue to decline until $\sim$2020, the expected minimum of the ongoing solar Cycle 24. In addition, interplanetary scintillation (IPS) observations of the inner heliosphere for the period 1983--2013 and in the distance range 0.2--0.8 AU, have also shown a similar and steady decline in solar wind micro-turbulence levels, in sync with the declining photospheric fields. Using the correlation between the polar field and heliospheric magnetic field (HMF) at solar minimum, we have estimated the value of the HMF in 2020 to be 3.9 ($\pm$0.6) and a floor value of the HMF of $\sim$3.2 ($\pm$0.4) nT. Given this floor value for the HMF, our analysis suggests that the estimated peak sunspot number for solar Cycle 25 is likely to be $\sim$ 62 ($\pm$12).

We report the electron density in a plasma tail of Comet ISON (C/2012 S1) derived from interplanetary scintillation (IPS) observations during November 1--28, 2013. Comet ISON showed a well-developed plasma tail (longer than 2.98 x 10^{7} km) before its perihelion passage on November 28. We identified a radio source whose line-of-sight approached the ISON's plasma tail in the above period and obtained its IPS data using the Solar Wind Imaging Facility at 327 MHz. We used the Heliospheric Imager onboard the Solar-Terrestrial Relation Observatory to distinguish between the cometary tail and solar eruption origins of their enhanced scintillation. From our examinations, we confirmed three IPS enhancements of a radio source 1148-00 on November 13, 16, and 17, which could be attributed to the disturbance in the cometary tail. Power spectra of 1148-00 had the steeper slope than normal ones during its occultation by the plasma tail. We estimated the electron density in the ISON's plasma tail and found 84 cm^{-3} around the tail axis at a distance of 3.74 x 10^{7} km from the cometary nucleus and an unexpected variation of the electron density in the vicinity of the tail boundary.

We report radial speed evolution of interplanetary coronal mass ejections (ICMEs) detected by the SOHO/LASCO coronagraph, interplanetary scintillation (IPS) at 327 MHz, and in-situ observations. In this study, we analyze solar wind disturbance factor (g-value) data derived from IPS observations during 1997-2009 covering nearly the whole period of Solar Cycle 23. By comparing observations from the SOHO/LASCO, IPS, and in situ, we then identify 39 ICMEs that could be analyzed carefully. Here, we define two speeds V_SOHO and V_bg, which are initial speed of ICME and the speed of the background solar wind, respectively. Examinations of these speeds yield the following results; i) Fast ICMEs (with V_SOHO - V_bg > 500 km/s) rapidly decelerate, moderate ICMEs (with 0 km/s < V_SOHO - V_bg < 500 km/s) show either gradually decelerating or uniform motion, and slow ICMEs (with V_SOHO - V_bg < 0 km/s) accelerate. The radial speeds converge on the speed of background solar wind during their outward propagation. We subsequently find; ii) both the acceleration and deceleration are nearly complete by 0.79 (+/- 0.04) AU, and those are ended when the ICME speed reaches a given speed. We find the value of that to be 480 (+/- 21) km/s. iii) For the fast and moderate ICMEs, a linear equation with a constant coefficient is more appropriate than a quadratic equation to describe their kinematics, because the chi-square for the linear equation satisfies the statistical significance level of 0.05, while the quadratic one is not. These results support the hypothesis that the radial motion of ICMEs is governed by a drag force due to interaction with the background solar wind. These findings also suggest that ICMEs propagating faster than background solar wind are controlled mainly by the hydrodynamic Stokes drag.

The ratio of the rms electron density fluctuations to the background density in the solar wind (density modulation index, $\epsilon_{N} \equiv \Delta{N}/N$) is of vital importance in understanding several problems in heliospheric physics related to solar wind turbulence. In this paper, we have investigated the behavior of $\epsilon_{N}$ in the inner-heliosphere from 0.26 to 0.82 AU. The density fluctuations $\Delta{N}$ have been deduced using extensive ground-based observations of interplanetary scintillation (IPS) at 327 MHz, which probe spatial scales of a few hundred km. The background densities ($N$) have been derived using near-Earth observations from the Advanced Composition Explorer ($\it{ACE}$). Our analysis reveals that $0.001 \lesssim \epsilon_{N} \lesssim 0.02$ and does not vary appreciably with heliocentric distance. We also find that $\epsilon_{N}$ declines by 8% from 1998 to 2008. We discuss the impact of these findings on problems ranging from our understanding of Forbush decreases to the behavior of the solar wind dynamic pressure over the recent peculiar solar minimum at the end of cycle 23..

We report kinematic properties of slow interplanetary coronal mass ejections (ICMEs) identified by SOHO/LASCO, interplanetary scintillation, and in-situ observations, and propose a modified equation for the ICME motion. We identify seven ICMEs between 2010 and 2011, and examine them with 39 events reported in our previous work. We examine 15 fast (V_SOHO - V_bg > 500 km/s), 25 moderate (0 km/s <= V_SOHO - V_bg <= 500 km/s), and 6 slow (V_SOHO - V_bg < 0 km/s) ICMEs, where V_SOHO and V_bg are the initial speed of ICMEs and the speed of the background solar wind, respectively. For slow ICMEs, we found the following results: i) They accelerate toward the speed of the background solar wind during their propagation, and reach their final speed by 0.34 (+-) 0.03 AU. ii) The acceleration ends when they reach 479 (+-) 126 km/s; this is close to the typical speed of the solar wind during the period of this study. iii) When k_1 and k_2 are assumed to be constants, a quadratic equation for the acceleration a = -k_2(V - V_bg)|V - V_bg| is more appropriate than a linear one a = -k_1(V - V_bg), where V is the propagation speed of ICMEs, while the latter gives a smaller chi-square value than the former. For the motion of the fast and moderate ICMEs, we found a modified drag equation. From the viewpoint of fluid dynamics, we interpret this equation as indicating that ICMEs with 0 km/s <= V - V_bg <= 2300 km/s are controlled mainly by the hydrodynamic Stokes drag force, while the aerodynamic drag force is a predominant factor for the propagation of ICME with V - V_bg > 2300 km/s.

The 3D structure of solar wind and its evolution in time is needed for heliospheric modeling and interpretation of energetic neutral atoms observations. We present a model to retrieve the solar wind structure in heliolatitude and time using all available and complementary data sources. We determine the heliolatitude structure of solar wind speed on a yearly time grid over the past 1.5 solar cycles based on remote-sensing observations of interplanetary scintillations, in situ out-of-ecliptic measurements from Ulysses, and in situ in-ecliptic measurements from the OMNI-2 database. Since the in situ information on the solar wind density structure out of ecliptic is not available apart from the Ulysses data, we derive correlation formulae between solar wind speed and density and use the information on the solar wind speed from interplanetary scintillation observations to retrieve the 3D structure of solar wind density. With the variations of solar wind density and speed in time and heliolatitude available we calculate variations in solar wind flux, dynamic pressure and charge exchange rate in the approximation of stationary H atoms.

The goal of the Fully Online Datacenter of Ultraviolet Emissions (FONDUE) Working Team of the International Space Science Institute in Bern, Switzerland, was to establish a common calibration of various UV and EUV heliospheric observations, both spectroscopic and photometric. Realization of this goal required an up-to-date model of spatial distribution of neutral interstellar hydrogen in the heliosphere, and to that end, a credible model of the radiation pressure and ionization processes was needed. This chapter describes the solar factors shaping the distribution of neutral interstellar H in the heliosphere. Presented are the solar Lyman-alpha flux and the solar Lyman-alpha resonant radiation pressure force acting on neutral H atoms in the heliosphere, solar EUV radiation and the photoionization of heliospheric hydrogen, and their evolution in time and the still hypothetical variation with heliolatitude. Further, solar wind and its evolution with solar activity is presented in the context of the charge exchange ionization of heliospheric hydrogen, and in the context of dynamic pressure variations. Also the electron ionization and its variation with time, heliolatitude, and solar distance is presented. After a review of all of those topics, we present an interim model of solar wind and the other solar factors based on up-to-date in situ and remote sensing observations of solar wind. Results of this effort will further be utilised to improve on the model of solar wind evolution, which will be an invaluable asset in all heliospheric measurements, including, among others, the observations of Energetic Neutral Atoms by the Interstellar Boundary Explorer (IBEX).

Extensive interplanetary scintillation (IPS) observations at 327 MHz obtained between 1983 and 2009 clearly show a steady and significant drop in the turbulence levels in the entire inner heliosphere starting from around ~1995. We believe that this large-scale IPS signature, in the inner heliosphere, coupled with the fact that solar polar fields have also been declining since ~1995, provide a consistent result showing that the buildup to the deepest minimum in 100 years actually began more than a decade earlier.

Observations of the polar region of the Sun are critically important for understanding the solar dynamo and the acceleration of solar wind. We carried out precise magnetic observations on both the North polar region and the quiet Sun at the East limb with the Spectro-Polarimeter of the Solar Optical Telescope aboard Hinode to characterize the polar region with respect to the quiet Sun. The average area and the total magnetic flux of the kG magnetic concentrations in the polar region appear to be larger than those of the quiet Sun. The magnetic field vectors classified as vertical in the quiet Sun have symmetric histograms around zero in the strengths, showing balanced positive and negative flux, while the histogram in the North polar region is clearly asymmetric, showing a predominance of the negative polarity. The total magnetic flux of the polar region is larger than that of the quiet Sun. In contrast, the histogram of the horizontal magnetic fields is exactly the same between the polar region and the quiet Sun. This is consistent with the idea that a local dynamo process is responsible for the horizontal magnetic fields. A high-resolution potential field extrapolation shows that the majority of magnetic field lines from the kG-patches in the polar region are open with a fanning-out structure very low in the atmosphere, while in the quiet Sun, almost all the field lines are closed.

We analyze the precursor of a Forbush Decrease (FD) observed with the Global Muon Detector Network on December 14, 2006. An intense geomagnetic storm is also recorded during this FD with the peak Kp index of 8+. By using the "two-dimensional map" of the cosmic ray intensity produced after removing the contribution from the diurnal anisotropy, we succeed in extracting clear signatures of the precursor. A striking feature of this event is that a weak loss-cone signature is first recorded more than a day prior to the Storm Sudden Commencement (SSC) onset. This suggests that the loss-cone precursor appeared only 7 hours after the Coronal Mass Ejection (CME) eruption from the Sun, when the Interplanetary (IP) shock driven by the Interplanetary Coronal Mass Ejection (ICME) located at 0.4 AU from the Sun. We find the precursor being successively observed with multiple detectors in the network according to the Earth's spin and confirmed that the precursor continuously exists in space. The long lead time (15.6 hours) of this precursor which is almost twice the typical value indicates that the IMF was more quiet in this event than a typical power spectrum assumed for the IMF turbulence. The amplitude (-6.45 %) of the loss-cone anisotropy at the SSC onset is more than twice the FD size, indicating that the maximum intensity depression behind the IP shock is much larger than the FD size recorded at the Earth in this event. We also find the excess intensity from the sunward IMF direction clearly observed during ~10 hours preceding the SSC onset. It is shown that this excess intensity is consistent with the measurement of the particles accelerated by the head-on collisions with the approaching shock. This is the first detailed observation of the precursor due to the shock reflected particles with muon detectors.

We deduce on hourly basis the spatial gradient of the cosmic ray density in three dimensions from the directional anisotropy of high-energy (~50 GeV) galactic cosmic ray (GCR) intensity observed with a global network of muon detectors on the Earth's surface. By analyzing the average features of the gradient in the corotational interaction regions (CIRs) recorded in successive two solar activity minimum periods, we find that the observed latitudinal gradient (Gz) changes its sign from negative to positive on the Earth's heliospheric current sheet (HCS) crossing from the northern to the southern hemisphere in A<0 epoch, while it changes from positive to negative in A>0 epoch. This is in accordance with the drift prediction. We also find a negative enhancement in Gx after the HCS crossing in both A<0 and A>0 epochs, but not in Gy. This asymmetrical feature of Gx and Gy indicates significant contributions from the parallel and perpendicular diffusions to the the gradient in CIRs in addition to the contribution from the drift effect.

Topics of the Colloquium: a) Interplanetary scintillation b) Interstellar scintillation c) Modeling and physical origin of the interplanetary and the interstellar plasma turbulence d) Scintillation as a tool for investigation of radio sources e) Seeing through interplanetary and interstellar turbulent media Ppt-presentations are available on the Web-site: http://www.prao.ru/conf/Colloquium/main.html