The study of ambient solar wind, a continuous pressure-driven plasma flow emanating from our Sun, is an important component of space weather research. The ambient solar wind flows in interplanetary space determine how solar storms evolve through the heliosphere before reaching Earth, and especially during solar minimum are themselves a driver of activity in the Earth's magnetic field. Accurately forecasting the ambient solar wind flow is therefore imperative to space weather awareness. Here we present a novel machine learning approach in which solutions from magnetic models of the solar corona are used to output the solar wind conditions near the Earth. The results are compared to observations and existing models in a comprehensive validation analysis, and the new model outperforms existing models in almost all measures. The final model discussed here represents an extremely fast, well-validated and open-source approach to the forecasting of ambient solar wind at Earth.

Visible coronal structure, in particular the spatial evolution of coronal streamers, provides indirect information about solar magnetic activity and the underlying solar dynamo. Their apparent absence of structure observed during the total eclipses of throughout the Maunder Minimum has been interpreted as evidence of a significant change in the solar magnetic field from that during modern cycles. Eclipse observations available from the more recent Dalton Minimum may be able to provide further information, sunspot activity being between the levels seen during recent cycles and in the Maunder Minimum. Here, we show and examine two graphical records of the total solar eclipse on 1806 June 16, during the Dalton Minimum. These records show significant rays and streamers around an inner ring. The ring is estimated to be ~ 0.44 R_S in width and the streamers in excess of 11.88 R_S in length. In combination with records of spicules or prominences, these eclipse records visually contrast the Dalton Minimum with the Maunder Minimum in terms of their coronal structure and support the existing discussions based on the sunspot observations. These eclipse records are broadly consistent with the solar cycle phase in the modelled open solar flux and the reconstructed slow solar wind at most latitudes.

Local inversions are often observed in the heliospheric magnetic field (HMF), but their origins and evolution are not yet fully understood.Parker Solar Probe has recently observed rapid, Alfvenic, HMF inversions in the inner heliosphere, known as 'switchbacks', which have been interpreted as the possible remnants of coronal jets. It has also been suggested that inverted HMF may be produced by near-Sun interchange reconnection; a key process in mechanisms proposed for slow solar wind release. These cases suggest that the source of inverted HMF is near the Sun, and it follows that these inversions would gradually decay and straighten as they propagate out through the heliosphere. Alternatively, HMF inversions could form during solar wind transit, through phenomena such velocity shears, draping over ejecta, or waves and turbulence. Such processes are expected to lead to a qualitatively radial evolution of inverted HMF structures. Using Helios measurements spanning 0.3-1 AU, we examine the occurrence rate of inverted HMF, as well as other magnetic field morphologies, as a function of radial distance A, and find that it continually increases. This trend may be explained by inverted HMF observed between 0.3-1 AU being primarily driven by one or more of the above in-transit processes, rather than created at the Sun. We make suggestions as to the relative importance of these different processes based on the evolution of the magnetic field properties associated with inverted HMF. We also explore alternative explanations outside of our suggested driving processes which may lead to the observed trend.

The rate at which the solar wind extracts angular momentum from the Sun has been predicted by theoretical models for many decades, and yet we lack a conclusive measurement from in-situ observations. In this letter we present a new estimate of the time-varying angular momentum flux in the equatorial solar wind, as observed by the \textit{Wind} spacecraft from 1994-2019. We separate the angular momentum flux into contributions from the protons, alpha particles, and magnetic stresses, showing that the mechanical flux in the protons is $\sim$3 times larger than the magnetic field stresses. We observe the tendency for the angular momentum flux of fast wind streams to be oppositely signed to the slow wind streams, as noted by previous authors. From the average total flux, we estimate the global angular momentum loss rate of the Sun to be $3.3\times10^{30}$erg, which lies within the range of various MHD wind models in the literature. This angular momentum loss rate is a factor of $\sim$2 weaker than required for a Skumanich-like rotation period evolution ($\Omega_*\propto $ stellar age$^{-1/2}$), which should be considered in studies of the rotation period evolution of Sun-like stars.

