Solar physics: shining a light on solar impacts

Comparing changes in temperature and solar radiation on centennial timescales can help to constrain the Sun’s impact on climate. New findings regarding the minimum activity level of the Sun reveal that comparisons made so far may have been too simplistic.

The variation of solar wind correlation lengths over three solar cycles

We present the results of a study of solar wind velocity and magnetic field correlation lengths over the last 35 years. The correlation length of the magnetic field magnitude λ | B| increases on average by a factor of two at solar maxima compared to solar minima. The correlation lengths of the components of the magnetic field λ_{B_{XYZ}} and of the velocity λ_{V_{YZ}} do not show this change and have similar values, indicating a continual turbulent correlation length of around 1.4×106 km. We conclude that a linear relation between λ | B|, VB 2, and Kp suggests that the former is related to the total magnetic energy in the solar wind and an estimate of the average size of geoeffective structures, which is, in turn, proportional to VB 2. By looking at the distribution of daily correlation lengths we show that the solar minimum values of λ | B| correspond to the turbulent outer scale. A tail of larger λ | B| values is present at solar maximum causing the increase in mean value.

A multispacecraft analysis of a small-scale transient entrained by solar wind streams

The images taken by the Heliospheric Imagers (HIs), part of the SECCHI imaging package onboard the pair of STEREO spacecraft, provide information on the radial and latitudinal evolution of the plasma compressed inside corotating interaction regions (CIRs). A plasma density wave imaged by the HI instrument onboard STEREO-B was found to propagate towards STEREO-A, enabling a comparison between simultaneous remotesensing and in situ observations of its structure to be performed. In situ measurements made by STEREO-A show that the plasma density wave is associated with the passage of a CIR. The magnetic field compressed after the CIR stream interface (SI) is found to have a planar distribution. Minimum variance analysis of the magnetic field vectors shows that the SI is inclined at 54° to the orbital plane of the STEREO-A spacecraft. This inclination of the CIR SI is comparable to the inclination of the associated plasma density wave observed by HI. A small-scale magnetic cloud with a flux rope topology and radial extent of 0.08 AU is also embedded prior to the SI. The pitch-angle distribution of suprathermal electrons measured by the STEREO-A SWEA instrument shows that an open magnetic field topology in the cloud replaced the heliospheric current sheet locally. These observations confirm that HI observes CIRs in difference images when a small-scale transient is caught up in the compression region.

Magnetic discontinuities in the near-earth solar wind: evidence of in-transit turbulence or remnants of coronal structure?

Fluctuations in the solar wind plasma and magnetic field are well described by the sum of two power law distributions. It has been postulated that these distributions are the result of two independent processes: turbulence, which contributes mainly to the smaller fluctuations, and crossing the boundaries of flux tubes of coronal origin, which dominates the larger variations. In this study we explore the correspondence between changes in the magnetic field with changes in other solar wind properties. Changes in density and temperature may result from either turbulence or coronal structures, whereas changes in composition, such as the alpha-to-proton ratio are unlikely to arise from in-transit effects. Observations spanning the entire ACE dataset are compared with a null hypothesis of no correlation between magnetic field discontinuities and changes in other solar wind parameters. Evidence for coronal structuring is weaker than for in-transit turbulence, with only ∼ 25% of large magnetic field discontinuities associated with a significant change in the alpha-to-proton ratio, compared to ∼ 40% for significant density and temperature changes. However, note that a lack of detectable alpha-to-proton signature is not sufficient to discount a structure as having a solar origin.

Reconstruction and prediction of variations in the open solar magnetic flux and interplanetary conditions

Historic geomagnetic activity observations have been used to reveal centennial variations in the open solar flux and the near-Earth heliospheric conditions (the interplanetary magnetic field and the solar wind speed). The various methods are in very good agreement for the past 135 years when there were sufficient reliable magnetic observatories in operation to eliminate problems due to site-specific errors and calibration drifts. This review underlines the physical principles that allow these reconstructions to be made, as well as the details of the various algorithms employed and the results obtained. Discussion is included of: the importance of the averaging timescale; the key differences between “range” and “interdiurnal variability” geomagnetic data; the need to distinguish source field sector structure from heliospherically-imposed field structure; the importance of ensuring that regressions used are statistically robust; and uncertainty analysis. The reconstructions are exceedingly useful as they provide calibration between the in-situ spacecraft measurements from the past five decades and the millennial records of heliospheric behaviour deduced from measured abundances of cosmogenic radionuclides found in terrestrial reservoirs. Continuity of open solar flux, using sunspot number to quantify the emergence rate, is the basis of a number of models that have been very successful in reproducing the variation derived from geomagnetic activity. These models allow us to extend the reconstructions back to before the development of the magnetometer and to cover the Maunder minimum. Allied to the radionuclide data, the models are revealing much about how the Sun and heliosphere behaved outside of grand solar maxima and are providing a means of predicting how solar activity is likely to evolve now that the recent grand maximum (that had prevailed throughout the space age) has come to an end.

