Differences between the CME fronts tracked by an expert, an automated algorithm, and the Solar Stormwatch project

Observations from the Heliospheric Imager (HI) instruments aboard the twin STEREO spacecraft have enabled the compilation of several catalogues of coronal mass ejections (CMEs), each characterizing the propagation of CMEs through the inner heliosphere. Three such catalogues are the Rutherford Appleton Laboratory (RAL)-HI event list, the Solar Stormwatch CME catalogue, and, presented here, the J-tracker catalogue. Each catalogue uses a different method to characterize the location of CME fronts in the HI images: manual identification by an expert, the statistical reduction of the manual identifications of many citizen scientists, and an automated algorithm. We provide a quantitative comparison of the differences between these catalogues and techniques, using 51 CMEs common to each catalogue. The time-elongation profiles of these CME fronts are compared, as are the estimates of the CME kinematics derived from application of three widely used single-spacecraft-fitting techniques. The J-tracker and RAL-HI profiles are most similar, while the Solar Stormwatch profiles display a small systematic offset. Evidence is presented that these differences arise because the RAL-HI and J-tracker profiles follow the sunward edge of CME density enhancements, while Solar Stormwatch profiles track closer to the antisunward (leading) edge. We demonstrate that the method used to produce the time-elongation profile typically introduces more variability into the kinematic estimates than differences between the various single-spacecraft-fitting techniques. This has implications for the repeatability and robustness of these types of analyses, arguably especially so in the context of space weather forecasting, where it could make the results strongly dependent on the methods used by the forecaster.

Tests of sunspot number sequences: 3. Effects of regression procedures on the calibration of historic sunspot data

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 RB 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 RB are then re-scaled to the full observed RGO group number RA using a variety of regression techniques. It is found that a very high correlation between RA and RB (rAB > 0.98) does not prevent large errors in the intercalibration (for example sunspot maximum values can be over 30 % too large even for such levels of rAB). In generating the backbone sunspot number (RBB), Svalgaard and Schatten (2015, this issue) 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.

A new calibrated sunspot group series since 1749: statistics of active day fractions

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.

Lightning as a space-weather hazard: UK thunderstorm activity modulated by the passage of the heliospheric current sheet

Lightning flash rates, RL, are modulated by corotating interaction regions (CIRs) and the polarity of the heliospheric magnetic field (HMF) in near-Earth space. As the HMF polarity reverses at the heliospheric current sheet (HCS), typically within a CIR, these phenomena are likely related. In this study, RL is found to be significantly enhanced at the HCS and at 27 days prior/after. The strength of the enhancement depends on the polarity of the HMF reversal at the HCS. Near-Earth solar and galactic energetic particle fluxes are also ordered by HMF polarity, though the variations qualitatively differ from RL, with the main increase occurring prior to the HCS crossing. Thus, the CIR effect on lightning is either the result of compression/amplification of the HMF (and its subsequent interaction with the terrestrial system) or that energetic particle preconditioning of the Earth system prior to the HMF polarity change is central to solar wind lightning coupling mechanism.

Do the legs of magnetic clouds contain twisted flux-rope magnetic fields?

Magnetic clouds (MCs) are a subset of interplanetary coronal mass ejections (ICMEs) characterised primarily by a smooth rotation in the magnetic field direction indicative of the presence of a magnetic flux rope. Energetic particle signatures suggest MC flux ropes remain magnetically connected to the Sun at both ends, leading to widely used model of global MC structure as an extended flux rope, with a loop-like axis stretching out from the Sun into the heliosphere and back to the Sun. The time of flight of energetic particles, however, suggests shorter magnetic field line lengths than such a continuous twisted flux rope would produce. In this study, two simple models are compared with observed flux rope axis orientations of 196 MCs to show that the flux rope structure is confined to the MC leading edge. The magnetic cloud “legs,” which magnetically connect the flux rope to the Sun, are not recognisable as MCs and thus are unlikely to contain twisted flux rope fields. Spacecraft encounters with these non-flux rope legs may provide an explanation for the frequent observation of non-magnetic cloud ICMEs.

I killed my father, I ate human flesh, I quiver with joy: an obsession with Pier Paolo Pasolini

Group show, curated by Invisible Exports Gallery. Featuring work by Michael Bilsborough, Lizzi Bougatsos, BREYER P-ORRIDGE, Asger Carlsen, Troels Carlsen, Walt Cassidy, Andy Coolquitt, Vaginal Davis, Carlton DeWoody, Joey Frank, Paul Gabrielli, Ludovica Gioscia, Luis Gispert, Terence Hannum, Karen Heagle, Timothy Hull, Doug Ischar, Brian Kenny, Jeremy Kost, Aaron Krach, Yeni Mao, Leigha Mason, Mark McCoy, Robert Melee, Lucas Michael, Jennifer Needleman, Brent Owens, Paul P., Paolo Di Paolo, Franklin Preston, John Russell, Xaviera Simmons, Duston Spear, Scott Treleaven, Ramon Vega, Jordan Wolfson, Dustin Yellin

Tests of sunspot number sequences: 1. Using ionosonde data

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. This test is carried out for the original composite of the Wolf/Zürich/International sunspot number [R], the new “backbone” group sunspot number [RBB] and the proposed “corrected sunspot number” [RC]. 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, R, RBB, and 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.

On the origins and timescales of geoeffective IMF

Southward Interplanetary Magnetic Field (IMF) in the Geocentric Solar Magnetospheric (GSM) reference frame is the key element that controls the level of space-weather disturbance in Earth’s magnetosphere, ionosphere and thermosphere. We discuss the relation of this geoeffective IMF component to the IMF in the Geocentric Solar Ecliptic (GSE) frame and, using the almost continuous interplanetary data for 1996-2015 (inclusive), we show that large geomagnetic storms are always associated with strong southward, out-of-ecliptic field in the GSE frame: dipole tilt effects, that cause the difference between the southward field in the GSM and GSE frames, generally make only a minor contribution to these strongest storms. The time-of-day/time-of-year response patterns of geomagnetic indices and the optimum solar wind coupling function are both influenced by the timescale of the index response. We also study the occurrence spectrum of large out-of-ecliptic field and show that for one-hour averages it is, surprisingly, almost identical in ICMEs (Interplanetary Coronal Mass Ejections), around CIRs/SIRs (Corotating and Stream Interaction Regions) and in the “quiet” solar wind (which is shown to be consistent with the effect of weak SIRs). However, differences emerge when the timescale over which the field remains southward is considered: for longer averaging timescales the spectrum is broader inside ICMEs, showing that these events generate longer intervals of strongly southward average IMF and consequently stronger geomagnetic storms. The behavior of out-of-ecliptic field with timescale is shown to be very similar to that of deviations from the predicted Parker spiral orientation, suggesting the two share common origins.