The Sun and Sun-like stars lose angular momentum to their magnetised stellar winds. This braking torque is coupled to the stellar magnetic field, such that changes in the strength and/or geometry of the field modifies the efficiency of this process. Since the space-age, we have been able to directly measure solar wind properties using in-situ spacecraft. Furthermore, indirect proxies such as sunspot number, geomagnetic indices, and cosmogenic radionuclides, constrain the variation of solar wind properties on centennial, and millennial timescales. We use near-Earth measurements of the solar wind plasma and magnetic field to calculate the torque on the Sun throughout the space-age. Then, reconstructions of the solar open magnetic flux are used to estimate the time-varying braking torque during the last nine millennia. We assume a relationship for the solar mass loss rate based on observations during the space-age which, due to the weak dependence of the torque on mass loss rate, does not strongly affect our predicted torque. The average torque during the last nine millennia is found to be $2.2\times10^{30}$erg, which is comparable to the average value from the last two decades. Our dataset includes grand minima (such as the Maunder Minimum), and maxima in solar activity, where the torque varies from $\sim1-5\times10^{30}$erg (averaged on decadal timescales), respectively. We find no evidence for any secular variation of the torque on timescales of less than $9000$ years.

Paper 1 [Lockwood et al., 2018] generated annual means of a new version of the $aa$ geomagnetic activity index which includes corrections for secular drift in the geographic coordinates of the auroral oval, thereby resolving the difference between the centennial-scale change in the northern and southern hemisphere indices, $aa_N$ and $aa_S$. However, other hemispheric asymmetries in the $aa$ index remain: in particular, the distributions of 3-hourly $aa_N$ and $aa_S$ values are different and the correlation between them is not high on this timescale ($r = 0.66$). In the present paper, a location-dependant station sensitivity model is developed using the $am$ index (derived from a much more extensive network of stations in both hemispheres) and used to reduce the difference between the hemispheric $aa$ indices and improve their correlation (to $r = 0.79$) by generating corrected 3-hourly hemispheric indices, $aa_{HN}$ and $aa_{HS}$, which also include the secular drift corrections detailed in Paper 1. These are combined into a new, 'homogeneous' $aa$ index, $aa_H$. It is shown that $aa_H$, unlike $aa$, reveals the 'equinoctial'-like time-of-day/time-of-year pattern that is found for the $am$ index.

Originally complied for 1868-1967 and subsequently continued so that it now covers 150 years, the $aa$ index has become a vital resource for studying space climate change. However, there have been debates about the inter-calibration of data from the different stations. In addition, the effects of secular change in the geomagnetic field have not previously been allowed for. As a result, the components of the 'classical' $aa$ index for the southern and northern hemispheres ($aa_S$ and $aa_N$) have drifted apart. We here separately correct both $aa_S$ and $aa_N$ for both these effects using the same method as used to generate the classic $aa$ values but allowing ${\delta}$, the minimum angular separation of each station from a nominal auroral oval, to vary as calculated using the IGRF-12 and gufm1 models of the intrinsic geomagnetic field. Our approach is to correct the quantized aK-values for each station, originally scaled on the assumption that ${\delta}$ values are constant, with time-dependent scale factors that allow for the drift in ${\delta}$. This requires revisiting the intercalibration of successive stations used in making the $aa_S$ and $aa_N$ composites. These intercalibrations are defined using independent data and daily averages from 11 years before and after each station change and it is shown that they depend on the time of year. This procedure produces new homogenized hemispheric aa indices, $aa_{HS}$ and $aa_{HN}$, which show centennial-scale changes that are in very close agreement. Calibration problems with the classic $aa$ index are shown to have arisen from drifts in ${\delta}$ combined with simpler corrections which gave an incorrect temporal variation and underestimate the rise in $aa$ during the 20th century by about 15%.

We produce light curves for all ~34,000 targets observed with K2 in Campaign 17 (C17), identifying 34 planet candidates, 184 eclipsing binaries, and 222 other periodic variables. The location of the C17 field means follow-up can begin immediately now that the campaign has concluded and interesting targets have been identified. The C17 field has a large overlap with C6, so this latest campaign also offers a rare opportunity to study a large number of targets already observed in a previous K2 campaign. The timing of the C17 data release, shortly before science operations begin with the Transiting Exoplanet Survey Satellite (TESS), also lets us exercise some of the tools and methods developed for identification and dissemination of planet candidates from TESS. We find excellent agreement between these results and those identified using only K2-based tools. Among our planet candidates are several planet candidates with sizes < 4 R_E and orbiting stars with KepMag < 10 (indicating good RV targets of the sort TESS hopes to find) and a Jupiter-sized single-transit event around a star already hosting a 6 d planet candidate.