The evolution of the sun’s open magnetic flux: II. full solar cycle simulations

In this paper the origin and evolution of the Sun’s open magnetic flux is considered by conducting magnetic flux transport simulations over many solar cycles. The simulations include the effects of differential rotation, meridional flow and supergranular diffusion on the radial magnetic field at the surface of the Sun as new magnetic bipoles emerge and are transported poleward. In each cycle the emergence of roughly 2100 bipoles is considered. The net open flux produced by the surface distribution is calculated by constructing potential coronal fields with a source surface from the surface distribution at regular intervals. In the simulations the net open magnetic flux closely follows the total dipole component at the source surface and evolves independently from the surface flux. The behaviour of the open flux is highly dependent on meridional flow and many observed features are reproduced by the model. However, when meridional flow is present at observed values the maximum value of the open flux occurs at cycle minimum when the polar caps it helps produce are the strongest. This is inconsistent with observations by Lockwood, Stamper and Wild (1999) and Wang, Sheeley, and Lean (2000) who find the open flux peaking 1–2 years after cycle maximum. Only in unrealistic simulations where meridional flow is much smaller than diffusion does a maximum in open flux consistent with observations occur. It is therefore deduced that there is no realistic parameter range of the flux transport variables that can produce the correct magnitude variation in open flux under the present approximations. As a result the present standard model does not contain the correct physics to describe the evolution of the Sun’s open magnetic flux over an entire solar cycle. Future possible improvements in modeling are suggested.

The formation of large-scale current sheets within magnetic clouds

Magnetic clouds are a class of interplanetary coronal mass ejections (CME) predominantly characterised by a smooth rotation in the magnetic field direction, indicative of a magnetic flux rope structure. Many magnetic clouds, however, also contain sharp discontinuities within the smoothly varying magnetic field, suggestive of narrow current sheets. In this study we present observations and modelling of magnetic clouds with strong current sheet signatures close to the centre of the apparent flux rope structure. Using an analytical magnetic flux rope model, we demonstrate how such current sheets can form as a result of a cloud’s kinematic propagation from the Sun to the Earth, without any external forces or influences. This model is shown to match observations of four particular magnetic clouds remarkably well. The model predicts that current sheet intensity increases for increasing CME angular extent and decreasing CME radial expansion speed. Assuming such current sheets facilitate magnetic reconnection, the process of current sheet formation could ultimately lead a single flux rope becoming fragmented into multiple flux ropes. This change in topology has consequences for magnetic clouds as barriers to energetic particle propagation.

From the sun to the Earth: The 13 May 2005 coronal mass ejection

We report the results of a multi-instrument, multi-technique, coordinated study of the solar eruptive event of 13 May 2005. We discuss the resultant Earth-directed (halo) coronal mass ejection (CME), and the effects on the terrestrial space environment and upper Earth atmosphere. The interplanetary CME (ICME) impacted the Earth’s magnetosphere and caused the most-intense geomagnetic storm of 2005 with a Disturbed Storm Time (Dst) index reaching −263 nT at its peak. The terrestrial environment responded to the storm on a global scale. We have combined observations and measurements from coronal and interplanetary remote-sensing instruments, interplanetary and near-Earth in-situ measurements, remote-sensing observations and in-situ measurements of the terrestrial magnetosphere and ionosphere, along with coronal and heliospheric modelling. These analyses are used to trace the origin, development, propagation, terrestrial impact, and subsequent consequences of this event to obtain the most comprehensive view of a geo-effective solar eruption to date. This particular event is also part of a NASA-sponsored Living With a Star (LWS) study and an on-going US NSF-sponsored Solar, Heliospheric, and INterplanetary Environment (SHINE) community investigation.

Assessing the accuracy of CME speed and trajectory estimates from STEREO observations through a comparison of independent methods

