sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Gallagher, Peter T.

Normalized to: Gallagher, P.

89 article(s) in total. 391 co-authors, from 1 to 27 common article(s). Median position in authors list is 3,0.

[1]  [pdf] - 2083535
Statistical Study of GOES X-ray Quasi-Periodic Pulsations in Solar Flares
Comments: 17 pages, 13 figures, accepted in The Astrophysical Journal
Submitted: 2020-04-24
Small amplitude quasi-periodic pulsations (QPPs) detected in soft X-ray emission are commonplace in many flares. To date, the underpinning processes resulting in the QPPs are unknown. In this paper, we attempt to constrain the prevalence of \textit{stationary} QPPs in the largest statistical study to date, including a study of the relationship of QPP periods to the properties of the flaring active region, flare ribbons, and CME affiliation. We build upon the work of \cite{inglis2016} and use a model comparison test to search for significant power in the Fourier spectra of lightcurves of the GOES 1--8~\AA\ channel. We analyze all X-, M- and C- class flares of the past solar cycle, a total of 5519 flares, and search for periodicity in the 6-300~s timescale range. Approximately 46\% of X-class, 29\% of M-class and 7\% of C-class flares show evidence of stationary QPPs, with periods that follow a log-normal distribution peaked at 20~s. The QPP periods were found to be independent of flare magnitude, however a positive correlation was found between QPP period and flare duration. No dependence of the QPP periods to the global active region properties was identified. A positive correlation was found between QPPs and ribbon properties including unsigned magnetic flux, ribbon area and ribbon separation distance. We found that both flares with and without an associated CME can host QPPs. Furthermore, we demonstrate that for X- and M- class flares, decay phase QPPs have statistically longer periods than impulsive phase QPPs.
[2]  [pdf] - 2081172
A new facility for airborne solar astronomy: NASA's WB-57 at the 2017 total solar eclipse
Comments: 17 pages, 10 figures; accepted for publication by the Astrophysical Journal
Submitted: 2020-04-20
NASA's WB-57 High Altitude Research Program provides a deployable, mobile, stratospheric platform for scientific research. Airborne platforms are of particular value for making coronal observations during total solar eclipses because of their ability both to follow the Moon's shadow and to get above most of the atmospheric airmass that can interfere with astronomical observations. We used the 2017 Aug 21 eclipse as a pathfinding mission for high-altitude airborne solar astronomy, using the existing high-speed visible-light and near-/mid-wave infrared imaging suite mounted in the WB-57 nose cone. In this paper, we describe the aircraft, the instrument, and the 2017 mission; operations and data acquisition; and preliminary analysis of data quality from the existing instrument suite. We describe benefits and technical limitations of this platform for solar and other astronomical observations. We present a preliminary analysis of the visible-light data quality and discuss the limiting factors that must be overcome with future instrumentation. We conclude with a discussion of lessons learned from this pathfinding mission and prospects for future research at upcoming eclipses, as well as an evaluation of the capabilities of the WB-57 platform for future solar astronomy and general astronomical observation.
[3]  [pdf] - 2054142
A Comparison of Flare Forecasting Methods. IV. Evaluating Consecutive-Day Forecasting Patterns
Comments: 33 pages, 13 figures, accepted for publication in ApJ
Submitted: 2020-01-08, last modified: 2020-01-21
A crucial challenge to successful flare prediction is forecasting periods that transition between "flare-quiet" and "flare-active". Building on earlier studies in this series (Barnes et al. 2016; Leka et al. 2019a,b) in which we describe methodology, details, and results of flare forecasting comparison efforts, we focus here on patterns of forecast outcomes (success and failure) over multi-day periods. A novel analysis is developed to evaluate forecasting success in the context of catching the first event of flare-active periods, and conversely, of correctly predicting declining flare activity. We demonstrate these evaluation methods graphically and quantitatively as they provide both quick comparative evaluations and options for detailed analysis. For the testing interval 2016-2017, we determine the relative frequency distribution of two-day dichotomous forecast outcomes for three different event histories (i.e., event/event, no-event/event and event/no-event), and use it to highlight performance differences between forecasting methods. A trend is identified across all forecasting methods that a high/low forecast probability on day-1 remains high/low on day-2 even though flaring activity is transitioning. For M-class and larger flares, we find that explicitly including persistence or prior flare history in computing forecasts helps to improve overall forecast performance. It is also found that using magnetic/modern data leads to improvement in catching the first-event/first-no-event transitions. Finally, 15% of major (i.e., M-class or above) flare days over the testing interval were effectively missed due to a lack of observations from instruments away from the Earth-Sun line.
[4]  [pdf] - 2030663
Evolution of the Alfv\'en Mach number associated with a coronal mass ejection shock
Comments: 8 pages, 4 figures
Submitted: 2019-12-04, last modified: 2019-12-06
The Sun regularly produces large-scale eruptive events, such as coronal mass ejections (CMEs) that can drive shock waves through the solar corona. Such shocks can result in electron acceleration and subsequent radio emission in the form of a type II radio burst. However, the early-phase evolution of shock properties and its relationship to type II burst evolution is still subject to investigation. Here we study the evolution of a CME-driven shock by comparing three commonly used methods of calculating the Alfv\'en Mach number ($M_A$), namely: shock geometry, a comparison of CME speed to a model of the coronal Alfv\'en speed, and the type II band-splitting method. We applied the three methods to the 2017 September 2 event, focusing on the shock wave observed in extreme ultraviolet (EUV) by the Solar Ultraviolet Imager (SUVI) on board GOES-16, in white-light by the Large Angle and Spectrometric Coronagraph (LASCO) on board SOHO, and the type II radio burst observed by the Irish Low Frequency Array (I-LOFAR). We show that the three different methods of estimating shock $M_A$ yield consistent results and provide a means of relating shock property evolution to the type II emission duration. The type II radio emission emerged from near the nose of the CME when $M_A$ was in the range 1.4-2.4 at a heliocentric distance of $\sim$1.6 $R_\odot$. The emission ceased when the CME nose reached $\sim$2.4 $R_\odot$, despite an increasing Alfv\'en Mach number (up to 4). We suggest the radio emission cessation is due to the lack of quasi-perpendicular geometry at this altitude, which inhibits efficient electron acceleration and subsequent radio emission.
[5]  [pdf] - 1951626
Multiple Regions of Shock-accelerated Particles during a Solar Coronal Mass Ejection
Comments: 31 pages, 6 figures
Submitted: 2019-08-30
The Sun is an active star that can launch large eruptions of magnetised plasma into the heliosphere, called coronal mass ejections (CMEs). These ejections can drive shocks that accelerate particles to high energies, often resulting in radio emission at low frequencies (<200 MHz). To date, the relationship between the expansion of CMEs, shocks and particle acceleration is not well understood, partly due to the lack of radio imaging at low frequencies during the onset of shock-producing CMEs. Here, we report multi-instrument radio, white-light and ultraviolet imaging of the second largest flare in Solar Cycle 24 (2008-present) and its associated fast CME (3038+/-288 km/s). We identify the location of a multitude of radio shock signatures, called herringbones, and find evidence for shock accelerated electron beams at multiple locations along the expanding CME. These observations support theories of non-uniform, rippled shock fronts driven by an expanding CME in the solar corona.
[6]  [pdf] - 1928169
The Breakthrough Listen Search for Extraterrestrial Intelligence
Comments: Astronomy and Astrophysics Decadal Survey 2020, APC white paper, To be published in BAAS
Submitted: 2019-07-11, last modified: 2019-08-02
The discovery of the ubiquity of habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, have ushered in a renaissance in the millenia-old quest to answer our most profound question about the Universe and our place within it - Are we alone? The Breakthrough Listen Initiative, announced in July 2015 as a 10-year 100M USD program, is the most comprehensive effort in history to quantify the distribution of advanced, technologically capable life in the universe. In this white paper, we outline the status of the on-going observing campaign with our primary observing facilities, as well as planned activities with these instruments over the next few years. We also list collaborative facilities which will conduct searches for technosignatures in either primary observing mode, or commensally. We highlight some of the novel analysis techniques we are bringing to bear on multi-petabyte data sets, including machine learning tools we are deploying to search for a broader range of technosignatures than was previously possible.
[7]  [pdf] - 1953630
A Comparison of Flare Forecasting Methods. III. Systematic Behaviors of Operational Solar Flare Forecasting Systems
Comments: 23 pages, 6 figures, accepted for publication in The Astrophysical Journal
Submitted: 2019-07-05
A workshop was recently held at Nagoya University (31 October - 02 November 2017), sponsored by the Center for International Collaborative Research, at the Institute for Space-Earth Environmental Research, Nagoya University, Japan, to quantitatively compare the performance of today's operational solar flare forecasting facilities. Building upon Paper I of this series (Barnes et al. 2016), in Paper II (Leka et al. 2019) we described the participating methods for this latest comparison effort, the evaluation methodology, and presented quantitative comparisons. In this paper we focus on the behavior and performance of the methods when evaluated in the context of broad implementation differences. Acknowledging the short testing interval available and the small number of methods available, we do find that forecast performance: 1) appears to improve by including persistence or prior flare activity, region evolution, and a human "forecaster in the loop"; 2) is hurt by restricting data to disk-center observations; 3) may benefit from long-term statistics, but mostly when then combined with modern data sources and statistical approaches. These trends are arguably weak and must be viewed with numerous caveats, as discussed both here and in Paper II. Following this present work, we present in Paper IV a novel analysis method to evaluate temporal patterns of forecasting errors of both types (i.e., misses and false alarms; Park et al. 2019). Hence, most importantly, with this series of papers we demonstrate the techniques for facilitating comparisons in the interest of establishing performance-positive methodologies.
[8]  [pdf] - 1953629
A Comparison of Flare Forecasting Methods. II. Benchmarks, Metrics and Performance Results for Operational Solar Flare Forecasting Systems
Comments: 26 pages, 5 figures, accepted for publication in the Astrophysical Journal Supplement Series
Submitted: 2019-07-05
Solar flares are extremely energetic phenomena in our Solar System. Their impulsive, often drastic radiative increases, in particular at short wavelengths, bring immediate impacts that motivate solar physics and space weather research to understand solar flares to the point of being able to forecast them. As data and algorithms improve dramatically, questions must be asked concerning how well the forecasting performs; crucially, we must ask how to rigorously measure performance in order to critically gauge any improvements. Building upon earlier-developed methodology (Barnes et al, 2016, Paper I), international representatives of regional warning centers and research facilities assembled in 2017 at the Institute for Space-Earth Environmental Research, Nagoya University, Japan to - for the first time - directly compare the performance of operational solar flare forecasting methods. Multiple quantitative evaluation metrics are employed, with focus and discussion on evaluation methodologies given the restrictions of operational forecasting. Numerous methods performed consistently above the "no skill" level, although which method scored top marks is decisively a function of flare event definition and the metric used; there was no single winner. Following in this paper series we ask why the performances differ by examining implementation details (Leka et al. 2019, Paper III), and then we present a novel analysis method to evaluate temporal patterns of forecasting errors in (Park et al. 2019, Paper IV). With these works, this team presents a well-defined and robust methodology for evaluating solar flare forecasting methods in both research and operational frameworks, and today's performance benchmarks against which improvements and new methods may be compared.
[9]  [pdf] - 1888315
Properties and magnetic origins of solar S-bursts
Comments: Added acknowledgement to the Solar team from the Institute of Radio Astronomy of the National Academy of Sciences of Ukraine
Submitted: 2019-01-22, last modified: 2019-05-23
Context. Solar activity is often accompanied by solar radio emission, consisting of numerous types of solar radio bursts. At low frequencies (<100 MHz) radio bursts with short durations of milliseconds, such as solar S-bursts, have been identified. To date, their origin and many of their characteristics remain unclear. Aims. We report observations from the Ukrainian T-shaped Radio telescope, (UTR-2), and the LOw Frequency ARray (LOFAR) which give us new insight into their nature. Methods. Over 3000 S-bursts were observed on 9 July 2013 at frequencies of 17.4-83.1 MHz during a period of low solar activity. Leading models of S-burst generation were tested by analysing the spectral properties of S-bursts and estimating coronal magnetic field strengths. Results. S-bursts were found to have short durations of 0.5-0.9 s. Multiple instruments were used to measure the dependence of drift rate on frequency which is represented by a power law with an index of 1.57. For the first time, we show a linear relation between instantaneous bandwidth and frequency over a wide frequency band. The flux calibration and high sensitivity of UTR-2 enabled measurements of their fluxes, which yielded 11+3 SFU. The source particle velocities of S-bursts were found to be ~0.07 c. S-burst source heights were found to range from 1.3 R_Sun to 2 R_Sun. Furthermore, a contemporary theoretical model of S-burst generation was used to conduct remote sensing of the coronal magnetic field at these heights which yielded values of 0.9-5.8 G. Within error, these values are comparable to those predicted by various relations between magnetic field strength and height in the corona.
[10]  [pdf] - 1868126
Persistent Quasi-Periodic Pulsations During a Large X-Class Solar Flare
Comments: Accepted for publication in ApJ - 13 pages, 9 figures, 1 accompanying movie found at -
Submitted: 2019-03-04
Solar flares often display pulsating and oscillatory signatures in the emission, known as quasi-periodic pulsations (QPP). QPP are typically identified during the impulsive phase of flares, yet in some cases, their presence is detected late into the decay phase. Here, we report extensive fine structure QPP that are detected throughout the large X8.2 flare from 2017 September 10. Following the analysis of the thermal pulsations observed in the GOES/XRS and the 131 A channel of SDO/AIA, we find a pulsation period of ~65 s during the impulsive phase followed by lower amplitude QPP with a period of ~150 s in the decay phase, up to three hours after the peak of the flare. We find that during the time of the impulsive QPP, the soft X-ray source observed with RHESSI rapidly rises at a velocity of approximately 17 km/s following the plasmoid/coronal mass ejection (CME) eruption. We interpret these QPP in terms of a manifestation of the reconnection dynamics in the eruptive event. During the long-duration decay phase lasting several hours, extended downward contractions of collapsing loops/plasmoids that reach the top of the flare arcade are observed in EUV. We note that the existence of persistent QPP into the decay phase of this flare are most likely related to these features. The QPP during this phase are discussed in terms of MHD wave modes triggered in the post-flaring loops.
[11]  [pdf] - 1795933
Expansion of High Speed Solar Wind Streams from Coronal Holes through the Inner Heliosphere
Comments: 8 pages, 4 figures, Accepted for publication in The Astrophysical Journal Letters (APJL) on 25th November 2018
Submitted: 2018-11-28, last modified: 2018-11-29
Coronal holes (CHs) are regions of open magnetic flux which are the source of high speed solar wind (HSSW) streams. To date, it is not clear which aspects of CHs are of most influence on the properties of the solar wind as it expands through the Heliosphere. Here, we study the relationship between CH properties extracted from AIA (Atmospheric Imaging Assembly) images using CHIMERA (Coronal Hole Identification via Multi-thermal Emission Recognition Algorithm) and HSSW measurements from ACE (Advanced Composition Explorer) at L1. For CH longitudinal widths $\Delta\theta_{CH}<$67$^{\circ}$, the peak SW velocity ($v_{max}$) is found to scale as $v_{max}~\approx~330.8~+~5.7~\Delta\theta_{CH}$~km~s$^{-1}$. For larger longitudinal widths ($\Delta\theta_{CH}>$67$^{\circ}$), $v_{max}$ is found to tend to a constant value ($\sim$710~km~s$^{-1}$). Furthermore, we find that the duration of HSSW streams ($\Delta t$) are directly related to the longitudinal width of CHs ($\Delta t_{SW}$~$\approx$~0.09$\Delta\theta_{CH}$) and that their longitudinal expansion factor is $f_{SW}~\approx 1.2~\pm 0.1$. We also derive an expression for the coronal hole flux-tube expansion factor, $f_{FT}$, which varies as $f_{SW} \gtrsim f_{FT} \gtrsim 0.8$. These results enable us to estimate the peak speeds and durations of HSSW streams at L1 using the properties of CHs identified in the solar corona.
[12]  [pdf] - 1774049
Photospheric Shear Flows in Solar Active Regions and Their Relation to Flare Occurrence
Comments: 19 pages, 8 figures, accepted for publication in Solar Physics
Submitted: 2018-07-20
Solar active regions (ARs) that produce major flares typically exhibit strong plasma shear flows around photospheric magnetic polarity inversion lines (MPILs). It is therefore important to quantitatively measure such photospheric shear flows in ARs for a better understanding of their relation to flare occurrence. Photospheric flow fields were determined by applying the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM) method to a large data set of 2,548 co-aligned pairs of AR vector magnetograms with 12-min separation over the period 2012-2016. From each AR flow-field map, three shear-flow parameters were derived corresponding to the mean (<S>), maximum (S_max) and integral (S_sum) shear-flow speeds along strong-gradient, strong-field MPIL segments. We calculated flaring rates within 24 hr as a function of each shear-flow parameter, and also investigated the relation between the parameters and the waiting time ({\tau}) until the next major flare (class M1.0 or above) after the parameter observation. In general, it is found that the larger S_sum an AR has, the more likely it is for the AR to produce flares within 24 hr. It is also found that among ARs which produce major flares, if one has a larger value of S_sum then {\tau} generally gets shorter. These results suggest that large ARs with widespread and/or strong shear flows along MPILs tend to not only be more flare productive, but also produce major flares within 24 hr or less.
[13]  [pdf] - 1747966
Three-dimensional magnetic reconnection in a collapsing coronal loop system
Comments: 7 pages, 4 figure
Submitted: 2018-06-25
Magnetic reconnection is believed to be the primary mechanism by which non-potential energy stored in coronal magnetic fields is rapidly released during solar eruptive events. Unfortunately, owing to the small spatial scales on which reconnection is thought to occur, it is not directly observable in the solar corona. However, larger scale processes, such as associated inflow and outflow, and signatures of accelerated particles have been put forward as evidence of reconnection. We explore the origin of a persistent Type I radio source that accompanies a coronal X--shaped structure during its passage across the disk. Of particular interest is the time range around a partial collapse of the structure that is associated with inflow, outflow, and signatures of particle acceleration. Imaging radio observations from the NRH were used to localise the radio source. SDO AIA EUV observations from the same time period were analysed, looking for evidence of inflows and outflows. Further \texttt{mpole} magnetic reconstructions using SDO HMI observations allowed the magnetic connectivity associated with the radio source to be determined. The Type I radio source was well aligned with a magnetic separator identified in the extrapolations. During the partial collapse, gradual (1 km/s) and fast (5 km/s) inflow phases and fast (30 km/s) and rapid (80-100 km/s) outflow phases were observed, resulting in an estimated reconnection rate of ~ 0.06. The radio source brightening and dimming was found to be co-temporal with increased soft x-ray emission in both RHESSI and GOES. We interpret the brightening and dimming of the radio emission as evidence for accelerated electrons in the reconnection region responding to a gradual fall and rapid rise in electric drift velocity, in response to the inflowing and outflowing field lines. These results present a comprehensive example of 3D null-point reconnection.
[14]  [pdf] - 1676429
Flare Forecasting Using the Evolution of McIntosh Sunspot Classifications
Comments: 21 pages, 9 figures
Submitted: 2018-05-02
Most solar flares originate in sunspot groups, where magnetic field changes lead to energy build-up and release. However, few flare-forecasting methods use information of sunspot-group evolution, instead focusing on static point-in-time observations. Here, a new forecast method is presented based upon the 24-hr evolution in McIntosh classification of sunspot groups. Evolution-dependent $\geqslant$C1.0 and $\geqslant$M1.0 flaring rates are found from NOAA-numbered sunspot groups over December 1988 to June 1996 (Solar Cycle 22; SC22) before converting to probabilities assuming Poisson statistics. These flaring probabilities are used to generate operational forecasts for sunspot groups over July 1996 to December 2008 (SC23), with performance studied by verification metrics. Major findings are: i) considering Brier skill score (BSS) for $\geqslant$C1.0 flares, the evolution-dependent McIntosh-Poisson method ($\text{BSS}_{\text{evolution}}=0.09$) performs better than the static McIntosh-Poisson method ($\text{BSS}_{\text{static}} = -0.09$); ii) low BSS values arise partly from both methods over-forecasting SC23 flares from the SC22 rates, symptomatic of $\geqslant$C1.0 rates in SC23 being on average $\approx$80% of those in SC22 (with $\geqslant$M1.0 being $\approx$50%); iii) applying a bias-correction factor to reduce the SC22 rates used in forecasting SC23 flares yields modest improvement in skill relative to climatology for both methods ($\mathrm{BSS}^{\mathrm{corr}}_{\mathrm{static}} = 0.09$ and $\mathrm{BSS}^{\mathrm{corr}}_{\mathrm{evolution}} = 0.20$) and improved forecast reliability diagrams.
[15]  [pdf] - 1682652
CMEs in the Heliosphere: I. A Statistical Analysis of the Observational Properties of CMEs Detected in the Heliosphere from 2007 to 2017 by STEREO/HI-1
Comments: 26 pages, 8 figures, in press at Solar Physics. The abstract has been shortened to comply with the allowed length at
Submitted: 2018-04-06, last modified: 2018-04-09
We present a statistical analysis of coronal mass ejections (CMEs) imaged by the Heliospheric Imager (HI) instruments aboard NASAs twin-spacecraft STEREO mission between April 2007 and August 2017 for STEREO-A and between April 2007 and September 2014 for STEREO-B. The analysis exploits a catalogue that was generated within the FP7 HELCATS project. Here, we focus on the observational characteristics of CMEs imaged in the heliosphere by the inner (HI-1) cameras. More specifically, in this paper we present distributions of the basic observational parameters - namely occurrence frequency, central position angle (PA) and PA span - derived from nearly 2000 detections of CMEs in the heliosphere by HI-1 on STEREO-A or STEREO-B from the minimum between Solar Cycles 23 and 24 to the maximum of Cycle 24; STEREO-A analysis includes a further 158 CME detections from the descending phase of Cycle 24, by which time communication with STEREO-B had been lost. We compare heliospheric CME characteristics with properties of CMEs observed at coronal altitudes, and with sunspot number. As expected, heliospheric CME rates correlate with sunspot number, and are not inconsistent with coronal rates once instrumental factors/differences in cataloguing philosophy are considered. As well as being more abundant, heliospheric CMEs, like their coronal counterparts, tend to be wider during solar maximum. Our results confirm previous coronagraph analyses suggesting that CME launch sites don't simply migrate to higher latitudes with increasing solar activity. At solar minimum, CMEs tend to be launched from equatorial latitudes while, at maximum, CMEs appear to be launched over a much wider latitude range; this has implications for understanding the CME/solar source association. Our analysis provides some supporting evidence for the systematic dragging of CMEs to lower latitude as they propagate outwards.
[16]  [pdf] - 1721063
Shock location and CME 3D reconstruction of a solar type II radio burst with LOFAR
Comments: 7 Figures, 9 Pages
Submitted: 2018-04-03
Type II radio bursts are evidence of shocks in the solar atmosphere and inner heliosphere that emit radio waves ranging from sub-meter to kilometer lengths. These shocks may be associated with CMEs and reach speeds higher than the local magnetosonic speed. Radio imaging of decameter wavelengths (20-90 MHz) is now possible with LOFAR, opening a new radio window in which to study coronal shocks that leave the inner solar corona and enter the interplanetary medium and to understand their association with CMEs. To this end, we study a coronal shock associated with a CME and type II radio burst to determine the locations at which the radio emission is generated, and we investigate the origin of the band-splitting phenomenon.
[17]  [pdf] - 1661814
Connecting Coronal Mass Ejections to their Solar Active Region Sources: Combining Results from the HELCATS and FLARECAST Projects
Comments: Accepted for publication in Solar Physics. 31 pages, 15 figures, 4 tables
Submitted: 2018-03-17
Coronal mass ejections (CMEs) and other solar eruptive phenomena can be physically linked by combining data from a multitude of ground-based and space-based instruments alongside models, however this can be challenging for automated operational systems. The EU Framework Package 7 HELCATS project provides catalogues of CME observations and properties from the Helio- spheric Imagers onboard the two NASA/STEREO spacecraft in order to track the evolution of CMEs in the inner heliosphere. From the main HICAT catalogue of over 2,000 CME detections, an automated algorithm has been developed to connect the CMEs observed by STEREO to any corresponding solar flares and active region (AR) sources on the solar surface. CME kinematic properties, such as speed and angular width, are compared with AR magnetic field properties, such as magnetic flux, area, and neutral line characteristics. The resulting LOWCAT catalogue is also compared to the extensive AR property database created by the EU Horizon 2020 FLARECAST project, which provides more complex magnetic field parameters derived from vector magnetograms. Initial statistical analysis has been undertaken on the new data to provide insight into the link between flare and CME events, and characteristics of eruptive ARs. Warning thresholds determined from analysis of the evolution of these parameters is shown to be a useful output for operational space weather purposes. Parameters of particular interest for further analysis include total unsigned flux, vertical current, and current helicity. The automated method developed to create the LOWCAT catalogue may also be useful for future efforts to develop operational CME forecasting.
[18]  [pdf] - 1698075
LOFAR observations of the quiet solar corona
Submitted: 2018-03-01
The quiet solar corona emits meter-wave thermal bremsstrahlung. Coronal radio emission can only propagate above that radius, $R_\omega$, where the local plasma frequency eqals the observing frequency. The radio interferometer LOw Frequency ARray (LOFAR) observes in its low band (10 -- 90 MHz) solar radio emission originating from the middle and upper corona. We present the first solar aperture synthesis imaging observations in the low band of LOFAR in 12 frequencies each separated by 5 MHz. From each of these radio maps we infer $R_\omega$, and a scale height temperature, $T$. These results can be combined into coronal density and temperature profiles. We derived radial intensity profiles from the radio images. We focus on polar directions with simpler, radial magnetic field structure. Intensity profiles were modeled by ray-tracing simulations, following wave paths through the refractive solar corona, and including free-free emission and absorption. We fitted model profiles to observations with $R_\omega$ and $T$ as fitting parameters. In the low corona, $R_\omega < 1.5$ solar radii, we find high scale height temperatures up to 2.2e6 K, much more than the brightness temperatures usually found there. But if all $R_\omega$ values are combined into a density profile, this profile can be fitted by a hydrostatic model with the same temperature, thereby confirming this with two independent methods. The density profile deviates from the hydrostatic model above 1.5 solar radii, indicating the transition into the solar wind. These results demonstrate what information can be gleaned from solar low-frequency radio images. The scale height temperatures we find are not only higher than brightness temperatures, but also than temperatures derived from coronograph or EUV data. Future observations will provide continuous frequency coverage, eliminating the need for local hydrostatic density models.
[19]  [pdf] - 1867942
Characteristics of type III radio bursts and solar S bursts
Comments: 10 pages, 4 figures. arXiv admin note: text overlap with arXiv:1507.07496
Submitted: 2018-02-27
The Sun is an active source of radio emission which is often associated with the acceleration of electrons arising from processes such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), numerous solar S bursts (where S stands for short) and storms of Type III radio bursts have been observed, that are not directly relates to flares and CMEs. Here, we expand our understanding on the spectral characteristic of these two different types of radio bursts based on observations from the Low Frequency Array (LOFAR). On 9 July 2013, over 3000 solar S bursts accompanied by over 800 Type III radio bursts were observed over a time period of ~8 hours. The characteristics of Type III radio bursts are consistent to previous studies, while S bursts show narrow bandwidths, durations and drift rates of about 1/2 the drift rate of Type III bursts. Type III bursts and solar S bursts occur in a region in the corona where plasma emission is the dominant emission mechanism as determined by data constrained density and magnetic field models.
[20]  [pdf] - 1622392
Active Region Photospheric Magnetic Properties Derived from Line-of-sight and Radial Fields
Submitted: 2017-12-19
The effect of using two representations of the normal-to-surface magnetic field to calculate photospheric measures that are related to active region (AR) potential for flaring is presented. Several AR properties were computed using line-of-sight ($B_{\rm los}$) and spherical-radial ($B_{r}$) magnetograms from the Space-weather HMI Active Region Patch (SHARP) products of the \emph{Solar Dynamics Observatory}, characterizing the presence and features of magnetic polarity inversion lines, fractality, and magnetic connectivity of the AR photospheric field. The data analyzed corresponds to $\approx$ 4,000 AR observations, achieved by randomly selecting 25 % of days between September 2012 and May 2016 for analysis at 6-hr cadence. Results from this statistical study include: i) the $B_{r}$ component results in a slight upwards shift of property values in a manner consistent with a field-strength underestimation by the $B_{\rm los}$ component; ii) using the $B_{r}$ component results in significantly lower inter-property correlation in one-third of the cases, implying more independent information about the state of the AR photospheric magnetic field; iii) flaring rates for each property vary between the field components in a manner consistent with the differences in property-value ranges resulting from the components; iv) flaring rates generally increase for higher values of properties, except Fourier spectral power index that has flare rates peaking around a value of $5/3$. These findings indicate that there may be advantages in using $B_{r}$ rather than $B_{\rm los}$ in calculating flare-related AR magnetic properties, especially for regions located far from central meridian.
[21]  [pdf] - 1849132
Automated Coronal Hole Identification via Multi-Thermal Intensity Segmentation
Submitted: 2017-11-30
Coronal holes (CH) are regions of open magnetic fields that appear as dark areas in the solar corona due to their low density and temperature compared to the surrounding quiet corona. To date, accurate identification and segmentation of CHs has been a difficult task due to their comparable intensity to local quiet Sun regions. Current segmentation methods typically rely on the use of single EUV passband and magnetogram images to extract CH information. Here, the Coronal Hole Identification via Multi-thermal Emission Recognition Algorithm (CHIMERA) is described, which analyses multi-thermal images from the Atmospheric Image Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) to segment coronal hole boundaries by their intensity ratio across three passbands (171 \AA, 193 \AA, and 211 \AA). The algorithm allows accurate extraction of CH boundaries and many of their properties, such as area, position, latitudinal and longitudinal width, and magnetic polarity of segmented CHs. From these properties, a clear linear relationship was identified between the duration of geomagnetic storms and coronal hole areas. CHIMERA can therefore form the basis of more accurate forecasting of the start and duration of geomagnetic storms.
[22]  [pdf] - 1589291
Pulsations in the Earth's Lower Ionosphere Synchronized with Solar Flare Emission
Comments: Accepted in JGR Space Physics
Submitted: 2017-10-04
Solar flare emission at X-ray and extreme ultraviolet (EUV) energies can cause substantial enhancements in the electron density in the Earth's lower ionosphere. It is now become clear that flares exhibit quasi-periodic pulsations with timescales of minutes at X-ray energies, but to date, it has not been known if the ionosphere is sensitive to this variability. Here, using a combination of Very Low Frequency (24 kHz) measurement together with space-based X-ray and EUV observations, we report pulsations of the ionospheric D-region, which are synchronized with a set of pulsating flare loops. Modeling of the ionosphere show that the D-region electron density varies by up to an order of magnitude over the timescale of the pulsations ($\sim$20 mins). Our results reveal that the Earth's ionosphere is more sensitive to small-scale changes in solar soft X-ray flux than previously thought, and implies that planetary ionospheres are closely coupled to small-scale changes in solar/stellar activity.
[23]  [pdf] - 1585811
The Association of a J-burst with a Solar Jet
Comments: 11 pages, 8 figures
Submitted: 2017-07-11, last modified: 2017-08-14
Context. The Sun is an active star that produces large-scale energetic events such as solar flares and coronal mass ejections and numerous smaller-scale events such as solar jets. These events are often associated with accelerated particles that can cause emission at radio wavelengths. The reconfiguration of the solar magnetic field in the corona is believed to be the cause of the majority of solar energetic events and accelerated particles. Aims. Here, we investigate a bright J-burst that was associated with a solar jet and the possible emission mechanism causing these two phenomena. Methods. We used data from the Solar Dynamics Observatory (SDO) to observe a solar jet, and radio data from the Low Frequency Array (LOFAR) and the Nan\c{c}ay Radioheliograph (NRH) to observe a J-burst over a broad frequency range (33-173 MHz) on 9 July 2013 at ~11:06 UT. Results. The J-burst showed fundamental and harmonic components and it was associated with a solar jet observed at extreme ultraviolet wavelengths with SDO. The solar jet occurred at a time and location coincident with the radio burst, in the northern hemisphere, and not inside a group of complex active regions in the southern hemisphere. The jet occurred in the negative polarity region of an area of bipolar plage. Newly emerged positive flux in this region appeared to be the trigger of the jet. Conclusions. Magnetic reconnection between the overlying coronal field lines and the newly emerged positive field lines is most likely the cause of the solar jet. Radio imaging provides a clear association between the jet and the J-burst which shows the path of the accelerated electrons.
[24]  [pdf] - 1581719
Modeling observations of solar coronal mass ejections with heliospheric imagers verified with the Heliophysics System Observatory
Comments: 27 pages, 7 figures, 3 tables, accepted for publication in the AGU journal "Space Weather" on 2017 July 1
Submitted: 2017-03-02, last modified: 2017-07-07
We present an advance towards accurately predicting the arrivals of coronal mass ejections (CMEs) at the terrestrial planets, including Earth. For the first time, we are able to assess a CME prediction model using data over 2/3 of a solar cycle of observations with the Heliophysics System Observatory. We validate modeling results of 1337 CMEs observed with the Solar Terrestrial Relations Observatory (STEREO) heliospheric imagers (HI) (science data) from 8 years of observations by 5 in situ observing spacecraft. We use the self-similar expansion model for CME fronts assuming 60 degree longitudinal width, constant speed and constant propagation direction. With these assumptions we find that 23%-35% of all CMEs that were predicted to hit a certain spacecraft lead to clear in situ signatures, so that for 1 correct prediction, 2 to 3 false alarms would have been issued. In addition, we find that the prediction accuracy does not degrade with the HI longitudinal separation from Earth. Predicted arrival times are on average within 2.6 +/- 16.6 hours difference of the in situ arrival time, similar to analytical and numerical modeling, and a true skill statistic of 0.21. We also discuss various factors that may improve the accuracy of space weather forecasting using wide-angle heliospheric imager observations. These results form a first order approximated baseline of the prediction accuracy that is possible with HI and other methods used for data by an operational space weather mission at the Sun-Earth L5 point.
[25]  [pdf] - 1584553
Detection and Interpretation Of Long-Lived X-Ray Quasi-Periodic Pulsations in the X-Class Solar Flare On 2013 May 14
Submitted: 2017-06-12
Quasi-periodic pulsations (QPP) seen in the time derivative of the GOES soft X-ray light curves are analyzed for the near-limb X3.2 event on 14 May 2013. The pulsations are apparent for a total of at least two hours from the impulsive phase to well into the decay phase, with a total of 163 distinct pulses evident to the naked eye. A wavelet analysis shows that the characteristic time scale of these pulsations increases systematically from $\sim$25 s at 01:10 UT, the time of the GOES peak, to $\sim$100 s at 02:00 UT. A second ridge in the wavelet power spectrum, most likely associated with flaring emission from a different active region, shows an increase from $\sim$40 s at 01:40 UT to $\sim$100 s at 03:10 UT. We assume that the QPP that produced the first ridge result from vertical kink-mode oscillations of the newly formed loops following magnetic reconnection in the coronal current sheet. This allows us to estimate the magnetic field strength as a function of altitude given the density, loop length, and QPP time scale as functions of time determined from the GOES light curves and RHESSI images. The calculated magnetic field strength of the newly formed loops ranges from about $\sim$500 G at an altitude of 24 Mm to a low value of $\sim$10 G at 60 Mm, in general agreement with the expected values at these altitudes. Fast sausage mode oscillations are also discussed and cannot be ruled out as an alternate mechanism for producing the QPP.
[26]  [pdf] - 1522672
Geomagnetically Induced Currents in the Irish Power Network during Geomagnetic Storms
Comments: 14 pages, 11 Figures, 4 Tables
Submitted: 2016-11-23
Geomagnetically induced currents (GICs) are a well-known terrestrial space weather hazard. They occur in power transmission networks and are known to have adverse effects in both high and mid-latitude countries. Here, we study GICs in the Irish power transmission network (geomagnetic latitude 54.7--58.5$^{\circ}$ N) during five geomagnetic storms (06-07 March 2016, 20-21 December 2015, 17-18 March 2015, 29-31 October 2003 and 13-14 March 1989). We simulate electric fields using a plane wave method together with two ground resistivity models, one of which is derived from magnetotelluric measurements (MT model). We then calculate GICs in the 220, 275 and 400~kV transmission network. During the largest of the storm periods studied, the peak electric field was calculated to be as large as 3.8~V~km\textsuperscript{-1}, with associated GICs of up to 23~A using our MT model. Using our homogenous resistivity model, those peak values were 1.46~V~km\textsuperscript{-1} and 25.8~A. We find that three 400 and 275~kV substations are the most likely locations for the Irish transformers to experience large GICs.
[27]  [pdf] - 1532674
Understanding the Physical Nature of Coronal "EIT Waves"
Comments: 26 pages, 2 figures, accepted for publication in Solar Physics
Submitted: 2016-11-16
For almost 20 years the physical nature of globally propagating waves in the solar corona (commonly called "EIT waves") has been controversial and subject to debate. Additional theories have been proposed over the years to explain observations that did not fit with the originally proposed fast-mode wave interpretation. However, the incompatibility of observations made using the Extreme-ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory with the fast-mode wave interpretation was challenged by differing viewpoints from the twin Solar Terrestrial Relations Observatory spacecraft and higher spatial/temporal resolution data from the Solar Dynamics Observatory. In this article, we reexamine the theories proposed to explain "EIT waves" to identify measurable properties and behaviours that can be compared to current and future observations. Most of us conclude that "EIT waves" are best described as fast-mode large-amplitude waves/shocks that are initially driven by the impulsive expansion of an erupting coronal mass ejection in the low corona.
[28]  [pdf] - 1532314
Observing the formation of flare-driven coronal rain
Comments: 23 pages, 21 figures
Submitted: 2016-10-28
Flare-driven coronal rain can manifest from rapidly cooled plasma condensations near coronal loop-tops in thermally unstable post-flare arcades. We detect 5 phases that characterise the post-flare decay: heating, evaporation, conductive cooling dominance for ~120 s, radiative / enthalpy cooling dominance for ~4700 s and finally catastrophic cooling occurring within 35-124 s leading to rain strands with s periodicity of 55-70 s. We find an excellent agreement between the observations and model predictions of the dominant cooling timescales and the onset of catastrophic cooling. At the rain formation site we detect co-moving, multi-thermal rain clumps that undergo catastrophic cooling from ~1 MK to ~22000 K. During catastrophic cooling the plasma cools at a maximum rate of 22700 K s-1 in multiple loop-top sources. We calculated the density of the EUV plasma from the DEM of the multi-thermal source employing regularised inversion. Assuming a pressure balance, we estimate the density of the chromospheric component of rain to be 9.21x10^11 +-1.76x10^11 cm-3 which is comparable with quiescent coronal rain densities. With up to 8 parallel strands in the EUV loop cross section, we calculate the mass loss rate from the post-flare arcade to be as much as 1.98x10^12 +/-4.95x10^11 g s-1. Finally, we reveal a close proximity between the model predictions of 10^5.8 K and the observed properties between 10^5.9 K and 10^6.2 K, that defines the temperature onset of catastrophic cooling. The close correspondence between the observations and numerical models suggests that indeed acoustic waves (with a sound travel time of 68 s) could play an important role in redistributing energy and sustaining the enthalpy-based radiative cooling.
[29]  [pdf] - 1532204
A large-scale search for evidence of quasi-periodic pulsations in solar flares
Comments: 16 pages, 8 figures, 4 tables. Accepted for publication in The Astrophysical Journal
Submitted: 2016-10-24
The nature of quasi-periodic pulsations in solar flares is poorly constrained, and critically the general prevalence of such signals in solar flares is unknown. Therefore, we perform a large-scale search for evidence of signals consistent with quasi-periodic pulsations in solar flares, focusing on the 1 - 300s timescale. We analyse 675 M- and X-class flares observed by GOES in 1-8\AA\ soft X-rays between 2011 February 1 and 2015 December 31. Additionally, over the same era we analyse Fermi/GBM 15-25 keV X-ray data for each of these flares that was associated with a Fermi/GBM solar flare trigger, a total of 261 events. Using a model comparison method, we determine whether there is evidence for a substantial enhancement in the Fourier power spectrum that may be consistent with a QPP signature, based on three tested models; a power-law plus a constant, a broken power-law plus constant, and a power-law-plus-constant with an additional QPP signature component. From this, we determine that ~30% of GOES events and ~8% of Fermi/GBM events show strong signatures consistent with classical interpretations of QPP. For the remaining events either two or more tested models cannot be strongly distinguished from each other, or the events are well-described by single power-law or broken power-law Fourier power spectra. For both instruments, a preferred characteristic timescale of ~5-30 s was found in the QPP-like events, with no dependence on flare magnitude in either GOES or GBM data. We also show that individual events in the sample show similar characteristic timescales in both GBM and GOES datasets. We discuss the implications of these results for our understanding of solar flares and possible QPP mechanisms.
[30]  [pdf] - 1531479
Radio Diagnostics of Electron Acceleration Sites During the Eruption of a Flux Rope in the Solar Corona
Submitted: 2016-09-06
Electron acceleration in the solar corona is often associated with flares and the eruption of twisted magnetic structures known as flux ropes. However, the locations and mechanisms of such particle acceleration during the flare and eruption are still subject to much investigation. Observing the exact sites of particle acceleration can help confirm how the flare and eruption are initiated and how they evolve. Here we use the Atmospheric Imaging Assembly to analyse a flare and erupting flux rope on 2014-April-18, while observations from the Nancay Radio Astronomy Facility allows us to diagnose the sites of electron acceleration during the eruption. Our analysis shows evidence for a pre-formed flux rope which slowly rises and becomes destabilised at the time of a C-class flare, plasma jet and the escape of >75 keV electrons from rope center into the corona. As the eruption proceeds, continued acceleration of electrons with energies of ~5 keV occurs above the flux rope for a period over 5 minutes. At flare peak, one site of electron acceleration is located close to the flare site while another is driven by the erupting flux rope into the corona at speeds of up to 400 km/s. Energetic electrons then fill the erupting volume, eventually allowing the flux rope legs to be clearly imaged from radio sources at 150-445MHz. Following the analysis of Joshi et al. (2015), we conclude that the sites of energetic electrons are consistent with flux rope eruption via a tether-cutting or flux cancellation scenario inside a magnetic fan-spine structure. In total, our radio observations allow us to better understand the evolution of a flux rope eruption and its associated electron acceleration sites, from eruption initiation to propagation into the corona.
[31]  [pdf] - 1490851
A Comparison of Flare Forecasting Methods, I: Results from the "All-Clear" Workshop
Comments: Accepted for publication in ApJ
Submitted: 2016-08-22
Solar flares produce radiation which can have an almost immediate effect on the near-Earth environment, making it crucial to forecast flares in order to mitigate their negative effects. The number of published approaches to flare forecasting using photospheric magnetic field observations has proliferated, with varying claims about how well each works. Because of the different analysis techniques and data sets used, it is essentially impossible to compare the results from the literature. This problem is exacerbated by the low event rates of large solar flares. The challenges of forecasting rare events have long been recognized in the meteorology community, but have yet to be fully acknowledged by the space weather community. During the interagency workshop on "all clear" forecasts held in Boulder, CO in 2009, the performance of a number of existing algorithms was compared on common data sets, specifically line-of-sight magnetic field and continuum intensity images from MDI, with consistent definitions of what constitutes an event. We demonstrate the importance of making such systematic comparisons, and of using standard verification statistics to determine what constitutes a good prediction scheme. When a comparison was made in this fashion, no one method clearly outperformed all others, which may in part be due to the strong correlations among the parameters used by different methods to characterize an active region. For M-class flares and above, the set of methods tends towards a weakly positive skill score (as measured with several distinct metrics), with no participating method proving substantially better than climatological forecasts.
[32]  [pdf] - 1460966
Quasi-Periodic Pulsations during the Impulsive and Decay phases of an X-class Flare
Comments: 7 pages, 4 figures, 1 table
Submitted: 2016-07-23
Quasi-periodic pulsations (QPP) are often observed in X-ray emission from solar flares. To date, it is unclear what their physical origins are. Here, we present a multi-instrument investigation of the nature of QPP during the impulsive and decay phases of the X1.0 flare of 28 October 2013. We focus on the character of the fine structure pulsations evident in the soft X-ray time derivatives and compare this variability with structure across multiple wavelengths including hard X-ray and microwave emission. We find that during the impulsive phase of the flare, high correlations between pulsations in the thermal and non-thermal emissions are seen. A characteristic timescale of ~20s is observed in all channels and a second timescale of ~55s is observed in the non-thermal emissions. Soft X-ray pulsations are seen to persist into the decay phase of this flare, up to 20 minutes after the non-thermal emission has ceased. We find that these decay phase thermal pulsations have very small amplitude and show an increase in characteristic timescale from ~40s up to ~70s. We interpret the bursty nature of the co-existing multi-wavelength QPP during the impulsive phase in terms of episodic particle acceleration and plasma heating. The persistent thermal decay phase QPP are most likely connected with compressive MHD processes in the post-flare loops such as the fast sausage mode or the vertical kink mode.
[33]  [pdf] - 1441089
Long-Term Tracking of Corotating Density Structures using Heliospheric Imaging
Comments: 27 pages, 12 figures, published in Solar Physics
Submitted: 2016-06-03, last modified: 2016-07-20
The systematic monitoring of the solar wind in high-cadence and high-resolution heliospheric images taken by the Solar-Terrestrial Relation Observatory (STEREO) spacecraft permits the study of the spatial and temporal evolution of variable solar wind flows from the Sun out to 1~AU, and beyond. As part of the EU Framework 7 (FP7) Heliospheric Cataloguing, Analysis and Techniques Service (HELCATS) project, we have generated a catalogue listing the properties of 190 corotating structures well-observed in images taken by the Heliospheric Imager instruments on-board STEREO-A. We present here one of very few long-term analyses of solar wind structures advected by the background solar wind. This analysis confirms that most of the corotating density structures detected by the heliospheric imagers comprises a series of density inhomogeneities advected by the slow solar wind that eventually become entrained by stream interaction regions. We have derived the spatial-temporal evolution of each of these corotating density structures by using a well-established fitting technique. The mean radial propagation speed of the corotating structures is found to be $311 \pm 31$ km~s$^{-1}$. We predicted the arrival time of each corotating density structure at different probes. We show that the speeds of the corotating density structures derived using our fitting technique track well the long-term variation of the radial speed of the slow solar wind during solar minimum years (2007--2008). Furthermore, we demonstrate that these features originate near the coronal neutral line that eventually becomes the heliospheric current sheet.
[34]  [pdf] - 1432440
Flaring Rates and the Evolution of Sunspot Group McIntosh Classifications
Comments: 31 pages, 11 figures
Submitted: 2016-07-04
Sunspot groups are the main source of solar flares, with the energy to power them being supplied by magnetic-field evolution (e.g. flux emergence or twisting/shearing). To date, few studies have investigated the statistical relation between sunspot-group evolution and flaring, with none considering evolution in the McIntosh classification scheme. Here we present a statistical analysis of sunspot groups from Solar Cycle 22, focusing on 24-hour changes in the three McIntosh classification components. Evolution-dependent >C1.0, >M1.0, and >X1.0 flaring rates are calculated, leading to the following results: (i) flaring rates become increasingly higher for greater degrees of upward evolution through the McIntosh classes, with the opposite found for downward evolution; (ii) the highest flaring rates are found for upward evolution from larger, more complex, classes (e.g. Zurich D- and E-classes evolving upward to F-class produce >C1.0 rates of 2.66 +/- 0.28 and 2.31 +/- 0.09 flares per 24 hours, respectively); (iii) increasingly complex classes give higher rates for all flare magnitudes, even when sunspot groups do not evolve over 24 hours. These results support the hypothesis that injection of magnetic energy by flux emergence (i.e. increasing in Zurich or compactness classes) leads to a higher frequency and magnitude of flaring.
[35]  [pdf] - 1391948
Conditions for electron-cyclotron maser emission in the solar corona
Comments: 4 pages, 3 figures
Submitted: 2016-04-16
Context. The Sun is an active source of radio emission ranging from long duration radio bursts associated with solar flares and coronal mass ejections to more complex, short duration radio bursts such as solar S bursts, radio spikes and fibre bursts. While plasma emission is thought to be the dominant emission mechanism for most radio bursts, the electron-cyclotron maser (ECM) mechanism may be responsible for more complex, short-duration bursts as well as fine structures associated with long-duration bursts. Aims. We investigate the conditions for ECM in the solar corona by considering the ratio of the electron plasma frequency {\omega}p to the electron-cyclotron frequency {\Omega}e. The ECM is theoretically possible when {\omega}p/{\Omega}e < 1. Methods. Two-dimensional electron density, magnetic field, plasma frequency, and electron cyclotron frequency maps of the off- limb corona were created using observations from SDO/AIA and SOHO/LASCO, together with potential field extrapolations of the magnetic field. These maps were then used to calculate {\omega}p/{\Omega}e and Alfven velocity maps of the off-limb corona. Results. We found that the condition for ECM emission ({\omega}p/{\Omega}e < 1) is possible at heights < 1.07 R_sun in an active region near the limb; that is, where magnetic field strengths are > 40 G and electron densities are greater than 3x10^8 cm-3. In addition, we found comparatively high Alfv\'en velocities (> 0.02 c or > 6000 km s-1) at heights < 1.07 R_sun within the active region. Conclusions. This demonstrates that the condition for ECM emission is satisfied within areas of the corona containing large magnetic fields, such as the core of a large active region. Therefore, ECM could be a possible emission mechanism for high-frequency radio and microwave bursts.
[36]  [pdf] - 1342999
Performance of Major Flare Watches from the Max Millennium Program (2001-2010)
Comments: 19 pages, 2 figures, accepted for publication in Solar Physics
Submitted: 2015-12-14
The physical processes that trigger solar flares are not well understood and significant debate remains around processes governing particle acceleration, energy partition, and particle and energy transport. Observations at high resolution in energy, time, and space are required in multiple energy ranges over the whole course of many flares in order to build an understanding of these processes. Obtaining high-quality, co-temporal data from ground- and space- based instruments is crucial to achieving this goal and was the primary motivation for starting the Max Millennium program and Major Flare Watch (MFW) alerts, aimed at coordinating observations of all flares >X1 GOES X-ray classification (including those partially occulted by the limb). We present a review of the performance of MFWs from 1 February 2001 to 31 May 2010, inclusive, that finds: (1) 220 MFWs were issued in 3,407 days considered (6.5% duty cycle), with these occurring in 32 uninterrupted periods that typically last 2-8 days; (2) 56% of flares >X1 were caught, occurring in 19% of MFW days; (3) MFW periods ended at suitable times, but substantial gain could have been achieved in percentage of flares caught if periods had started 24 h earlier; (4) MFWs successfully forecast X-class flares with a true skill statistic (TSS) verification metric score of 0.500, that is comparable to a categorical flare/no-flare interpretation of the NOAA Space Weather Prediction Centre probabilistic forecasts (TSS = 0.488).
[37]  [pdf] - 1342883
Automatic Detection of Magnetic delta in Sunspot Groups
Comments: 14 pages, 5 figures, accepted for publication in Solarphysics
Submitted: 2015-10-21
Large and magnetically complex sunspot groups are known to be associated with flares. To date, the Mount Wilson scheme has been used to classify sunspot groups based on their morphological and magnetic properties. The most flare prolific class, the delta sunspot-group, is characterised by opposite polarity umbrae within a common penumbra, separated by less than 2 degrees. In this article, we present a new system, called the Solar Monitor Active Region Tracker - Delta Finder (SMART-DF), that can be used to automatically detect and classify magnetic deltas in near-realtime. Using continuum images and magnetograms from the Helioseismic and Magnetic Imager (HMI) onboard NASA's Solar Dynamics Observatory (SDO), we first estimate distances between opposite polarity umbrae. Opposite polarity pairs having distances of less that 2 degrees are then identified, and if these pairs are found to share a common penumbra, they are identified as a magnetic delta configuration. The algorithm was compared to manual delta detections reported by the Space Weather Prediction Center (SWPC), operated by the National Oceanic and Atmospheric Administration (NOAA). SMART-DF detected 21 out of 23 active regions (ARs) that were marked as delta spots by NOAA during 2011 - 2012 (within +/- 60 degrees longitude). SMART-DF in addition detected five ARs which were not announced as delta spots by NOAA. The near-relatime operation of SMART-DF resulted in many deltas being identified in advance of NOAA's daily notification. SMART-DF will be integrated with SolarMonitor ( and the near-realtime information will be available to the public.
[38]  [pdf] - 1280963
Low frequency radio observations of bi-directional electron beams in the solar corona
Submitted: 2015-08-05
The radio signature of a shock travelling through the solar corona is known as a type II solar radio burst. In rare cases these bursts can exhibit a fine structure known as `herringbones', which are a direct indicator of particle acceleration occurring at the shock front. However, few studies have been performed on herringbones and the details of the underlying particle acceleration processes are unknown. Here, we use an image processing technique known as the Hough transform to statistically analyse the herringbone fine structure in a radio burst at $\sim$20-90 MHz observed from the Rosse Solar-Terrestrial Observatory on 2011 September 22. We identify 188 individual bursts which are signatures of bi-directional electron beams continuously accelerated to speeds of 0.16$_{-0.10}^{+0.11} c$. This occurs at a shock acceleration site initially at a constant altitude of $\sim$0.6 R$_{\odot}$ in the corona, followed by a shift to $\sim$0.5 R$_{\odot}$. The anti-sunward beams travel a distance of 170$_{-97}^{+174}$ Mm (and possibly further) away from the acceleration site, while those travelling toward the sun come to a stop sooner, reaching a smaller distance of 112$_{-76}^{+84}$ Mm. We show that the stopping distance for the sunward beams may depend on the total number density and the velocity of the beam. Our study concludes that a detailed statistical analysis of herringbone fine structure can provide information on the physical properties of the corona which lead to these relatively rare radio bursts.
[39]  [pdf] - 1255327
LOFAR tied-array imaging and spectroscopy of solar S bursts
Comments: 6 pages, 4 figures
Submitted: 2015-07-27
Context. The Sun is an active source of radio emission that is often associated with energetic phenomena ranging from nanoflares to coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), numerous millisecond duration radio bursts have been reported, such as radio spikes or solar S bursts (where S stands for short). To date, these have neither been studied extensively nor imaged because of the instrumental limitations of previous radio telescopes. Aims. Here, Low Frequency Array (LOFAR) observations were used to study the spectral and spatial characteristics of a multitude of S bursts, as well as their origin and possible emission mechanisms. Methods. We used 170 simultaneous tied-array beams for spectroscopy and imaging of S bursts. Since S bursts have short timescales and fine frequency structures, high cadence (~50 ms) tied-array images were used instead of standard interferometric imaging, that is currently limited to one image per second. Results. On 9 July 2013, over 3000 S bursts were observed over a time period of ~8 hours. S bursts were found to appear as groups of short-lived (<1 s) and narrow-bandwidth (~2.5 MHz) features, the majority drifting at ~3.5 MHz/s and a wide range of circular polarisation degrees (2-8 times more polarised than the accompanying Type III bursts). Extrapolation of the photospheric magnetic field using the potential field source surface (PFSS) model suggests that S bursts are associated with a trans-equatorial loop system that connects an active region in the southern hemisphere to a bipolar region of plage in the northern hemisphere. Conclusions. We have identified polarised, short-lived solar radio bursts that have never been imaged before. They are observed at a height and frequency range where plasma emission is the dominant emission mechanism, however they possess some of the characteristics of electron-cyclotron maser emission.
[40]  [pdf] - 1240974
Solar and Heliospheric Physics with the Square Kilometre Array
Comments: 15 pages, Proceedings of Advancing Astrophysics with the Square Kilometre Array (AASKA14). 9 -13 June, 2014. Giardini Naxos, Italy. Online at, id.169
Submitted: 2015-07-02
The fields of solar radiophysics and solar system radio physics, or radio heliophysics, will benefit immensely from an instrument with the capabilities projected for SKA. Potential applications include interplanetary scintillation (IPS), radio-burst tracking, and solar spectral radio imaging with a superior sensitivity. These will provide breakthrough new insights and results in topics of fundamental importance, such as the physics of impulsive energy releases, magnetohydrodynamic oscillations and turbulence, the dynamics of post-eruptive processes, energetic particle acceleration, the structure of the solar wind and the development and evolution of solar wind transients at distances up to and beyond the orbit of the Earth. The combination of the high spectral, time and spatial resolution and the unprecedented sensitivity of the SKA will radically advance our understanding of basic physical processes operating in solar and heliospheric plasmas and provide a solid foundation for the forecasting of space weather events.
[41]  [pdf] - 1237823
Ellerman Bombs with Jets: Cause and Effect
Comments: 9 pages, 9 figures, 2 supplementary movies. Accepted to ApJ
Submitted: 2015-03-18
Ellerman Bombs (EBs) are thought to arise as a result of photospheric magnetic reconnection. We use data from the Swedish 1-m Solar Telescope (SST), to study EB events on the solar disk and at the limb. Both datasets show that EBs are connected to the foot-points of forming chromospheric jets. The limb observations show that a bright structure in the H$\alpha$ blue wing connects to the EB initially fuelling it, leading to the ejection of material upwards. The material moves along a loop structure where a newly formed jet is subsequently observed in the red wing of H$\alpha$. In the disk dataset, an EB initiates a jet which propagates away from the apparent reconnection site within the EB flame. The EB then splits into two, with associated brightenings in the inter-granular lanes (IGLs). Micro-jets are then observed, extending to 500 km with a lifetime of a few minutes. Observed velocities of the micro-jets are approximately 5-10 km s$^{-1}$, while their chromospheric counterparts range from 50-80 km s$^{-1}$. MURaM simulations of quiet Sun reconnection show that micro-jets with similar properties to that of the observations follow the line of reconnection in the photosphere, with associated H$\alpha$ brightening at the location of increased temperature.
[42]  [pdf] - 933054
Stable Umbral Chromospheric Structures
Comments: Accepted to Astronomy and Astrophysics. Online material ( and will be available at A&A
Submitted: 2014-12-18
Aims. To understand the morphology of the chromosphere in sunspot umbra. We investigate if the horizontal structures observed in the spectral core of the Ca II H line are ephemeral visuals caused by the shock dynamics of more stable structures, and examine their relationship with observables in the H-alpha line. Methods. Filtergrams in the core of the Ca II H and H-alpha lines as observed with the Swedish 1-m Solar Telescope are employed. We utilise a technique that creates composite images and tracks the flash propagation horizontally. Results. We find 0"15 wide horizontal structures, in all of the three target sunspots, for every flash where the seeing was moderate to good. Discrete dark structures are identified that are stable for at least two umbral flashes, as well as systems of structures that live for up to 24 minutes. We find cases of extremely extended structures with similar stability, with one such structure showing an extent of 5". Some of these structures have a correspondence in H-alpha but we were unable to find a one to one correspondence for every occurrence. If the dark streaks are formed at the same heights as umbral flashes then there are systems of structures with strong departures from the vertical for all three analysed sunspots. Conclusions. Long-lived Ca II H filamentary horizontal structures are a common and likely ever-present feature in the umbra of sunspots. If the magnetic field in the chromosphere of the umbra is indeed aligned with the structures, then the present theoretical understanding of the typical umbra needs to be revisited.
[43]  [pdf] - 1223564
The energetics of a global shock wave in the low solar corona
Comments: 12 pages, 6 figures. Accepted for publication in The Astrophysical Journal
Submitted: 2014-12-09
As the most energetic eruptions in the solar system, coronal mass ejections (CMEs) can produce shock waves at both their front and flanks as they erupt from the Sun into the heliosphere. However, the amount of energy produced in these eruptions, and the proportion of their energy required to produce the waves, is not well characterised. Here we use observations of a solar eruption from 2014 February 25 to estimate the energy budget of an erupting CME and the globally-propagating "EIT wave" produced by the rapid expansion of the CME flanks in the low solar corona. The "EIT wave" is shown using a combination of radio spectra and extreme ultraviolet images to be a shock front with a Mach number greater than one. Its initial energy is then calculated using the Sedov-Taylor blast-wave approximation, which provides an approximation for a shock front propagating through a region of variable density. This approach provides an initial energy estimate of $\approx$2.8$\times$10$^{31}$ ergs to produce the "EIT wave", which is approximately 10% the kinetic energy of the associated CME (shown to be $\approx$2.5$\times$10$^{32}$ ergs). These results indicate that the energy of the "EIT wave" may be significant and must be considered when estimating the total energy budget of solar eruptions.
[44]  [pdf] - 904081
Solar Hard X-ray Source Sizes in a Beam-Heated and Ionised Chromosphere
Submitted: 2014-11-19
Solar flare hard X-rays (HXRs) are produced as bremsstrahlung when an accelerated population of electrons interacts with the dense chromospheric plasma. HXR observations presented by using the Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) have shown that HXR source sizes are 3-6 times more extended in height than those predicted by the standard collisional thick target model (CTTM). Several possible explanations have been put forward including the multi-threaded nature of flare loops, pitch-angle scattering, and magnetic mirroring. However, the nonuniform ionisation (NUI) structure along the path of the electron beam has not been fully explored as a solution to this problem. Ionised plasma is known to be less effective at producing nonthermal bremsstrahlung HXRs when compared to neutral plasma. If the peak HXR emission was produced in a locally ionised region within the chromosphere, the intensity of emission will be preferentially reduced around this peak, resulting in a more extended source. Due to this effect, along with the associated density enhancement in the upper chromosphere, injection of a beam of electrons into a partially ionised plasma should result in a HXR source which is substantially more vertically extended relative to that for a neutral target. Here we present the results of a modification to the CTTM which takes into account both a localised form of chromospheric NUI and an increased target density. We find 50 keV HXR source widths, with and without the inclusion of a locally ionised region, of ~3 Mm and ~0.7 Mm, respectively. This helps to provide a theoretical solution to the currently open question of overly-extended HXR sources.
[45]  [pdf] - 1216916
Understanding CME and associated shock in the solar corona by merging multi wavelengths observation
Comments: 17 pages, 13 figures
Submitted: 2014-09-12
Using multi-wavelength imaging observations, in EUV, white light and radio, and radio spectral data over a large frequency range, we analyzed the triggering and development of a complex eruptive event. This one includes two components, an eruptive jet and a CME which interact during more than 30 min, and can be considered as physically linked. This was an unusual event. The jet is generated above a typical complex magnetic configuration which has been investigated in many former studies related to the build-up of eruptive jets; this configuration includes fan-field lines originating from a corona null point above a parasitic polarity, which is embedded in one polarity region of large Active Region (AR). The initiation and development of the CME, observed first in EUV, does not show usual signatures. In this case, the eruptive jet is the main actor of this event. The CME appears first as a simple loop system which becomes destabilized by magnetic reconnection between the outer part of the jet and the ambient medium. The progression of the CME is closely associated with the occurrence of two successive types II bursts from distinct origin. An important part of this study is the first radio type II burst for which the joint spectral and imaging observations allowed: i) to follow, step by step, the evolution of the spectrum and of the trajectory of the radio burst, in relationship with the CME evolution; ii) to obtain, without introducing an electronic density model, the B-field and the Alfven speed.
[46]  [pdf] - 863117
LOFAR tied-array imaging of Type III solar radio bursts
Comments: 9 pages, 6 figures
Submitted: 2014-07-16
The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.
[47]  [pdf] - 1209911
Quasiperiodic acceleration of electrons by a plasmoid-driven shock in the solar atmosphere
Submitted: 2014-06-03
Cosmic rays and solar energetic particles may be accelerated to relativistic energies by shock waves in astrophysical plasmas. On the Sun, shocks and particle acceleration are often associated with the eruption of magnetized plasmoids, called coronal mass ejections (CMEs). However, the physical relationship between CMEs and shock particle acceleration is not well understood. Here, we use extreme ultraviolet, radio and white-light imaging of a solar eruptive event on 22 September 2011 to show that a CME-induced shock (Alfv\'en Mach number 2.4$^{+0.7}_{-0.8}$) was coincident with a coronal wave and an intense metric radio burst generated by intermittent acceleration of electrons to kinetic energies of 2-46 keV (0.1-0.4 c). Our observations show that plasmoid-driven quasi-perpendicular shocks are capable of producing quasi-periodic acceleration of electrons, an effect consistent with a turbulent or rippled plasma shock surface.
[48]  [pdf] - 1208606
CorPITA: An Automated Algorithm for the Identification and Analysis of Coronal "EIT Waves"
Comments: 20 pages, 8 figures, accepted for publication in Solar Physics
Submitted: 2014-03-26
The continuous stream of data available from the Atmospheric Imaging Assembly (AIA) telescopes onboard the Solar Dynamics Observatory (SDO) spacecraft has allowed a deeper understanding of the Sun. However, the sheer volume of data has necessitated the development of automated techniques to identify and analyse various phenomena. In this article, we describe the Coronal Pulse Identification and Tracking Algorithm (CorPITA) for the identification and analysis of coronal "EIT waves". CorPITA uses an intensity-profile technique to identify the propagating pulse, tracking it throughout its evolution before returning estimates of its kinematics. The algorithm is applied here to a data-set from February 2011, allowing its capabilities to be examined and critiqued. This algorithm forms part of the SDO Feature Finding Team initiative and will be implemented as part of the Heliophysics Event Knowledgebase (HEK). This is the first fully automated algorithm to identify and track the propagating "EIT wave" rather than any associated phenomena and will allow a deeper understanding of this controversial phenomenon.
[49]  [pdf] - 820784
The formation heights of coronal shocks from 2D density and Alfv\'en speed maps
Submitted: 2014-02-17, last modified: 2014-03-10
Super-Alfv\'enic shock waves associated with coronal mass ejections (CMEs) can produce radio emission known as Type II bursts. In the absence of direct imaging, accurate estimates of coronal electron densities, magnetic field strengths and Alfv\'en speeds are required in order to calculate the kinematics of shocks. To date, 1D radial models have been used, but these are not appropriate for shocks propagating in non-radial directions. Here, we study a coronal shock wave associated with a CME and Type II radio burst using 2D electron density and Alfv\'en speed maps to determine the locations that shocks are excited as the CME expands through the corona. Coronal density maps were obtained from emission measures derived from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) and polarized brightness measurements from the Large Angle and Spectrometric Coronagraph (LASCO) on board the Solar and Heliospheric Observatory (SOHO). Alfv\'en speed maps were calculated using these density maps and magnetic field extrapolations from the Helioseismic and Magnetic Imager (SDO/HMI). The computed density and Alfv\'en speed maps were then used to calculate the shock kinematics in non-radial directions. Using the kinematics of the Type II burst and associated shock, we find our observations to be consistent with the formation of a shock located at the CME flanks where the Alfv\'en speed has a local minimum. 1D density models are not appropriate for shocks that propagate non-radially along the flanks of a CME. Rather, the 2D density, magnetic field and Alfv\'en speed maps described here give a more accurate method for determining the fundamental properties of shocks and their relation to CMEs.
[50]  [pdf] - 1202799
The Compatibility of Flare Temperatures Observed with AIA, GOES, and RHESSI
Comments: Accepted by Solar Physics Journal
Submitted: 2014-01-16
We test the compatibility and biases of multi-thermal flare DEM (differential emission measure) peak temperatures determined with AIA with those determined by GOES and RHESSI using the isothermal assumption. In a set of 149 M- and X-class flares observed during the first two years of the SDO mission, AIA finds DEM peak temperatures at the time of the peak GOES 1-8 A flux to have an average of Tp = 12.0+/-2.9 MK and Gaussian DEM widths of log10(sigma_T) = 0.50+/-0.13. From GOES observations of the same 149 events, a mean temperature of Tp = 15.6+/-2.4 MK is inferred, which is systematically higher by a factor of TGOES/TAIA = 1.4+/-0.4. We demonstrate that this discrepancy results from the isothermal assumption in the inversion of the GOES filter ratio. From isothermal fits to photon spectra at energies of E ~ 6-12 keV of 61 of these events, RHESSI finds the temperature to be higher still by a factor of TRHESSI/TAIA = 1.9+/-1.0. We find that this is partly a consequence of the isothermal assumption. However, RHESSI is not sensitive to the low-temperature range of the DEM peak, and thus RHESSI samples only the high-temperature tail of the DEM function. This can also contribute to the discrepancy between AIA and RHESSI temperatures. The higher flare temperatures found by GOES and RHESSI imply correspondingly lower emission measures. We conclude that self-consistent flare DEM temperatures and emission measures require simultaneous fitting of EUV (AIA) and soft X-ray (GOES and RHESSI) fluxes.
[51]  [pdf] - 1173090
Improved methods for determining the kinematics of coronal mass ejections and coronal waves
Comments: 12 pages, 12 figures
Submitted: 2013-07-30
The study of solar eruptive events and associated phenomena is of great importance in the context of solar and heliophysics. Coronal mass ejections (CMEs) and coronal waves are energetic manifestations of the restructuring of the solar magnetic field and mass motion of the plasma. Characterising this motion is vital for deriving the dynamics of these events and thus understanding the physics driving their initiation and propagation. The development and use of appropriate methods for measuring event kinematics is therefore imperative. Traditional approaches to the study of CME and coronal wave kinematics do not return wholly accurate nor robust estimates of the true event kinematics and associated uncertainties. We highlight the drawbacks of these approaches, and demonstrate improved methods for accurate and reliable determination of the kinematics. The Savitzky-Golay filter is demonstrated as a more appropriate fitting technique for CME and coronal wave studies, and a residual resampling bootstrap technique is demonstrated as a statistically rigorous method for the determination of kinematic error estimates and goodness-of-fit tests. It is shown that the scatter on distance-time measurements of small sample size can significantly limit the ability to derive accurate and reliable kinematics. This may be overcome by (i) increasing measurement precision and sampling cadence, and (ii) applying robust methods for deriving the kinematics and reliably determining their associated uncertainties. If a priori knowledge exists and a pre-determined model form for the kinematics is available (or indeed any justified fitting-form to be tested against the data), then its precision can be examined using a bootstrapping technique to determine the confidence interval associated with the model/fitting parameters.
[52]  [pdf] - 1166462
Solar Flare X-ray Source Motion as a Response to Electron Spectral Hardening
Comments: 8 pages, 5 figures
Submitted: 2013-05-07
Context: Solar flare hard X-rays (HXRs) are thought to be produced by nonthermal coronal electrons stopping in the chromosphere, or remaining trapped in the corona. The collisional thick target model (CTTM) predicts that sources produced by harder power-law injection spectra should appear further down the legs or footpoints of a flare loop. Therefore, hardening of the injected power-law electron spectrum during flare onset should be concurrent with a descending hard X-ray source. Aims: To test this implication of the CTTM by comparing its predicted HXR source locations with those derived from observations of a solar flare which exhibits a nonthermally-dominated spectrum before the peak in HXRs, known as an early impulsive event. Methods: HXR images and spectra of an early impulsive C-class flare were obtained using the Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). Images were reconstructed to produce HXR source height evolutions for three energy bands. Spatially-integrated spectral analysis was performed to isolate nonthermal emission, and to determine the power-law index of the electron injection spectrum. The observed height-time evolutions were then fit with CTTM-based simulated heights for each energy. Results: A good match between model and observed source heights was reached, requiring a density model that agreed well with previous studies of flare loop densities. Conclusions: The CTTM has been used to produce a descent of model HXR source heights that compares well with observations of this event. Based on this interpretation, downward motion of nonthermal sources should indeed occur in any flare where there is spectral hardening in the electron distribution during a flare. However, this would often be masked by thermal emission associated with flare plasma pre-heating.
[53]  [pdf] - 635488
The Thermal Properties of Solar Flares Over Three Solar Cycles Using GOES X-ray Observations
Submitted: 2012-06-05, last modified: 2013-03-06
Solar flare X-ray emission results from rapidly increasing temperatures and emission measures in flaring active region loops. To date, observations from the X-Ray Sensor (XRS) onboard the Geostationary Operational Environmental Satellite (GOES) have been used to derive these properties, but have been limited by a number of factors, including the lack of a consistent background subtraction method capable of being automatically applied to large numbers of flares. In this paper, we describe an automated temperature and emission measure-based background subtraction method (TEBBS), which builds on the methods of Bornmann (1990). Our algorithm ensures that the derived temperature is always greater than the instrumental limit and the pre-flare background temperature, and that the temperature and emission measure are increasing during the flare rise phase. Additionally, TEBBS utilizes the improved estimates of GOES temperatures and emission measures from White et al. (2005). TEBBS was successfully applied to over 50,000 solar flares occurring over nearly three solar cycles (1980-2007), and used to create an extensive catalog of the solar flare thermal properties. We confirm that the peak emission measure and total radiative losses scale with background subtracted GOES X-ray flux as power-laws, while the peak temperature scales logarithmically. As expected, the peak emission measure shows an increasing trend with peak temperature, although the total radiative losses do not. While these results are comparable to previous studies, we find that flares of a given GOES class have lower peak temperatures and higher peak emission measures than previously reported. The resulting TEBBS database of thermal flare plasma properties is publicly available on Solar Monitor ( and will be available on Heliophysics Integrated Observatory (
[54]  [pdf] - 945382
Evidence for Partial Taylor Relaxation from Changes in Magnetic Geometry and Energy during a Solar Flare
Comments: Accepted for Publication in Astronomy & Astrophysics. 11 pages, 11 figures
Submitted: 2012-12-24
Solar flares are powered by energy stored in the coronal magnetic field, a portion of which is released when the field reconfigures into a lower energy state. Investigation of sunspot magnetic field topology during flare activity is useful to improve our understanding of flaring processes. Here we investigate the deviation of the non-linear field configuration from that of the linear and potential configurations, and study the free energy available leading up to and after a flare. The evolution of the magnetic field in NOAA region 10953 was examined using data from Hinode/SOT-SP, over a period of 12 hours leading up to and after a GOES B1.0 flare. Previous work on this region found pre- and post-flare changes in photospheric vector magnetic field parameters of flux elements outside the primary sunspot. 3D geometry was thus investigated using potential, linear force-free, and non-linear force-free field extrapolations in order to fully understand the evolution of the field lines. Traced field line geometrical and footpoint orientation differences show that the field does not completely relax to a fully potential or linear force-free state after the flare. Magnetic and free magnetic energies increase significantly ~ 6.5-2.5 hours before the flare by ~ 10^31 erg. After the flare, the non-linear force-free magnetic energy and free magnetic energies decrease but do not return to pre-flare 'quiet' values. The post-flare non-linear force-free field configuration is closer (but not equal) to that of the linear force-free field configuration than a potential one. However, the small degree of similarity suggests that partial Taylor relaxation has occurred over a time scale of ~ 3-4 hours.
[55]  [pdf] - 1150842
The SWAP EUV Imaging Telescope Part I: Instrument Overview and Pre-Flight Testing
Comments: 26 pages, 9 figures, 1 movie
Submitted: 2012-08-22
The Sun Watcher with Active Pixels and Image Processing (SWAP) is an EUV solar telescope on board ESA's Project for Onboard Autonomy 2 (PROBA2) mission launched on 2 November 2009. SWAP has a spectral bandpass centered on 17.4 nm and provides images of the low solar corona over a 54x54 arcmin field-of-view with 3.2 arcsec pixels and an imaging cadence of about two minutes. SWAP is designed to monitor all space-weather-relevant events and features in the low solar corona. Given the limited resources of the PROBA2 microsatellite, the SWAP telescope is designed with various innovative technologies, including an off-axis optical design and a CMOS-APS detector. This article provides reference documentation for users of the SWAP image data.
[56]  [pdf] - 1117765
Observations of Low Frequency Solar Radio Bursts from the Rosse Solar-Terrestrial Observatory
Submitted: 2012-04-04, last modified: 2012-07-26
The Rosse Solar-Terrestrial Observatory (RSTO; was established at Birr Castle, Co. Offaly, Ireland (53 05'38.9", 7 55'12.7") in 2010 to study solar radio bursts and the response of the Earth's ionosphere and geomagnetic field. To date, three Compound Astronomical Low-cost Low-frequency Instrument for Spectroscopy and Transportable Observatory (CALLISTO) spectrometers have been installed, with the capability of observing in the frequency range 10-870 MHz. The receivers are fed simultaneously by biconical and log-periodic antennas. Nominally, frequency spectra in the range 10-400 MHz are obtained with 4 sweeps per second over 600 channels. Here, we describe the RSTO solar radio spectrometer set-up, and present dynamic spectra of a sample of Type II, III and IV radio bursts. In particular, we describe fine-scale structure observed in Type II bursts, including band splitting and rapidly varying herringbone features.
[57]  [pdf] - 542292
Automatic Detection and Tracking of CMEs II: Multiscale Filtering of Coronagraph Data
Comments: 18 pages, 10 figures
Submitted: 2012-07-25
Studying CMEs in coronagraph data can be challenging due to their diffuse structure and transient nature, and user-specific biases may be introduced through visual inspection of the images. The large amount of data available from the SOHO, STEREO, and future coronagraph missions, also makes manual cataloguing of CMEs tedious, and so a robust method of detection and analysis is required. This has led to the development of automated CME detection and cata- loguing packages such as CACTus, SEEDS and ARTEMIS. Here we present the development of a new CORIMP (coronal image processing) CME detection and tracking technique that overcomes many of the drawbacks of current catalogues. It works by first employing the dynamic CME separation technique outlined in a companion paper, and then characterising CME structure via a multiscale edge-detection algorithm. The detections are chained through time to determine the CME kinematics and morphological changes as it propagates across the plane-of-sky. The effectiveness of the method is demonstrated by its application to a selection of SOHO/LASCO and STEREO/SECCHI images, as well as to synthetic coronagraph images created from a model corona with a variety of CMEs. The algorithms described in this article are being applied to the whole LASCO and SECCHI datasets, and a catalogue of results will soon be available to the public.
[58]  [pdf] - 528814
A Bayesian approach to comparing theoretic models to observational data: A case study from solar flare physics
Submitted: 2011-02-01, last modified: 2012-06-24
Solar flares are large-scale releases of energy in the solar atmosphere, which are characterised by rapid changes in the hydrodynamic properties of plasma from the photosphere to the corona. Solar physicists have typically attempted to understand these complex events using a combination of theoretical models and observational data. From a statistical perspective, there are many challenges associated with making accurate and statistically significant comparisons between theory and observations, due primarily to the large number of free parameters associated with physical models. This class of ill-posed statistical problem is ideally suited to Bayesian methods. In this paper, the solar flare studied by Raftery et al. (2009) is reanalysed using a Bayesian framework. This enables us to study the evolution of the flare's temperature, emission measure and energy loss in a statistically self-consistent manner. The Bayesian-based model selection techniques imply that no decision can be made regarding which of the conductive or non-thermal beam heating play the most important role in heating the flare plasma during the impulsive phase of this event.
[59]  [pdf] - 1118151
Coronal Mass Ejection Mass, Energy, and Force Estimates Using STEREO
Submitted: 2012-04-20
Understanding coronal mass ejection (CME) energetics and dynamics has been a long-standing problem, and although previous observational estimates have been made, such studies have been hindered by large uncertainties in CME mass. Here, the two vantage points of the Solar Terrestrial Relations Observatory (STEREO) COR1 and COR2 coronagraphs were used to accurately estimate the mass of the 2008 December 12 CME. Acceleration estimates derived from the position of the CME front in 3-D were combined with the mass estimates to calculate the magnitude of the kinetic energy and driving force at different stages of the CME evolution. The CME asymptotically approaches a mass of 3.4\pm1.0x10^15 g beyond ~10 R_sun. The kinetic energy shows an initial rise towards 6.3\pm3.7x10^29 erg at ~3 R_sun, beyond which it rises steadily to 4.2\pm2.5x10^30 erg at ~18 R_sun. The dynamics are described by an early phase of strong acceleration, dominated by a force of peak magnitude of 3.4\pm2.2x10^14N at ~3 R_sun, after which a force of 3.8\pm5.4x10^13 N takes affect between ~7-18 R_sun. These results are consistent with magnetic (Lorentz) forces acting at heliocentric distances of <7 R_sun, while solar wind drag forces dominate at larger distances (>7 R_sun).
[60]  [pdf] - 481302
Toward Reliable Benchmarking of Solar Flare Forecasting Methods
Comments: 7 pages, 1 figure
Submitted: 2012-02-27
Solar flares occur in complex sunspot groups, but it remains unclear how the probability of producing a flare of a given magnitude relates to the characteristics of the sunspot group. Here, we use Geostationary Operational Environmental Satellite X-ray flares and McIntosh group classifications from solar cycles 21 and 22 to calculate average flare rates for each McIntosh class and use these to determine Poisson probabilities for different flare magnitudes. Forecast verification measures are studied to find optimum thresholds to convert Poisson flare probabilities into yes/no predictions of cycle 23 flares. A case is presented to adopt the true skill statistic (TSS) as a standard for forecast comparison over the commonly used Heidke skill score (HSS). In predicting flares over 24 hr, the maximum values of TSS achieved are 0.44 (C-class), 0.53 (M-class), 0.74 (X-class), 0.54 (>=M1.0), and 0.46 (>=C1.0). The maximum values of HSS are 0.38 (C-class), 0.27 (M-class), 0.14 (X-class), 0.28 (>=M1.0), and 0.41 (>=C1.0). These show that Poisson probabilities perform comparably to some more complex prediction systems, but the overall inaccuracy highlights the problem with using average values to represent flaring rate distributions.
[61]  [pdf] - 420116
Automated Solar Feature Detection for Space Weather Applications
Comments: 20 pages, 6 figures, Book chapter in "Applied Signal and Image Processing: Multidisciplinary Advancements" (ed. Rami Qahwaji, Roger Green and Evor L. Hines) pp. 207-225
Submitted: 2011-09-30, last modified: 2011-10-03
The solar surface and atmosphere are highly dynamic plasma environments, which evolve over a wide range of temporal and spatial scales. Large-scale eruptions, such as coronal mass ejections, can be accelerated to millions of kilometres per hour in a matter of minutes, making their automated detection and characterisation challenging. Additionally, there are numerous faint solar features, such as coronal holes and coronal dimmings, which are important for space weather monitoring and forecasting, but their low intensity and sometimes transient nature makes them problematic to detect using traditional image processing techniques. These difficulties are compounded by advances in ground- and space- based instrumentation, which have increased the volume of data that solar physicists are confronted with on a minute-by-minute basis; NASA's Solar Dynamics Observatory for example is returning many thousands of images per hour (~1.5 TB/day). This chapter reviews recent advances in the application of images processing techniques to the automated detection of active regions, coronal holes, filaments, CMEs, and coronal dimmings for the purposes of space weather monitoring and prediction.
[62]  [pdf] - 753371
An Observational Overview of Solar Flares
Comments: This is an article for a monograph on the physics of solar flares, inspired by RHESSI observations. The individual articles are to appear in Space Science Reviews (2011)
Submitted: 2011-09-27, last modified: 2011-09-28
We present an overview of solar flares and associated phenomena, drawing upon a wide range of observational data primarily from the RHESSI era. Following an introductory discussion and overview of the status of observational capabilities, the article is split into topical sections which deal with different areas of flare phenomena (footpoints and ribbons, coronal sources, relationship to coronal mass ejections) and their interconnections. We also discuss flare soft X-ray spectroscopy and the energetics of the process. The emphasis is to describe the observations from multiple points of view, while bearing in mind the models that link them to each other and to theory. The present theoretical and observational understanding of solar flares is far from complete, so we conclude with a brief discussion of models, and a list of missing but important observations.
[63]  [pdf] - 1084428
The Wave Properties of Coronal Bright Fronts Observed Using SDO/AIA
Comments: 14 pages, 2 figures. Accepted for publication in The Astrophysical Journal Letters
Submitted: 2011-09-27
Coronal bright fronts (CBFs) are large scale wavefronts that propagate though the solar corona at hundreds of kilometers per second. While their kinematics have been studied in detail, many questions remain regarding the temporal evolution of their amplitude and pulse width. Here, contemporaneous high cadence, multi-thermal observations of the solar corona from the Solar Dynamic Observatory (SDO) and Solar TErrestrial RElations Observatory (STEREO) spacecraft are used to determine the kinematics and expansion rate of a CBF wavefront observed on 2010 August 14. The CBF was found to have a lower initial velocity with weaker deceleration in STEREO observations compared to SDO (~340 km/s and -72 m/s/s as opposed to ~410 km/s and -279 m/s/s). The CBF kinematics from SDO were found to be highly passband-dependent, with an initial velocity ranging from 379+/-12 km/s to 460+/-28 km/s and acceleration ranging from -128+/-28 m/s/s to -431+/-86 m/s/s in the 335A and 304A passbands respectively. These kinematics were used to estimate a quiet coronal magnetic field strength range of ~1-2 G. Significant pulse broadening was also observed, with expansion rates of ~130 km/s (STEREO) and ~220 km/s (SDO). By treating the CBF as a linear superposition of sinusoidal waves within a Gaussian envelope, the resulting dispersion rate of the pulse was found to be ~8-13 Mm^2 s^-1. These results are indicative of a fast-mode magnetoacoustic wave pulse propagating through an inhomogeneous medium.
[64]  [pdf] - 1084244
LEMUR: Large European Module for solar Ultraviolet Research. European contribution to JAXA's Solar-C mission
Comments: 35 pages, 14 figures. To appear on Experimental Astronomy
Submitted: 2011-09-20, last modified: 2011-09-21
Understanding the solar outer atmosphere requires concerted, simultaneous solar observations from the visible to the vacuum ultraviolet (VUV) and soft X-rays, at high spatial resolution (between 0.1" and 0.3"), at high temporal resolution (on the order of 10 s, i.e., the time scale of chromospheric dynamics), with a wide temperature coverage (0.01 MK to 20 MK, from the chromosphere to the flaring corona), and the capability of measuring magnetic fields through spectropolarimetry at visible and near-infrared wavelengths. Simultaneous spectroscopic measurements sampling the entire temperature range are particularly important. These requirements are fulfilled by the Japanese Solar-C mission (Plan B), composed of a spacecraft in a geosynchronous orbit with a payload providing a significant improvement of imaging and spectropolarimetric capabilities in the UV, visible, and near-infrared with respect to what is available today and foreseen in the near future. The Large European Module for solar Ultraviolet Research (LEMUR), described in this paper, is a large VUV telescope feeding a scientific payload of high-resolution imaging spectrographs and cameras. LEMUR consists of two major components: a VUV solar telescope with a 30 cm diameter mirror and a focal length of 3.6 m, and a focal-plane package composed of VUV spectrometers covering six carefully chosen wavelength ranges between 17 and 127 nm. The LEMUR slit covers 280" on the Sun with 0.14" per pixel sampling. In addition, LEMUR is capable of measuring mass flows velocities (line shifts) down to 2 km/s or better. LEMUR has been proposed to ESA as the European contribution to the Solar C mission.
[65]  [pdf] - 1077154
STEREO Direct Imaging of a CME-driven Shock to 0.5AU
Comments: 13 pages, 5 figures (black and white), AAS LATEX v5.2
Submitted: 2011-06-08
Fast coronal mass ejections (CMEs) generate standing or bow shocks as they propagate through the corona and solar wind. Although CME shocks have previously been detected indirectly via their emission at radio frequencies, direct imaging has remained elusive due to their low contrast at optical wavelengths. Here we report the first images of a CME-driven shock as it propagates through interplanetary space from 8R_Sun to 120R_Sun (0.5AU), using observations from the STEREO Heliospheric Imager (HI). The CME was measured to have a velocity of ~1000 km/s and a Mach number of 4.1\pm1.2, while the shock front stand-off distance (Delta) was found to increase linearly to ~20 R_Sun at 0.5 AU. The normalised standoff distance (Delta/DO) showed reasonable agreement with semi-empirical relations, where D_O is the CME radius. However, when normalised using the radius of curvature, Delta/R_O did not agree well with theory, implying that R_O was under-estimated by a factor of ~3-8. This is most likely due to the difficulty in estimating the larger radius of curvature along the CME axis from the observations, which provide only a cross-sectional view of the CME.
[66]  [pdf] - 470214
The Evolution of Sunspot Magnetic Fields Associated with a Solar Flare
Comments: Accepted for Publication in Solar Physics. 16 pages, 4 figures
Submitted: 2011-05-10
Solar flares occur due to the sudden release of energy stored in active-region magnetic fields. To date, the pre-cursors to flaring are still not fully understood, although there is evidence that flaring is related to changes in the topology or complexity of an active region's magnetic field. Here, the evolution of the magnetic field in active region NOAA 10953 was examined using Hinode/SOT-SP data, over a period of 12 hours leading up to and after a GOES B1.0 flare. A number of magnetic-field properties and low-order aspects of magnetic-field topology were extracted from two flux regions that exhibited increased Ca II H emission during the flare. Pre-flare increases in vertical field strength, vertical current density, and inclination angle of ~ 8degrees towards the vertical were observed in flux elements surrounding the primary sunspot. The vertical field strength and current density subsequently decreased in the post-flare state, with the inclination becoming more horizontal by ~7degrees. This behaviour of the field vector may provide a physical basis for future flare forecasting efforts.
[67]  [pdf] - 1053515
Deceleration and Dispersion of Large-scale Coronal Bright Fronts
Comments: 15 pages, 18 figures
Submitted: 2011-04-21
One of the most dramatic manifestations of solar activity are large-scale coronal bright fronts (CBFs) observed in extreme ultraviolet (EUV) images of the solar atmosphere. To date, the energetics and kinematics of CBFs remain poorly understood, due to the low image cadence and sensitivity of previous EUV imagers and the limited methods used to extract the features. In this paper, the trajectory and morphology of CBFs was determined in order to investigate the varying properties of a sample of CBFs, including their kinematics and pulse shape, dispersion, and dissipation. We have developed a semi-automatic intensity profiling technique to extract the morphology and accurate positions of CBFs in 2.5-10 min cadence images from STEREO/EUVI. The technique was applied to sequences of 171A and 195A images from STEREO/EUVI in order to measure the wave properties of four separate CBF events. Following launch at velocities of ~240-450kms^{-1} each of the four events studied showed significant negative acceleration ranging from ~ -290 to -60ms^{-2}. The CBF spatial and temporal widths were found to increase from ~50 Mm to ~200 Mm and ~100 s to ~1500 s respectively, suggesting that they are dispersive in nature. The variation in position-angle averaged pulse-integrated intensity with propagation shows no clear trend across the four events studied. These results are most consistent with CBFs being dispersive magnetoacoustic waves.
[68]  [pdf] - 1052711
Short-term evolution of coronal hole boundaries
Submitted: 2011-03-14
The interaction of open and closed field lines at coronal hole boundaries is widely accepted to be due to interchange magnetic reconnection. To date, it is unclear how the boundaries vary on short timescales and at what velocity this occurs. Here, we describe an automated boundary tracking method used to determine coronal hole boundary displacements on short timescales. The bound- ary displacements were found to be isotropic and to have typical expansion/contraction speeds of \leq2 km s^-1, which indicate magnetic reconnection rates of \leq 3 \times 10^-3. The observed displacements were used in conjunction with the interchange reconnection model to derive typical diffusion coeffi- cients of \leq 3 \times 10^13 cm^2 s^-1. These results are consistent with an interchange reconnection process in the low corona driven by the random granular motion of open and closed fields in the photosphere.
[69]  [pdf] - 348128
Coronal Mass Ejection Detection using Wavelets, Curvelets and Ridgelets: Applications for Space Weather Monitoring
Comments: Accepted for publication in Advances in Space Research (3 April 2010)
Submitted: 2010-12-08
Coronal mass ejections (CMEs) are large-scale eruptions of plasma and magnetic feld that can produce adverse space weather at Earth and other locations in the Heliosphere. Due to the intrinsic multiscale nature of features in coronagraph images, wavelet and multiscale image processing techniques are well suited to enhancing the visibility of CMEs and supressing noise. However, wavelets are better suited to identifying point-like features, such as noise or background stars, than to enhancing the visibility of the curved form of a typical CME front. Higher order multiscale techniques, such as ridgelets and curvelets, were therefore explored to characterise the morphology (width, curvature) and kinematics (position, velocity, acceleration) of CMEs. Curvelets in particular were found to be well suited to characterising CME properties in a self-consistent manner. Curvelets are thus likely to be of benefit to autonomous monitoring of CME properties for space weather applications.
[70]  [pdf] - 261770
Propagation of an Earth-directed coronal mass ejection in three dimensions
Comments: 5 figures, 2 supplementary movies
Submitted: 2010-10-04, last modified: 2010-11-12
Solar coronal mass ejections (CMEs) are the most significant drivers of adverse space weather at Earth, but the physics governing their propagation through the heliosphere is not well understood. While stereoscopic imaging of CMEs with the Solar Terrestrial Relations Observatory (STEREO) has provided some insight into their three-dimensional (3D) propagation, the mechanisms governing their evolution remain unclear due to difficulties in reconstructing their true 3D structure. Here we use a new elliptical tie-pointing technique to reconstruct a full CME front in 3D, enabling us to quantify its deflected trajectory from high latitudes along the ecliptic, and measure its increasing angular width and propagation from 2-46 solar radii (approximately 0.2 AU). Beyond 7 solar radii, we show that its motion is determined by an aerodynamic drag in the solar wind and, using our reconstruction as input for a 3D magnetohydrodynamic simulation, we determine an accurate arrival time at the Lagrangian L1 point near Earth.
[71]  [pdf] - 1032858
Large-scale Bright Fronts in the Solar Corona: A Review of "EIT waves"
Comments: 34 pages, 19 figures
Submitted: 2010-06-01, last modified: 2010-11-11
``EIT waves" are large-scale coronal bright fronts (CBFs) that were first observed in 195 \AA\ images obtained using the Extreme-ultraviolet Imaging Telescope (EIT) onboard the \emph{Solar and Heliospheric Observatory (SOHO)}. Commonly called ``EIT waves", CBFs typically appear as diffuse fronts that propagate pseudo-radially across the solar disk at velocities of 100--700 km s$^{-1}$ with front widths of 50-100 Mm. As their speed is greater than the quiet coronal sound speed ($c_s\leq$200 km s$^{-1}$) and comparable to the local Alfv\'{e}n speed ($v_A\leq$1000 km s$^{-1}$), they were initially interpreted as fast-mode magnetoacoustic waves ($v_{f}=(c_s^2 + v_A^2)^{1/2}$). Their propagation is now known to be modified by regions where the magnetosonic sound speed varies, such as active regions and coronal holes, but there is also evidence for stationary CBFs at coronal hole boundaries. The latter has led to the suggestion that they may be a manifestation of a processes such as Joule heating or magnetic reconnection, rather than a wave-related phenomena. While the general morphological and kinematic properties of CBFs and their association with coronal mass ejections have now been well described, there are many questions regarding their excitation and propagation. In particular, the theoretical interpretation of these enigmatic events as magnetohydrodynamic waves or due to changes in magnetic topology remains the topic of much debate.
[72]  [pdf] - 370695
Solar Wind Drag and the Kinematics of Interplanetary Coronal Mass Ejections
Submitted: 2010-10-01
Coronal mass ejections (CMEs) are large-scale ejections of plasma and magnetic field from the solar corona, which propagate through interplanetary space at velocities of $\sim$100--2500~km~s$^{-1}$. Although plane-of-sky coronagraph measurements have provided some insight into their kinematics near the Sun ($<$32~R$_\odot$), it is still unclear what forces govern their evolution during both their early acceleration and later propagation. Here, we use the dual perspectives of the Solar TErrestrial RElations Observatory (STEREO) spacecrafts to derive the three-dimensional kinematics of CMEs over a range of heliocentric distances ($\sim$2--250\,R$_{\odot}$). We find evidence for solar wind (SW) drag-forces acting in interplanetary space, with a fast CME decelerated and a slow CME accelerated towards typical SW velocities. We also find that the fast CME showed linear ($\delta=1$) dependence on the velocity difference between the CME and the SW, while the slow CME showed a quadratic ($\delta=2$) dependence. The differing forms of drag for the two CMEs indicate the forces and thus mechanism responsible for there acceleration may be different.
[73]  [pdf] - 197113
Solar Magnetic Feature Detection and Tracking for Space Weather Monitoring
Submitted: 2010-06-30
We present an automated system for detecting, tracking, and cataloging emerging active regions throughout their evolution and decay using SOHO Michelson Doppler Interferometer (MDI) magnetograms. The SolarMonitor Active Region Tracking (SMART) algorithm relies on consecutive image differencing to remove both quiet-Sun and transient magnetic features, and region-growing techniques to group flux concentrations into classifiable features. We determine magnetic properties such as region size, total flux, flux imbalance, flux emergence rate, Schrijver's R-value, R* (a modified version of R), and Falconer's measurement of non-potentiality. A persistence algorithm is used to associate developed active regions with emerging flux regions in previous measurements, and to track regions beyond the limb through multiple solar rotations. We find that the total number and area of magnetic regions on disk vary with the sunspot cycle. While sunspot numbers are a proxy to the solar magnetic field, SMART offers a direct diagnostic of the surface magnetic field and its variation over timescale of hours to years. SMART will form the basis of the active region extraction and tracking algorithm for the Heliophysics Integrated Observatory (HELIO).
[74]  [pdf] - 195856
Characterising Complexity in Solar Magnetogram Data using a Wavelet-based Segmentation Method
Comments: accepted for publication in ApJ
Submitted: 2010-05-10
The multifractal nature of solar photospheric magnetic structures are studied using the 2D wavelet transform modulus maxima (WTMM) method. This relies on computing partition functions from the wavelet transform skeleton defined by the WTMM method. This skeleton provides an adaptive space-scale partition of the fractal distribution under study, from which one can extract the multifractal singularity spectrum. We describe the implementation of a multiscale image processing segmentation procedure based on the partitioning of the WT skeleton which allows the disentangling of the information concerning the multifractal properties of active regions from the surrounding quiet-Sun field. The quiet Sun exhibits a average H\"older exponent $\sim -0.75$, with observed multifractal properties due to the supergranular structure. On the other hand, active region multifractal spectra exhibit an average H\"older exponent $\sim 0.38$ similar to those found when studying experimental data from turbulent flows.
[75]  [pdf] - 1026186
Constraining 3D Magnetic Field Extrapolations Using The Twin Perspectives of STEREO
Submitted: 2010-04-08
The 3D magnetic topology of a solar active region (NOAA 10956) was reconstructed using a linear force-free field extrapolation constrained using the twin perspectives of \emph{STEREO}. A set of coronal field configurations was initially generated from extrapolations of the photospheric magnetic field observed by the Michelson Doppler Imager (MDI) on \emph{SOHO}. Using an EUV intensity-based cost function, the extrapolated field lines that were most consistent with 171\AA\ passband images from the Extreme UltraViolet Imager (EUVI) on \emph{STEREO} were identified. This facilitated quantitative constraints to be placed on the twist ($\alpha$) of the extrapolated field lines, where $\nabla \times {\bf B} = \alpha {\bf B}$. Using the constrained values of $\alpha$, the evolution in time of twist, connectivity, and magnetic energy were then studied. A flux emergence event was found to result in significant changes in the magnetic topology and total magnetic energy of the region.
[76]  [pdf] - 1025954
Investigating the driving mechanisms of coronal mass ejections
Comments: 9 pages 7 figures
Submitted: 2010-03-26
The objective of this investigation was to first examine the kinematics of coronal mass ejections (CMEs) using EUV and coronagraph images, and then to make a comparison with theoretical models in the hope to identify the driving mechanisms of the CMEs. We have studied two CMEs which occurred on 2006 Dec. 17 (CME06) and 2007 Dec. 31 (CME07). The models studied in this work were catastrophe, breakout, and toroidal instability models. We found that after the eruption, the accelerations of both events exhibited a drop before increasing again. Our comparisons with the theories suggested that CME06 can be best described by a hybrid of the catastrophe and breakout models while CME07 is most consistent with the breakout model.
[77]  [pdf] - 901904
Turbulence, Complexity, and Solar Flares
Comments: 2 figures
Submitted: 2009-09-30
The issue of predicting solar flares is one of the most fundamental in physics, addressing issues of plasma physics, high-energy physics, and modelling of complex systems. It also poses societal consequences, with our ever-increasing need for accurate space weather forecasts. Solar flares arise naturally as a competition between an input (flux emergence and rearrangement) in the photosphere and an output (electrical current build up and resistive dissipation) in the corona. Although initially localised, this redistribution affects neighbouring regions and an avalanche occurs resulting in large scale eruptions of plasma, particles, and magnetic field. As flares are powered from the stressed field rooted in the photosphere, a study of the photospheric magnetic complexity can be used to both predict activity and understand the physics of the magnetic field. The magnetic energy spectrum and multifractal spectrum are highlighted as two possible approaches to this.
[78]  [pdf] - 24237
Reconstructing the 3-D Trajectories of CMEs in the Inner Heliosphere
Comments: 17 pages, 10 figures, 1 appendix
Submitted: 2009-05-13
A method for the full three-dimensional (3-D) reconstruction of the trajectories of coronal mass ejections (CMEs) using Solar TErrestrial RElations Observatory (STEREO) data is presented. Four CMEs that were simultaneously observed by the inner and outer coronagraphs (COR1 and 2) of the Ahead and Behind STEREO satellites were analysed. These observations were used to derive CME trajectories in 3-D out to ~15Rsun. The reconstructions using COR1/2 data support a radial propagation model. Assuming pseudo-radial propagation at large distances from the Sun (15-240Rsun), the CME positions were extrapolated into the Heliospheric Imager (HI) field-of-view. We estimated the CME velocities in the different fields-of-view. It was found that CMEs slower than the solar wind were accelerated, while CMEs faster than the solar wind were decelerated, with both tending to the solar wind velocity.
[79]  [pdf] - 24168
Automated Coronal Hole Detection using Local Intensity Thresholding Techniques
Submitted: 2009-05-12
We identify coronal holes using a histogram-based intensity thresholding technique and compare their properties to fast solar wind streams at three different points in the heliosphere. The thresholding technique was tested on EUV and X-ray images obtained using instruments onboard STEREO, SOHO and Hinode. The full-disk images were transformed into Lambert equal-area projection maps and partitioned into a series of overlapping sub-images from which local histograms were extracted. The histograms were used to determine the threshold for the low intensity regions, which were then classified as coronal holes or filaments using magnetograms from the SOHO/MDI. For all three instruments, the local thresholding algorithm was found to successfully determine coronal hole boundaries in a consistent manner. Coronal hole properties extracted using the segmentation algorithm were then compared with in situ measurements of the solar wind at 1 AU from ACE and STEREO. Our results indicate that flux tubes rooted in coronal holes expand super-radially within 1 AU and that larger (smaller) coronal holes result in longer (shorter) duration high-speed solar wind streams.
[80]  [pdf] - 20576
The kinematics of coronal mass ejections using multiscale methods
Comments: 10 pages, 13 figures, accepted for publication
Submitted: 2009-01-21
The diffuse morphology and transient nature of coronal mass ejections (CMEs) make them difficult to identify and track using traditional image processing techniques. We apply multiscale methods to enhance the visibility of the faint CME front. This enables an ellipse characterisation to objectively study the changing morphology and kinematics of a sample of events imaged by the Large Angle Spectrometric Coronagraph (LASCO) onboard the Solar and Heliospheric Observatory (SOHO) and the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) onboard the Solar Terrestrial Relations Observatory (STEREO). The accuracy of these methods allows us to test the CMEs for non-constant acceleration and expansion. We exploit the multiscale nature of CMEs to extract structure with a multiscale decomposition, akin to a Canny edge detector. Spatio-temporal filtering highlights the CME front as it propagates in time. We apply an ellipse parameterisation of the front to extract the kinematics (height, velocity, acceleration) and changing morphology (width, orientation). The kinematic evolution of the CMEs discussed in this paper have been shown to differ from existing catalogues. These catalogues are based upon running-difference techniques that can lead to over-estimating CME heights. Our resulting kinematic curves are not well-fitted with the constant acceleration model. It is shown that some events have high acceleration below $\sim$5 R$_{\sun}$. Furthermore, we find that the CME angular widths measured by these catalogues are over-estimated, and indeed for some events our analysis shows non-constant CME expansion across the plane-of-sky.
[81]  [pdf] - 19048
Multi-wavelength observations and modelling of a canonical solar flare
Comments: 10 pages, 7 figures, 2 tables. Accepted for publication in A&A
Submitted: 2008-12-01
This paper investigates the temporal evolution of temperature, emission measure, energy loss and velocity in a C-class solar flare from both an observational and theoretical perspective. The properties of the flare were derived by following the systematic cooling of the plasma through the response functions of a number of instruments -- RHESSI (>5 MK), GOES-12 (5-30 MK), TRACE 171 A (1 MK) and SOHO/CDS (~0.03-8 MK). These measurements were studied in combination with simulations from the 0-D EBTEL model. At the flare on-set, upflows of ~90 km s-1 and low level emission were observed in Fe XIX, consistent with pre-flare heating and gentle chromospheric evaporation. During the impulsive phase, upflows of ~80 km s-1 in Fe XIX and simultaneous downflows of 20 km s-1 in He I and O V were observed, indicating explosive chromospheric evaporation. The plasma was subsequently found to reach a peak temperature of ~13 MK in approximately 10 minutes. Using EBTEL, conduction was found to be the dominant loss mechanism during the initial ~300s of the decay phase. It was also found to be responsible for driving gentle chromospheric evaporation during this period. As the temperature fell below ~8 MK, and for the next ~4,000s, radiative losses were determined to dominate over conductive losses. The radiative loss phase was accompanied by significant downflows of <40 km s-1 in O V. This is the first extensive study of the evolution of a canonical solar flare using both spectroscopic and broad-band instruments in conjunction with a hydrodynamic model. While our results are in broad agreement with the standard flare model, the simulations suggest that both conductive and non-thermal beam heating play important roles in heating the flare plasma during the impulsive phase of at least this event.
[82]  [pdf] - 12612
The Kinematics of a Globally Propagating Disturbance in the Solar Corona
Comments: 4 pages 4 figures
Submitted: 2008-05-14
The kinematics of a globally propagating disturbance (also known as an ``EIT wave") is discussed using Extreme UltraViolet Imager (EUVI) data Solar Terrestrial Relations Observatory (STEREO). We show for the first time that an impulsively generated propagating disturbance has similar kinematics in all four EUVI passbands (304, 171, 195, and 284 A). In the 304 A passband the disturbance shows a velocity peak of 238+/-20 kms-1 within ~28 minutes of its launch, varying in acceleration from 76 ms-2 to -102 ms-2. This passband contains a strong contribution from a Si XI line (303.32 A) with a peak formation temperature of ~1.6 MK. The 304 A emission may therefore be coronal rather than chromospheric in origin. Comparable velocities and accelerations are found in the coronal 195 A passband, while lower values are found in the lower cadence 284 A passband. In the higher cadence 171 A passband the velocity varies significantly, peaking at 475+/-47 kms-1 within ~20 minutes of launch, with a variation in acceleration from 816 ms-2 to -413 ms-2. The high image cadence of the 171 A passband (2.5 minutes compared to 10 minutes for the similar temperature response 195 A passband) is found to have a major effect on the measured velocity and acceleration of the pulse, which increase by factors of ~2 and ~10, respectively. This implies that previously measured values (e.g., using EIT) may have been underestimated. We also note that the disturbance shows strong reflection from a coronal hole in both the 171 and 195 A passbands. The observations are consistent with an impulsively generated fast-mode magnetoacoustic wave.
[83]  [pdf] - 12241
Multiresolution analysis of active region magnetic structure and its correlation with the Mt. Wilson classification and flaring activity
Comments: Accepted for publication in Solar Physics
Submitted: 2008-05-01
Two different multi-resolution analyses are used to decompose the structure of active region magnetic flux into concentrations of different size scales. Lines separating these opposite polarity regions of flux at each size scale are found. These lines are used as a mask on a map of the magnetic field gradient to sample the local gradient between opposite polarity regions of given scale sizes. It is shown that the maximum, average and standard deviation of the magnetic flux gradient for alpha, beta, beta-gamma and beta-gamma-delta active regions increase in the order listed, and that the order is maintained over all length-scales. This study demonstrates that, on average, the Mt. Wilson classification encodes the notion of activity over all length-scales in the active region, and not just those length-scales at which the strongest flux gradients are found. Further, it is also shown that the average gradients in the field, and the average length-scale at which they occur, also increase in the same order. Finally, there are significant differences in the gradient distribution, between flaring and non-flaring active regions, which are maintained over all length-scales. It is also shown that the average gradient content of active regions that have large flares (GOES class 'M' and above) is larger than that for active regions containing flares of all flare sizes; this difference is also maintained at all length-scales.
[84]  [pdf] - 11710
Multiscale Edge Detection in the Corona
Comments: 14 pages, 7 figures, In Solar Physics Topical Issue "Solar Image Analysis and Visualization"
Submitted: 2008-04-11
Coronal Mass Ejections (CMEs) are challenging objects to detect using automated techniques, due to their high velocity and diffuse, irregular morphology. A necessary step to automating the detection process is to first remove the subjectivity introduced by the observer used in the current, standard, CME detection and tracking method. Here we describe and demonstrate a multiscale edge detection technique that addresses this step and could serve as one part of an automated CME detection system. This method provides a way to objectively define a CME front with associated error estimates. These fronts can then be used to extract CME morphology and kinematics. We apply this technique to a CME observed on 18 April 2000 by the Large Angle Solar COronagraph experiment (LASCO) C2/C3 and a CME observed on 21 April 2002 by LASCO C2/C3 and the Transition Region and Coronal Explorer (TRACE). For the two examples in this work, the heights determined by the standard manual method are larger than those determined with the multiscale method by approximately 10% using LASCO data and approximately 20% using TRACE data.
[85]  [pdf] - 80861
Observational Evidence for Gentle Chromospheric Evaporation During the Impulsive Phase of a Solar Flare
Comments: 5 pages
Submitted: 2006-03-23
Observational evidence for gentle chromospheric evaporation during the impulsive phase of a C9.1 solar flare is presented using data from the Reuven Ramaty High-Energy Solar Spectroscopic Imager and the Coronal Diagnostic Spectrometer on board the Solar and Heliospheric Observatory. Until now, evidence for gentle evaporation has often been reported during the decay phase of solar flares, where thermal conduction is thought to be the driving mechanism. Here we show that the chromospheric response to a low flux of nonthermal electrons (>=5x10^9 ergs cm^-2 s^-1) results in plasma upflows of 13+/-16, 16+/-18, and 110+/-58 km s^-1 in the cool He I and O V emission lines and the 8 MK Fe XIX line. These findings, in conjunction with other recently reported work, now confirm that the dynamic response of the solar atmosphere is sensitively dependent on the flux of incident electrons.
[86]  [pdf] - 76245
RHESSI and SOHO/CDS Observations of Explosive Chromospheric Evaporation
Comments: 5 pages, 4 figures, ApJL (In Press)
Submitted: 2005-09-22, last modified: 2005-12-21
Simultaneous observations of explosive chromospheric evaporation are presented using data from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and the Coronal Diagnostic Spectrometer (CDS) onboard SOHO. For the first time, co-spatial imaging and spectroscopy have been used to observe explosive evaporation within a hard X-ray emitting region. RHESSI X-ray images and spectra were used to determine the flux of non-thermal electrons accelerated during the impulsive phase of an M2.2 flare. Assuming a thick-target model, the injected electron spectrum was found to have a spectral index of ~7.3, a low energy cut-off of ~20 keV, and a resulting flux of >4x10^10 ergs cm^-2 s^-1. The dynamic response of the atmosphere was determined using CDS spectra, finding a mean upflow velocity of 230+/-38 km s^-1 in Fe XIX (592.23A), and associated downflows of 36+/-16 km s^-1 and 43+/-22 km s^-1 at chromospheric and transition region temperatures, respectively, relative to an averaged quiet-Sun spectra. The errors represent a 1 sigma dispersion. The properties of the accelerated electron spectrum and the corresponding evaporative velocities were found to be consistent with the predictions of theory.
[87]  [pdf] - 75801
Plasma Diagnostics of Active Region Evolution and Implications for Coronal Heating
Comments: 9 pages, 11 figures. MNRAS, In Press
Submitted: 2005-09-08
A detailed study is presented of the decaying solar active region NOAA 10103 observed with the Coronal Diagnostic Spectrometer (CDS), the Michelson Doppler Imager (MDI) and the Extreme-ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory (SOHO). Electron density maps formed using Si X (356.03A/347.41A) show that the density varies from ~10^10 cm^-3 in the active region core, to ~7x10^8 cm^-3 at the region boundaries. Over the five days of observations, the average electron density fell by ~30%. Temperature maps formed using Fe XVI(335.41A)/Fe XIV(334.18A) show electron temperatures of \~2.34x10^6 K in the active region core, and ~2.10x10^6 K at the region boundaries. Similarly to the electron density, there was a small decrease in the average electron temperature over the five day period. The radiative, conductive, and mass flow losses were calculated and used to determine the resultant heating rate (P_H). Radiative losses were found to dominate the active region cooling process. As the region decayed, the heating rate decreased by almost a factor of five between the first and last day of observations. The heating rate was then compared to the total unsigned magnetic flux (Phi_tot), yielding a power-law of the form P_H ~ Phi_tot^(0.81 +/- 0.32). This result suggests that waves rather than nanoflares may be the dominant heating mechanism in this active region.
[88]  [pdf] - 61114
Initial results from SECIS observations of the 2001 eclipse
Comments: 6 pages, 3 figures, oral contribution to the "SOHO 13: Waves, oscillations and small scale events in the Solar Atmosphere" conference, proceedings ESA SP-546
Submitted: 2003-11-24, last modified: 2003-11-26
SECIS observations of the June 2001 total solar eclipse were taken using an Fe xiv 5303 A filter. Existing software was modified and new code was developed for the reduction and analysis of these data. The observations, data reduction, study of the atmospheric and instrumental effects, together with some preliminary results are discussed. Emphasis is given to the techniques used for the automated alignment of the 8000 images, the application of the a Trous algorithm for noise filtering and the software developed for the automated detection of intensity oscillations using wavelet analysis. In line with findings from the 1999 SECIS total eclipse observations, intensity oscillations with periods in the range of 20-30 s, both inside and just outside coronal loops are also presented.
[89]  [pdf] - 56698
Eclipse observations of high-frequency oscillations in active region coronal loops
Comments: 6 pages, 6 figures
Submitted: 2003-05-13
One of the mechanisms proposed for heating the corona above solar active regions is the damping of magnetohydrodynamic (MHD) waves. Continuing on previous work, we provide observational evidence for the existence of high-frequency MHD waves in coronal loops observed during the August 1999 total solar eclipse. A wavelet analysis is used to identify twenty 4x4 arcsec2 areas showing intensity oscillations. All detections lie in the frequency Hz (5-3 s), last for at least 3 periods at a confidence level of more than 99% and arise just outside known coronal loops. This leads us to suggest that they occur in low emission-measure or different temperature loops associated with the active region.