Both direct observations and reconstructions from various datasets, suggest that conditions were radically different during the Maunder Minimum (MM) than during the space era. Using an MHD model, we develop a set of feasible solutions to infer the properties of the solar wind during this interval. Additionally, we use these results to drive a global magnetospheric model. Finally, using the 2008/2009 solar minimum as an upper limit for MM conditions, we use results from the International Reference Ionosphere (ILI) model to speculate on the state of the ionosphere. The results describe interplanetary, magnetospheric, and ionospheric conditions that were substantially different than today. For example: (1) the solar wind density and magnetic field strength were an order of magnitude lower; (2) the Earth's magnetopause and shock standoff distances were a factor of two larger; and (3) the maximum electron density in the ionosphere was substantially lower.

Using the reconstruction of power input to the magnetosphere given in Paper 1 (arXiv:1708.04904), we reconstruct annual means of geomagnetic indices over the past 400 years to within a 1-sigma error of +/-20 pc. In addition, we study the behaviour of the lognormal distribution of daily and hourly values about these annual means and show that we can also reconstruct the fraction of geomagnetically-active (storm-like) days and (substorm-like) hours in each year to accuracies of 50-60 pc. The results are the first physics-based quantification of the space weather conditions in both the Dalton and Maunder minima. We predict terrestrial disturbance levels in future repeats of these minima, allowing for the weakening of Earth's dipole moment.

The heliospheric magnetic field is of pivotal importance in solar and space physics. The field is rooted in the Sun's photosphere, where it has been observed for many years. Global maps of the solar magnetic field based on full disk magnetograms are commonly used as boundary conditions for coronal and solar wind models. Two primary observational constraints on the models are (1) the open field regions in the model should approximately correspond to coronal holes observed in emission, and (2) the magnitude of the open magnetic flux in the model should match that inferred from in situ spacecraft measurements. In this study, we calculate both MHD and PFSS solutions using fourteen different magnetic maps produced from five different types of observatory magnetograms, for the time period surrounding July, 2010. We have found that for all of the model/map combinations, models that have coronal hole areas close to observations underestimate the interplanetary magnetic flux, or, conversely, for models to match the interplanetary flux, the modeled open field regions are larger than coronal holes observed in EUV emission. In an alternative approach, we estimate the open magnetic flux entirely from solar observations by combining automatically detected coronal holes for Carrington rotation 2098 with observatory synoptic magnetic maps. This approach also underestimates the interplanetary magnetic flux. Our results imply that either typical observatory maps underestimate the Sun's magnetic flux, or a significant portion of the open magnetic flux is not rooted in regions that are obviously dark in EUV and X-ray emission.

Using information on geomagnetic activity, sunspot numbers and cosmogenic isotopes, supported by historic eclipse images and in conjunction with models, it has been possible to reconstruct annual means of solar wind speed and number density and heliospheric magnetic field (HMF) intensity since 1611, when telescopic observations of sunspots began. These models are developed and tuned using data recorded by near-Earth interplanetary spacecraft and by solar magnetograms over the past 53 years. In this paper, we use these reconstructions to quantify power input into the magnetosphere over the past 400 years. For each year, both the annual mean power input is computed and its distribution in daily means. This is possible because the distribution of daily values divided by the annual mean is shown to maintain the same lognormal form with a constant variance. This study is another important step towards the development of a physics-based, long-term climatology of space weather conditions.

We use 5 test data series to quantify putative discontinuities around 1946 in 5 annual-mean sunspot number or group number sequences. The series tested are: the original and new versions of the Wolf/Zurich/International sunspot number composite [$R_{ISNv1}$ and $R_{ISNv2}$] ; the corrected version of $R_{ISNv1}$ [$R_C$]; the backbone group number [$R_{BB}$]; and the group number composite [$R_{UEA}$]. The test data are: the group number $N_G$ and total sunspot area $A_G$ from the RGO photoheliographic data; the CaK index from re-analysis of MWO CaII K spectroheliograms; the group number from the MWO sunspot drawings, $N_{MWO}$; and ionospheric critical frequencies at Slough [$foF2$]. The test data all vary with sunspot numbers, in some cases non-linearly. Tests use both before-and-after fit-residual comparison and correlation methods, applied to intervals iterated to minimise errors and eliminate the effect of the discontinuity date. It is not assumed that the correction required is by a constant factor, nor even linear in sunspot number. A non-linear correction is required by $R_C$, $R_{BB}$ and $R_{ISNv1}$, but not by $R_{ISNv2}$ or $R_{UEA}$. The test datasets give very similar results in all cases. By multiplying the probability distribution functions together we obtain the optimum correction for each data series that must be applied to pre-discontinuity data to make them consistent with the post-discontinuity data. It is shown that, on average, values for 1932-1943 are too small (relative to later values) by about 12.3% for $R_{ISNv1}$ but are too large for $R_{ISNv2}$ and $R_{BB}$ by 3.8% and 5.2%. The correction applied to generate $R_C$ from $R_{ISNv1}$ reduces this average factor to 0.5% but does not remove the non-linear variation, and other errors remain uncorrected. A test is provided by $R_{UEA}$, which is identical to the RGO $N_G$ values over the interval used.