We have estimated the speed and direction of propagation of a number of Coronal Mass Ejections (CMEs) using single-spacecraft data from the STEREO Heliospheric Imager (HI) wide-field cameras. In general, these values are in good agreement with those predicted by Thernisien, Vourlidas, and Howard in Solar Phys. 256, 111 -aEuro parts per thousand 130 (2009) using a forward modelling method to fit CMEs imaged by the STEREO COR2 coronagraphs. The directions of the CMEs predicted by both techniques are in good agreement despite the fact that many of the CMEs under study travel in directions that cause them to fade rapidly in the HI images. The velocities estimated from both techniques are in general agreement although there are some interesting differences that may provide evidence for the influence of the ambient solar wind on the speed of CMEs. The majority of CMEs with a velocity estimated to be below 400 km s(-1) in the COR2 field of view have higher estimated velocities in the HI field of view, while, conversely, those with COR2 velocities estimated to be above 400 km s(-1) have lower estimated HI velocities. We interpret this as evidence for the deceleration of fast CMEs and the acceleration of slower CMEs by interaction with the ambient solar wind beyond the COR2 field of view. We also show that the uncertainties in our derived parameters are influenced by the range of elongations over which each CME can be tracked. In order to reduce the uncertainty in the predicted arrival time of a CME at 1 Astronomical Unit (AU) to within six hours, the CME needs to be tracked out to at least 30 degrees elongation. This is in good agreement with predictions of the accuracy of our technique based on Monte Carlo simulations.

The Heliospheric Imagers onboard the STEREO mission

Mounted on the sides of two widely separated spacecraft, the two Heliospheric Imager (HI) instruments onboard NASA’s STEREO mission view, for the first time, the space between the Sun and Earth. These instruments are wide-angle visible-light imagers that incorporate sufficient baffling to eliminate scattered light to the extent that the passage of solar coronal mass ejections (CMEs) through the heliosphere can be detected. Each HI instrument comprises two cameras, HI-1 and HI-2, which have 20° and 70° fields of view and are off-pointed from the Sun direction by 14.0° and 53.7°, respectively, with their optical axes aligned in the ecliptic plane. This arrangement provides coverage over solar elongation angles from 4.0° to 88.7° at the viewpoints of the two spacecraft, thereby allowing the observation of Earth-directed CMEs along the Sun – Earth line to the vicinity of the Earth and beyond. Given the two separated platforms, this also presents the first opportunity to view the structure and evolution of CMEs in three dimensions. The STEREO spacecraft were launched from Cape Canaveral Air Force Base in late October 2006, and the HI instruments have been performing scientific observations since early 2007. The design, development, manufacture, and calibration of these unique instruments are reviewed in this paper. Mission operations, including the initial commissioning phase and the science operations phase, are described. Data processing and analysis procedures are briefly discussed, and ground-test results and in-orbit observations are used to demonstrate that the performance of the instruments meets the original scientific requirements.

Pre-CME onset fuses: do the STEREO heliospheric imagers hold the clues to the CME onset process?

Understanding the onset of coronal mass ejections (CMEs) is surely one of the holy grails of solar physics today. Inspection of data from the Heliospheric Imagers (HI), which are part of the SECCHI instrument suite aboard the two NASA STEREO spacecraft, appears to have revealed pre-eruption signatures which may provide valuable evidence for identifying the CME onset mechanism. Specifically, an examination of the HI images has revealed narrow rays comprised of a series of outward-propagating plasma blobs apparently forming near the edge of the streamer belt prior to many CME eruptions. In this pilot study, we inspect a limited dataset to explore the significance of this phenomenon, which we have termed a pre-CME ‘fuse’. Although, the enhanced expulsion of blobs may be consistent with an increase in the release of outward-propagating blobs from the streamers themselves, it could also be interpreted as evidence for interchange reconnection in the period leading to a CME onset. Indeed, it is argued that the latter could even have implications for the end-of-life of CMEs. Thus, the presence of these pre-CME fuses provides evidence that the CME onset mechanism is either related to streamer reconnection processes or the reconnection between closed field lines in the streamer belt and adjacent, open field lines. We investigate the nature of these fuses, including their timing and location with respect to CME launch sites, as well as their speed and topology.

In situ signatures of interchange reconnection between magnetic clouds and open magnetic fields: a mechanism for the erosion of polar coronal holes?