We compare four sunspot-number data sequences against geomagnetic and terrestrial auroral observations. The comparisons are made for the original SIDC composite of Wolf-Zurich-International sunspot number [$R_{ISNv1}$], the group sunspot number [$R_{G}$] by Hoyt and Schatten (Solar Phys., 1998), the new "backbone" group sunspot number [$R_{BB}$] by Svalgaard and Schatten (Solar Phys., 2016), and the "corrected" sunspot number [$R_{C}$] by Lockwood at al. (J.G.R., 2014). Each sunspot number is fitted with terrestrial observations, or parameters derived from terrestrial observations to be linearly proportional to sunspot number, over a 30-year calibration interval of 1982-2012. The fits are then used to compute test sequences, which extend further back in time and which are compared to $R_{ISNv1}$, $R_{G}$, $R_{BB}$, and $R_{C}$. To study the long-term trends, comparisons are made using averages over whole solar cycles (minimum-to-minimum). The test variations are generated in four ways: i) using the IDV(1d) and IDV geomagnetic indices (for 1845-2013) fitted over the calibration interval using the various sunspot numbers and the phase of the solar cycle; ii) from the open solar flux (OSF) generated for 1845 - 2013 from four pairings of geomagnetic indices by Lockwood et al. (Ann. Geophys., 2014) and analysed using the OSF continuity model of Solanki at al. (Nature, 2000) which employs a constant fractional OSF loss rate; iii) the same OSF data analysed using the OSF continuity model of Owens and Lockwood (J.G.R., 2012) in which the fractional loss rate varies with the tilt of the heliospheric current sheet and hence with the phase of the solar cycle; iv) the occurrence frequency of low-latitude aurora for 1780-1980 from the survey of Legrand and Simon (Ann. Geophys., 1987). For all cases, $R_{BB}$ exceeds the test terrestrial series by an amount that increases as one goes back in time.

New sunspot data composites, some of which are radically different in the character of their long-term variation, are evaluated over the interval 1845-2014. The method commonly used to calibrate historic sunspot data, relative to modern-day data, is "daisy-chaining", whereby calibration is passed from one data subset to the neighbouring one, usually using regressions of the data subsets for the intervals of their overlap. Recent studies have illustrated serious pitfalls in these regressions and the resulting errors can be compounded by their repeated use as the data sequence is extended back in time. Hence the recent composite data series by Usoskin et al. (2016), $R_{UEA}$, is a very important advance because it avoids regressions, daisy-chaining and other common, but invalid, assumptions: this is achieved by comparing the statistics of "active day" fractions to those for a single reference dataset. We study six sunspot data series including $R_{UEA}$ and the new "backbone" data series $R_{BB}$, recently generated by Svalgaard and Schatten (2016) by employing both regression and daisy-chaining. We show that all six can be used with a continuity model to reproduce the main features of the open solar flux variation for 1845-2014, as reconstructed from geomagnetic activity data. However, some differences can be identified that are consistent with tests using a basket of other proxies for solar magnetic fields. Using data from a variety of sunspot observers, we illustrate problems with the method employed in $R_{BB}$ which cause it to increasingly overestimate sunspot numbers going back in time and we recommend using $R_{UEA}$ because it employs more robust procedures that avoid such problems.

More than 70 years ago it was recognised that ionospheric F2-layer critical frequencies $foF2$ had a strong relationship to sunspot number. Using historic datasets from the Slough and Washington ionosondes, we evaluate the best statistical fits of $foF2$ to sunspot numbers (at each Universal Time [UT] separately) in order to search for drifts and abrupt changes in the fit residuals over Solar Cycles 17 - 21. Polynomial fits are made both with and without allowance for the white-light facular area, which has been reported as being associated with cycle-to-cycle changes in the sunspot number - $foF2$ relationship. Over the interval studied here, the ISN, $R$, the backbone group number $Rbb$, and the corrected number $Rc$ largely differ in their allowance for the 'Waldmeier discontinuity' around 1945 (the correction factor for which for $R$, $Rbb$ and $Rc$ is, respectively, zero, effectively over 20%, and explicitly 11.6%). It is shown that for Solar Cycles 18 - 21, all three sunspot data sequences perform well, but that the fit residuals are lowest and most uniform for $Rbb$. We here use $foF2$ for those UTs for which $R$, $Rbb$, and $Rc$ all give correlations exceeding 0.99 for intervals both before and after the Waldmeier discontinuity. The error introduced by the Waldmeier discontinuity causes $R$ to underestimate the fitted values based on the $foF2$ data for 1932 - 1945 but $Rbb$ overestimates them by almost the same factor, implying that the correction for the Waldmeier discontinuity inherent in $Rbb$ is too large by a factor of two. Fit residuals are smallest and most uniform for $Rc$ and the ionospheric data support the optimum discontinuity multiplicative correction factor derived from the independent Royal Greenwich Observatory (RGO) sunspot group data for the same interval.

We use sunspot group observations from the Royal Greenwich Observatory (RGO) to investigate the effects of intercalibrating data from observers with different visual acuities. The tests are made by counting the number of groups $R_B$ above a variable cut-off threshold of observed total whole-spot area (uncorrected for foreshortening) to simulate what a lower acuity observer would have seen. The synthesised annual means of $R_B$ are then re-scaled to the observed RGO group number $R_A$ using a variety of regression techniques. It is found that a very high correlation between $R_A$ and $R_B$ ($r_{AB}$ > 0.98) does not prevent large errors in the intercalibration (e.g. sunspot maximum values can be over 30% too large even for such levels of $r_{AB}$). In generating the backbone sunspot number, Svalgaard and Schatten [2015] force regression fits to pass through the scatter plot origin which generates unreliable fits (the residuals do not form a normal distribution) and causes sunspot cycle amplitudes to be exaggerated in the intercalibrated data. It is demonstrated that the use of Quantile-Quantile (Q-Q) plots to test for a normal distribution is a useful indicator of erroneous and misleading regression fits. Ordinary least squares linear fits, not forced to pass through the origin, are sometimes reliable (although the optimum method used is shown to be different when matching peak and average sunspot group numbers). However other fits are only reliable if non-linear regression is used. From these results it is entirely possible that the inflation of solar cycle amplitudes in the backbone group sunspot number as one goes back in time, relative to related solar-terrestrial parameters, is entirely caused by the use of inappropriate and non-robust regression techniques to calibrate the sunspot data.

Although the sunspot-number series have existed since the mid-19th century, they are still the subject of intense debate, with the largest uncertainty being related to the "calibration" of the visual acuity of individual observers in the past. Daisy-chain regression methods are applied to inter-calibrate the observers which may lead to significant bias and error accumulation. Here we present a novel method to calibrate the visual acuity of the key observers to the reference data set of Royal Greenwich Observatory sunspot groups for the period 1900-1976, using the statistics of the active-day fraction. For each observer we independently evaluate their observational thresholds [S_S] defined such that the observer is assumed to miss all of the groups with an area smaller than S_S and report all the groups larger than S_S. Next, using a Monte-Carlo method we construct, from the reference data set, a correction matrix for each observer. The correction matrices are significantly non-linear and cannot be approximated by a linear regression or proportionality. We emphasize that corrections based on a linear proportionality between annually averaged data lead to serious biases and distortions of the data. The correction matrices are applied to the original sunspot group records for each day, and finally the composite corrected series is produced for the period since 1748. The corrected series displays secular minima around 1800 (Dalton minimum) and 1900 (Gleissberg minimum), as well as the Modern grand maximum of activity in the second half of the 20th century. The uniqueness of the grand maximum is confirmed for the last 250 years. It is shown that the adoption of a linear relationship between the data of Wolf and Wolfer results in grossly inflated group numbers in the 18th and 19th centuries in some reconstructions.