We outline a method to determine the direction of solar open flux transport that results from the opening of magnetic clouds (MCs) by interchange reconnection at the Sun based solely on in-situ observations. This method uses established findings about i) the locations and magnetic polarities of emerging MC footpoints, ii) the hemispheric dependence of the helicity of MCs, and iii) the occurrence of interchange reconnection at the Sun being signaled by uni-directional suprathermal electrons inside MCs. Combining those observational facts in a statistical analysis of MCs during solar cycle 23 (period 1995 – 2007), we show that the time of disappearance of the northern polar coronal hole (1998 – 1999), permeated by an outward-pointing magnetic field, is associated with a peak in the number of MCs originating from the northern hemisphere and connected to the Sun by outward-pointing magnetic field lines. A similar peak is observed in the number of MCs originating from the southern hemisphere and connected to the Sun by inward-pointing magnetic field lines. This pattern is interpreted as the result of interchange reconnection occurring between MCs and the open field lines of nearby polar coronal holes. This reconnection process closes down polar coronal hole open field lines and transports these open field lines equatorward, thus contributing to the global coronal magnetic field reversal process. These results will be further constrainable with the rising phase of solar cycle 24.

Implications of non-cylindrical flux ropes for magnetic cloud reconstruction techniques and the interpretation of double flux rope events

Magnetic clouds (MCs) are a subset of interplanetary coronal mass ejections (ICMEs) which exhibit signatures consistent with a magnetic flux rope structure. Techniques for reconstructing flux rope orientation from single-point in situ observations typically assume the flux rope is locally cylindrical, e.g., minimum variance analysis (MVA) and force-free flux rope (FFFR) fitting. In this study, we outline a non-cylindrical magnetic flux rope model, in which the flux rope radius and axial curvature can both vary along the length of the axis. This model is not necessarily intended to represent the global structure of MCs, but it can be used to quantify the error in MC reconstruction resulting from the cylindrical approximation. When the local flux rope axis is approximately perpendicular to the heliocentric radial direction, which is also the effective spacecraft trajectory through a magnetic cloud, the error in using cylindrical reconstruction methods is relatively small (≈ 10∘). However, as the local axis orientation becomes increasingly aligned with the radial direction, the spacecraft trajectory may pass close to the axis at two separate locations. This results in a magnetic field time series which deviates significantly from encounters with a force-free flux rope, and consequently the error in the axis orientation derived from cylindrical reconstructions can be as much as 90∘. Such two-axis encounters can result in an apparent ‘double flux rope’ signature in the magnetic field time series, sometimes observed in spacecraft data. Analysing each axis encounter independently produces reasonably accurate axis orientations with MVA, but larger errors with FFFR fitting.

Effects of Thomson-scattering geometry on white-light imaging of an interplanetary shock: synthetic observations from forward magnetohydrodynamic modelling

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 polarisation. Therefore, measurement of the interplanetary polarised brightness could significantly reduce the ambiguity in performing three-dimensional reconstruction of local electron density from white-light imaging.

Observational tracking of the 2D structure of coronal mass ejections between the sun and 1 AU

The Solar TErrestrial RElations Observatory (STEREO) provides high cadence and high resolution images of the structure and morphology of coronal mass ejections (CMEs) in the inner heliosphere. CME directions and propagation speeds have often been estimated through the use of time-elongation maps obtained from the STEREO Heliospheric Imager (HI) data. Many of these CMEs have been identified by citizen scientists working within the SolarStormWatch project ( www.solarstormwatch.com ) as they work towards providing robust real-time identification of Earth-directed CMEs. The wide field of view of HI allows scientists to directly observe the two-dimensional (2D) structures, while the relative simplicity of time-elongation analysis means that it can be easily applied to many such events, thereby enabling a much deeper understanding of how CMEs evolve between the Sun and the Earth. For events with certain orientations, both the rear and front edges of the CME can be monitored at varying heliocentric distances (R) between the Sun and 1 AU. Here we take four example events with measurable position angle widths and identified by the citizen scientists. These events were chosen for the clarity of their structure within the HI cameras and their long track lengths in the time-elongation maps. We show a linear dependency with R for the growth of the radial width (W) and the 2D aspect ratio (χ) of these CMEs, which are measured out to ≈ 0.7 AU. We estimated the radial width from a linear best fit for the average of the four CMEs. We obtained the relationships W=0.14R+0.04 for the width and χ=2.5R+0.86 for the aspect ratio (W and R in units of AU).