Aims: Although the time of the Maunder minimum (1645--1715) is widely known as a period of extremely low solar activity, claims are still debated that solar activity during that period might still have been moderate, even higher than the current solar cycle #24. We have revisited all the existing pieces of evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum. Methods: We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum. Results: The level of solar activity during the Maunder minimum is reassessed on the basis of all available data sets. Conclusions: We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely below that during the Dalton minimum around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected at a high confidence level.

Coronal mass ejections (CMEs) can be continuously tracked through a large portion of the inner heliosphere by direct imaging in visible and radio wavebands. White-light (WL) signatures of solar wind transients, such as CMEs, result from Thomson scattering of sunlight by free electrons, and therefore depend on both the viewing geometry and the electron density. The Faraday rotation (FR) of radio waves from extragalactic pulsars and quasars, which arises due to the presence of such solar wind features, depends on the line-of-sight magnetic field component $B_\parallel$, and the electron density. To understand coordinated WL and FR observations of CMEs, we perform forward magnetohydrodynamic modelling of an Earth-directed shock and synthesise the signatures that would be remotely sensed at a number of widely distributed vantage points in the inner heliosphere. Removal of the background solar wind contribution reveals the shock-associated enhancements in WL and FR. While the efficiency of Thomson scattering depends on scattering angle, WL radiance $I$ decreases with heliocentric distance $r$ roughly according to the expression $I \propto r^{-3}$. The sheath region downstream of the Earth-directed shock is well viewed from the L4 and L5 Lagrangian points, demonstrating the benefits of these points in terms of space weather forecasting. The spatial position of the main scattering site $\mathbf{r}_{\rm sheath}$ and the mass of plasma at that position $M_{\rm sheath}$ can be inferred from the polarisation of the shock-associated enhancement in WL radiance. From the FR measurements, the local $B_{\parallel {\rm sheath}}$ at $\mathbf{r}_{\rm sheath}$ can then be estimated. Simultaneous observations in polarised WL and FR can not only be used to detect CMEs, but also to diagnose their plasma and magnetic field properties.

The ability to predict times of greater galactic cosmic ray (GCR) fluxes is important for reducing the hazards caused by these particles to satellite communications, aviation, or astronauts. The 11-year solar cycle variation in cosmic rays is highly correlated with the strength of the heliospheric magnetic field. Differences in GCR flux during alternate solar cycles yield a 22-year cycle, known as the Hale Cycle, which is thought to be due to different particle drift patterns when the northern solar pole has predominantly positive (denoted a qA>0 cycle) or negative (qA<0) polarities. This results in the onset of the peak cosmic ray flux at Earth occurring earlier during qA>0 cycles than for qA<0 cycles and hence the peak being more domed for qA>0 and more sharply peaked for qA<0. In this study, we demonstrate that properties of the large-scale heliospheric magnetic field are different during the declining phase of the qA<0 and qA>0 solar cycles, when the difference in GCR flux is most apparent. This suggests that particle drifts may not be the sole mechanism responsible for the Hale Cycle in GCR flux at Earth. However, it is also demonstrated that these polarity-dependent heliospheric differences are evident during the space-age but much less clear in earlier data: using geomagnetic reconstructions, it is shown that for the period of 1905 - 1965, alternate polarities do not give as significant a difference during the declining phase of the solar cycle. Thus we suggest that the 22-year cycle in cosmic ray flux is at least partly the result of direct modulation by the heliospheric magnetic field and that this effect may be primarily limited to the grand solar maximum of the space-age.

We investigate a CME propagating towards Earth on 29 March 2011. This event is specifically chosen for its predominately northward directed magnetic field, so that the influence from the momentum flux onto Earth can be isolated. We focus our study on understanding how a small Earth-directed segment propagates. Mass images are created from the white-light cameras onboard STEREO which are also converted into mass height-time maps (mass J-maps). The mass tracks on these J-maps correspond to the sheath region between the CME and its associated shock front as detected by in situ measurements at L1. A time-series of mass measurements from the STEREO COR-2A instrument are made along the Earth propagation direction. Qualitatively, this mass time-series shows a remarkable resemblance to the L1 in situ density series. The in situ measurements are used as inputs into a 3D magnetospheric space weather simulation from CCMC. These simulations display a sudden compression of the magnetosphere from the large momentum flux at the leading edge of the CME and predictions are made for the time-derivative of the magnetic field (dB/dt) on the ground. The predicted dB/dt were then compared with observations from specific equatorially-located ground stations and show notable similarity. This study of the momentum of a CME from the Sun down to its influence on magnetic ground stations on Earth is presented as preliminary proof of concept, such that future attempts may try to use remote sensing to create density and velocity time-series as inputs to magnetospheric simulations.

Stereoscopic white-light imaging of a large portion of the inner heliosphere has been used to track interplanetary coronal mass ejections. At large elongations from the Sun, the white-light brightness depends on both the local electron density and the efficiency of the Thomson-scattering process. To quantify the effects of the Thomson-scattering geometry, we study an interplanetary shock using forward magnetohydrodynamic simulation and synthetic white-light imaging. Identifiable as an inclined streak of enhanced brightness in a time-elongation map, the travelling shock can be readily imaged by an observer located within a wide range of longitudes in the ecliptic. Different parts of the shock front contribute to the imaged brightness pattern viewed by observers at different longitudes. Moreover, even for an observer located at a fixed longitude, a different part of the shock front will contribute to the imaged brightness at any given time. The observed brightness within each imaging pixel results from a weighted integral along its corresponding ray-path. It is possible to infer the longitudinal location of the shock from the brightness pattern in an optical sky map, based on the east-west asymmetry in its brightness and degree of polarization. Therefore, measurement of the interplanetary polarized brightness could significantly reduce the ambiguity in performing three-dimensional reconstruction of local electron density from white-light imaging.

Recent coordinated observations of interplanetary scintillation (IPS) and stereoscopic heliospheric imagers (HIs) are significant to continuously track the propagation and evolution of solar eruptions throughout interplanetary space. In order to obtain a better understanding of the observational signatures in these two remote-sensing techniques, the magnetohydrodynamics of the macro-scale interplanetary disturbance and the radio-wave scattering of the micro-scale electron-density fluctuation are coupled and investigated using a newly-constructed multi-scale numerical model. This model is then applied to a case of an interplanetary shock propagation within the ecliptic plane. The shock could be nearly invisible to an HI, once entering the Thomson-scattering sphere of the HI. The asymmetry in the optical images between the western and eastern HIs suggests the shock propagation off the Sun-Earth line. Meanwhile, an IPS signal, strongly dependent on the local electron density, is insensitive to the density cavity far downstream of the shock front. When this cavity (or the shock nose) is cut through by an IPS ray-path, a single speed component at the flank (or the nose) of the shock can be recorded; when an IPS ray-path penetrates the sheath between the shock nose and this cavity, two speed components at the sheath and flank can be detected. Moreover, once a shock front touches an IPS ray-path, the derived position and speed at the irregularity source of this IPS signal, together with an assumption of a radial and constant propagation of the shock, can be used to estimate the later appearance of the shock front in the elongation of the HI field of view. The results of synthetic measurements from forward modelling are helpful in inferring the in-situ properties of coronal mass ejection from real observational data via an inverse approach.

We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the May 12, 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux rope onto the steady-state solution. During the initial stage of CME expansion, the core flux rope reconnects with the neighboring field, which facilitates lateral expansion of the CME footprint in the low corona. The flux rope field also reconnects with the oppositely orientated overlying magnetic field in the manner of the breakout model. During this stage of the eruption, the simulated CME rotates counter-clockwise to achieve an orientation that is in agreement with the interplanetary flux rope observed at 1 AU. A significant component of the CME that expands into interplanetary space comprises one of the side lobes created mainly as a result of reconnection with the overlying field. Within 3 hours, reconnection effectively modifies the CME connectivity from the initial condition where both footpoints are rooted in the active region to a situation where one footpoint is displaced into the quiet Sun, at a significant distance ($\approx 1R_\odot$) from the original source region. The expansion and rotation due to interaction with the overlying magnetic field stops when the CME reaches the outer edge of the helmet streamer belt, where the field is organized on a global scale. The simulation thus offers a new view of the role reconnection plays in rotating a CME flux rope and transporting its footpoints while preserving its core structure.

An update of Owens et al. (2008) shows that the relationship between the coronal mass ejection (CME) rate and the heliospheric magnetic field strength predicts a field floor of less than 4 nT at 1 AU. This implies that the record low values measured during this solar minimum do not necessarily contradict the idea that open flux is conserved. The results are consistent with the hypothesis that CMEs add flux to the heliosphere and interchange reconnection between open flux and closed CME loops subtracts flux. An existing model embracing this hypothesis, however, overestimates flux during the current minimum, even though the CME rate has been low. The discrepancy calls for reasonable changes in model assumptions.